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Abstract 



Qh! We employ non-perturbative flow equations to compute the equation of state for 



two flavor QCD within an effective quark meson model. This yields the temperature 
and quark mass dependence of quantities like the chiral condensate or the pion mass. 
Our treatment covers both the chiral perturbation theory domain of validity and the 
domain of validity of universality associated with critical phenomena: We explicitly 
connect the physics at zero temperature and realistic quark mass with the universal 
behavior near the critical temperature T c and the chiral limit. For realistic quark 
masses the pion correlation length near T c turns out to be smaller than its zero 
temperature value. In the vicinity of T c and zero quark mass we obtain a precision 
estimate of the universal critical equation of state. It belongs to the universality class 
of the three dimensional 0(4) symmetric Heisenberg model. We also investigate 
scalar field theories directly in three dimensions near second order and first order 
phase transitions. In particular, we obtain the universal equation of state for weak 
first order phase transitions in scalar matrix models. We also study the coarse grained 
free energy. Its dependence on the coarse graining scale gives a quantitative criterion 
for the validity of Langer's theory of spontaneous bubble nucleation. 
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1 Introduction and overview 

Quantum chromodynamics (QCD) describes qualitatively different physics at different 
length scales. The theory is asymptotically free [1] and at short distances or high en- 
ergies the strong interaction dynamics of quarks and gluons can be determined from 
perturbation theory. On the other hand, at scales of a few hundred MeV confinement 
sets in and the spectrum of the theory consists only of color neutral states. A further 
important aspect of strong interaction dynamics is spontaneous chiral symmetry break- 
ing. In the limit of vanishing current quark masses the classical or "short distance" 
QCD action does not couple left- and right-handed quarks. As a consequence, for Nf 
massless quark flavors the classical QCD action exhibits a global chiral invariance under 
U L (N f ) x U R (N f ) = SU L (N f ) x SU R (Nf) x U v (l) x U A (1), with different unitary trans- 
formations acting on the left- and right-handed quark fields. However, this symmetry is 
not fully observed in the hadron spectrum. It is phenomenologically well established that 
in the quantum theory the symmetry SU^Nf) x SUu(Nf) is spontaneously broken to the 
diagonal SUv{Nf) vector-like subgroup. For Nf = 2 the isospin symmetry is realized in 
the hadron spectrum to a very good approximation. The comparably stronger explicit 
symmetry breaking due to non-zero current quark masses makes the realization of the 
SUv{3) symmetry less accurate. In addition, the axial abelian subgroup Ua(1) is broken 
in the quantum theory by an anomaly of the axial-vector current [2]. (The abelian Uy(l) 
subgroup corresponds to baryon number conservation.) 

Extrapolating QCD from short distance to long distance scales is clearly a non-pertur- 
bative problem. The strong gauge coupling a s changes with scale, growing with increased 
distance. However, not only effective couplings but also the relevant degrees of freedom 
can change with scale. Indeed, at low energies an essential part of strong interaction 
dynamics can be encoded in the masses and interactions of mesons. Prominent examples 
of systematic effective descriptions are the successes of chiral perturbation theory [3, 4] 
based on the non-linear sigma model or a description within the linear sigma model [5] 1 . 
The success of any analytical approach to bridge the gap between short distance and long 
distance QCD may depend on the ability to introduce field variables for composite objects 
such as mesons at a certain "compositeness" scale (cf. section 5). 

An even more challenging subject is finite temperature QCD [8]. In view of present and 
future heavy-ion collision experiments [9], but also for its cosmological relevance [10], there 
is an increased interest in the thermal properties of QCD. It was realized early [11] that in 
a thermal equilibrium situation at sufficiently high temperature or density QCD may differ 
in important aspects from the corresponding zero temperature or vacuum properties. If 
one considers the phase structure of QCD as a function of temperature one expects at a 
critical temperature around T c ~ 100 — 200 MeV two types of phase transitions or crossover 
phenomena 2 . First, there is the expectation that the confined quarks and gluons at zero 

1 It has been demonstrated recently [6, 7] that the results of chiral perturbation theory can be reproduced 
within the linear meson model once certain higher dimensional operators in its effective action are taken 
into account. 

transitions are expected for both increasing density and increasing temperature [12, 13]. We will 
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temperature become deconfined at sufficiently high temperature where the hot state is 
conventionally called the quark-gluon plasma. In the pure gauge theory, or equivalently 
in the limit of infinite quark masses, a high temperature deconfinement phase transition 
is well-established [14, 15]. Second, the spontaneously broken chiral symmetry in the 
limit of vanishing quark masses and zero temperature is expected to be restored in a high 
temperature phase. It is not yet clear whether these phase transitions occur at the same 
temperature or whether they persist for realistic quark masses or rather turn into smooth 
crossover. 

Such a transition must have occured during the early stages of the evolution of the 
universe. According to the standard hot big bang cosmology the transition took place 
about a microsecond after the big bang, where the temperature dropped to the order of 
100 MeV. Three orders of magnitude higher the previous transition occured, related to the 
dynamics of the electroweak interactions 3 , and for most of its evolution the early universe 
was to a good approximation in thermal equilibrium [10]. 

A promising prospect is to reproduce the QCD transition in heavy-ion collisions in the 
laboratory at sufficiently high energy. Present fixed target heavy-ion collision experiments 
at the AGS accelerator (BNL) and at the CERN SPS accelerator have produced heavy- 
ion beams with energies that could yield initial temperatures above T c [20]. Here the 
initial temperature refers to a frequently used scenario in which the central rapidity region 
reaches a stage of local thermal equilibrium, preceded by an intermediate pre-equilibrium 
stage after the first instant of nuclear contact [21, 22]. In this scenario the thermalized 4 
central region is characterized by a large number ratio of mesons to baryons [22], thus 
QCD at finite temperature and zero baryon number density would be of relevance for its 
description. Whereas SPS realizes for very heavy ion Pb beams a center-of-mass energy of 
160 GeV per nucleon, future colliders, RHIC with 200 GeV per nucleon (gold on gold) and 
the LHC with 6300 GeV per nucleon (lead on lead), will allow to produce temperatures 
substantially larger that the typical time scales of 10~ 23 sec. or distances of 1 fm associated 
with the QCD transition [25]. 

Of particular interest is the order of a possible phase transition which may have far- 
reaching phenomenological consequences. One reason is that a first order phase transition 
may create an out of equilibrium situation. For a first order transition the (coarse grained) 
free energy as a function of a suitable order parameter exhibits two distinct minima near 
the critical temperature T c . If in a heavy-ion collision of sufficiently high energy a region 
of the high temperature phase is created, it then cools through the phase transition. The 
initial high temperature phase becomes a local minimum of the free energy and the lowest 
minimum corresponds to the low temperature phase. However, a barrier separating both 

restrict the discussion here to QCD at non-zero temperature. 

3 Existence and properties of the electroweak phase transition depend strongly on the mass of the Higgs 
scalar Mh ■ For Higgs boson masses around Mh ^ 80 GeV there is no phase transition but rather a smooth 
crossover [16, 17, 18, 19]. 

4 Thermalization enters as an assumption that remains to be justified. There are results from the parton 
cascade model [23, 24] which give some support to the idea that there will be enough time in a RHIC 
collision for local thermal equilibrium to be established. 
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minima prevents a smooth transition at T c . At a first order phase transition, typically 
the conversion proceeds in a short dramatic period through the nucleation and growth of 
"bubbles" (droplets) of the low temperature phase (cf. section 4). 

In contrast to the discontinuous behavior at a first order transition, at a second order 
phase transition the order parameter changes continuously (but non-analytically) from one 
phase to the other. In particular, there is no mass scale present at a second order phase 
transition which corresponds to infinite correlation lengths. Large correlation lengths may 
yield distinctive, qualitative signatures in a relativistic heavy-ion collision. It has been 
argued [26, 27] that in the central rapidity region a large correlation length may be respon- 
sible for the creation of large domains, in which the pion field has a non-zero expectation 
value pointing in a fixed direction in isospin space ("disoriented chiral condensate" [28]). 
The observational consequences would be strong fluctuations in the number ratio of neu- 
tral to charged pions. Such a phenomenon was even related [26, 27] to the prominent 
Centauro events [29] in high energy cosmic ray experiments. Apart from the phenomeno- 
logical implications, a large correlation length opens the spectacular possibility that the 
QCD phase transition is characterized by universal properties! The notion of universality 
for critical phenomena is well-established in statistical physics [30]. Universal properties 
are independent of the details, like short distance couplings, of the model under investiga- 
tion. They only depend on the symmetries, the dimensionality of space and the number 
of field degrees of freedom. As a consequence a whole class of models is described by the 
same universal form of the equation of state in the vicinity of the critical temperature. The 
range of applicability typically covers very different physical systems in condensed matter 
physics and high temperature quantum field theory (cf. section 3). 

Finally there is the possibility of no true phase transition at all. For such a crossover 
phenomenon quantities may change rapidly in the transition region, and the phenomeno- 
logical consequences may be almost indistinguishable from a situation with a true phase 
transition. However, for a crossover there occurs no non-analytical behavior and all cor- 
relation lengths remain finite. 

Large correlation lengths near a second order (or weak first order) phase transition are 
appealing both from the theoretical and the experimental point of view. For the experi- 
mentalist they may yield a unique, distinctive signature that a QCD phase transition has 
occured. The theoretical challenge is to find a method that describes both the known low- 
energy zero temperature properties of QCD and the possible universal behavior associated 
with a large correlation length near the critical temperature of a second order phase transi- 
tion. The link between the zero temperature and the universal critical properties is crucial. 
One may think of separating the problem into two parts: First, one considers a tractable, 
easier model than the original one where both models belong to the same universality class. 
The simplified description can develop enormous predictive power because of the universal 
properties of the equation of state near the critical temperature. For instance, a number of 
stimulating contributions [31, 26, 27] pointed out that for sufficiently small up and down 
quark masses and sufficiently large strange quark mass the chiral phase transition is ex- 
pected to be in the universality class of the three dimensional 0(4) symmetric Heisenberg 
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model 5 . This model is known to exhibit a second order phase transition, where universal 
properties can be observed. Apart from the expected observational consequences, a num- 
ber of interesting conclusions have been drawn from universality arguments [26, 27], e.g. 
that the pion mass increases with temperature around T c . In particular, near the critical 
temperature only the pions and the sigma resonance play a role for the behavior of the 
chiral condensate and long-distance correlation functions. Here the chiral condensate, i.e. 
the expectation value of the quark bilinear (ipip), plays the role of an order parameter for 
the chiral transition in the limit of vanishing current quark masses. However, the questions 
how small the up and down quark mass m u and m d would have to be in order to obtain a 
large correlation length near T c and if the drawn conclusions are valid for realistic quark 
mass values remained unanswered. Nothing can be infered from universality arguments 
about the absolute magnitude of correlation lengths. Likewise the question of the value 
of the critical temperature is non-universal in nature. The non-universal amplitudes cru- 
cially depend on the zero temperature properties of the theory under investigation. In a 
second step one therefore finally has to link the universal behavior near T c and zero current 
quark mass on one hand with the known physical properties at T = for realistic quark 
masses on the other hand. 

One reason for the "missing link" is the breakdown of chiral perturbation theory for 
temperatures exceeding about one third of the zero temperature pion mass m n = 135 MeV 
(cf. section 5). Exploring the universal region in lattice QCD, on the other hand, is limited 
by present computer resources, though the use of improved actions seems promising [34]. 
Since decreasing quark masses increase the computational cost, simulations are typically 
done with unphysically heavy quarks where extrapolation to realistic values or to the chiral 
limit is needed. 

Equation of state for two flavor QCD 

It is a major result of this work to provide the equation of state for two flavor QCD 
which covers both the chiral perturbation theory domain of validity and the domain of 
validity of universality associated with critical phenomena [33]. The equation of state 
connects in a temperature range below ~ 170 MeV the derivative of the free energy or 
the effective potential with the average light current quark mass m = (m u + m<i)/2. It 
therefore allows to study the temperature and quark mass dependence of various quantities, 
like the pion mass or the chiral condensate (ipip) or the pion decay constant f n etc. 
In addition, we obtain a precision estimate of the universal critical equation of state of 
the three dimensional 0(4) symmetric Heisenberg model in the vicinity of the critical 
temperature and the chiral limit. From the scaling form of the equation of state we extract 
universal critical exponents, amplitude ratios and couplings. Recent lattice simulations for 
the universal critical equation of state of the three dimensional 0(4) model show agreement 
with our results at the few per cent level (cf. section 5). In addition, we find all non- 

5 We note that for two flavors the relevant chiral symmetry SU{2) x SU{2) is locally isomorphic to 
0(4). This does not hold in the case of a speculative "effective restoration" [31, 32] of the axial Ua{1) at 
temperatures around T c . 
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universal information to be encoded in two independent amplitudes that are properly 
computed for the two flavor case. 

Though these results will be presented in the text in more detail we point out two 
important answers one obtains from this study: First of all, for a thermal equilibrium 
situation the chiral transition gives no indication for strong fluctuations of pions with long 
wavelength 6 . In the presence of light quark masses the second order phase transition is 
turned into a smooth crossover: The expectation value of the quark bilinear (ipip) de- 
creases smoothly as the temperature is increased. The longest correlation lengths near the 
"crossover temperature" T pC) which turns out to be around 130 MeV for realistic quark 
masses, even becomes smaller than at T = 0. (It should be emphasized, however, that 
a tricritical behavior with a massless excitation remains possible for three flavors. This 
would not be characterized by the universal behavior of the 0(4) model; see the discus- 
sion below.) The second answer concerns the applicability of universal, i.e. almost model 
independent, arguments for a description of the chiral phase transition. Though the quark 
masses are fixed in nature, varying the average current quark mass to zero allows us to 
observe critical behavior independent from the realistic physical situation. In this way one 
can probe to what extent universality approximately is realized for realistic or even heavier 
quark masses. Universality has been assumed e.g. in two flavor lattice QCD simulations to 
guide extrapolation to realistic quark mass values or even to the chiral limit [36]. Despite 
the observed comparably short correlation lengths at non-zero temperature, we find the 
approximate validity of the 0(4) scaling behavior over a large temperature interval near 
and above T c = 100.7 MeV even for quark masses slightly larger than the realistic ones (cf. 
section 5). 

Our approach is based on the use of an exact non-perturbative flow equation for a 
scale dependent effective action T k [37, 38], which is the generating functional of the IP I 
Green functions in the presence of an infrared cutoff ~ k. In a thermal equilibrium context 
the so-called effective average action 1^ becomes temperature dependent and describes a 
coarse grained free energy with a coarse graining length scale ~ fc -1 . The coarse grained 
description allows us to consider the relevant physics at a given momentum-like scale k. 
The approach is closely connected to the Wilsonian renormalization group [30, 39, 40, 41, 
42, 43], often also called exact renormalization group. An introduction to the method is 
given in section 2. 

The flow equation for 1^ can be supplemented by an exact formalism for the introduc- 
tion of composite field variables or, more generally, a change of degrees of freedom at some 
given scale k [44] . We employ for scales below a "compositeness scale" of k$ ~ 600 MeV a 
description in terms of quark and scalar mesonic degrees of freedom. This effective quark 
meson model can be obtained from QCD in principle by "integrating out" the gluon degrees 
of freedom and by introducing fields for composite operators [44, 45] (cf. the discussion 
in sections 5 and 6). In this picture the scale kq, is associated to the scale at which the 

6 We have restricted the discussion to the thermal equilibrium properties of QCD. The question of 
strong fluctuations of pions with long wavelength has been investigated in the context of a chiral QCD 
phase transition far from equilibrium [26, 35]. 
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formation of mesonic bound states can be observed [44] in the flow of the momentum de- 
pendent four-quark interaction. Though considerable progress has been made, a reliable 
quantitative derivation of the effective quark meson model from QCD is still missing. We 
emphasize, however, that the quantitative aspects of this derivation will be of minor rele- 
vance for our practical calculations in the mesonic sector: If the effective Yukawa coupling 
between the quarks and the mesons turns out to be strong at the compositeness scale we 
observe a fast approach of the scale dependent effective couplings to approximate partial 
infrared fixed points [46, 33]. This behavior is discussed in section 5.4. As a consequence, 
the detailed form of the meson potential at /c$ becomes unimportant, except for the value 
of one relevant scalar mass term . In this work we fix from phenomenological input 
such that f n = 92.4 MeV (for m n = 135 MeV) which sets our unit of mass for two flavor 
QCD. The only other input parameter we use is the constituent quark mass M q to deter- 
mine the scale k^. We consider a range 300 MeV ^ M q ^ 350 MeV and find a rather weak 
dependence of our results on the precise value of M q . We point out that though a strong 
Yukawa coupling at is phenomenologically suggested by the comparably large value of 
the constituent quark mass M q it enters our description as a (consistent) assumption. On 
the present level of approximation the quark meson model for two flavors can be reduced 
to the well known 0(4) symmetric Gell-Mann-Levy linear sigma model [5] for the three 
pions and the "sigma resonance", however, it is coupled to quarks now. The meson de- 
grees of freedom can be understood as quark-antiquark bound states. Spontaneous chiral 
symmetry breaking is described by a non-vanishing expectation value of the meson field 
in absence of quark masses. The non-perturbative flow equations for the quark-meson 
model are obtained from a suitable truncation of the exact flow equation for IV We solve 
them numerically. Further details about the model and the applied approximations can be 
found in section 5. 

Quark mass dependent phase structure 

The details of a possible QCD phase transition, like the order or the "strength" of 
the transition, crucially depend on the number of flavors and the current quark masses. 
Though the number of flavors and the quark masses are fixed in nature, on the theoretical 
side it is possible to adjust them in order to get more insight into the phase structure 
of QCD. A considerable amount of information about the phase structure of non-zero 
temperature QCD has been provided by lattice Monte Carlo simulations whose results will 
be instructive in the following. The majority of lattice simulations incorporating dynamical 
fermions have been done with two or four flavors of equal mass quarks. There are a few 
simulations carried out for two light and one heavier quark which comes closest to the 
physical light quark mass ratios. To save computational cost, however, simulations are 
typically done with unphysically heavy quarks. 

Lattice QCD simulations in the limit of two light quark flavors show no signs of a 
first order deconfinement phase transition [47, 36] . The expectation value of the Polyakov 
loop 7 increases smoothly with increased temperature when dynamical fermions are present. 

7 We note that the Polyakov loop represents no order parameter for a possible deconfinement phase 
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There is no sharply defined transition, but rather a smooth crossover from one regime to 
the other. In accordance with our results, the chiral transition for two small, non-zero 
quark masses is found to be described by a crossover phenomenon [47, 36]. The lattice 
results indicate that both crossover, associated to the deconfinement and to the chiral 
transition, occur at the same temperature approximately around 8 140 — 150 MeV. This 
qualitative picture almost remains the same for the more realistic case of two light and 
one heavier quark. Simulations with "2+1" flavors carried out by the Columbia group in 
the staggered fermion scheme find a first order signal from the chiral condensate for three 
light flavors. However, the first order transition disappears for an increased strange quark 
mass, even before the strange quark adopts its physical mass value [48]. In contrast, from a 
Wilson fermion approach the Tsukuba group finds results that indicate a first order chiral 
phase transition for realistic quark masses[49]. However, simulations with Wilson fermions 
in this context are difficult because they break chiral symmetry completely away from the 
continuum limit. Though no final conclusion can be drawn from the lattice results by now, 
if there is a first order phase transition at all it is very likely to be weak. In addition, 
lattice data seem to indicate that QCD exhibits a (smooth) transition which is dominated 
by the approximate restoration of chiral symmetry. 

Our results and the results from lattice simulations are compatible with the following 
picture proposed from effective sigma model considerations [31, 50, 26, 51]: For three mass- 
less quark flavors the chiral phase transition is first order in nature. If the strange quark 
mass m s is raised the phase space shows a line of first order transitions which ends in a 
tricritical point with a massless excitation. For larger strange quark masses the transition 
remains of second order. If one takes into account that the up and down quark masses 
are non-zero this second order phase transition is smoothed into a crossover phenomenon. 
However, there remains a line of tricritical points separating the first order region for small 
m s and the crossover region for large m s . It has been pointed out [51, 52] that the QCD 
chiral phase transition may be close to a tricritical point with a massless excitation. 

Effective three dimensional behavior at high temperature 

Our investigations have ruled out the appearance of a large correlation length at the 
equilibrium chiral transition in the two flavor case, corresponding to an infinite mass limit 
for the remaining quark flavors. However, we have seen that a chiral transition with a 
large correlation length remains possible through the influence of the strange quark. In 
particular, along the mentioned line of tricritical points the transition is expected to be in 
the universality class of the Ising model [51]. In this has been discussed above, 

in a first step one may learn as much as possible about the phase transition only rely- 
ing on universality arguments. This alone typically poses a complicated non-perturbative 
problem. The problems are related to the effectively three dimensional behavior of high 

transition in the presence of light quarks. Only in the limit of infinite quark masses, or equivalcntly in the 
pure gauge theory, it is an order parameter for the spontaneous breaking of the center Z(N) of the gauge 
group SU(N) at high temperature [14]. 

8 Here the zero temperature p meson mass is used to set the energy scale. 
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temperature quantum field theory in four dimensions. Quantum field theory at non-zero 
temperature T can be formulated in terms of an Euclidean functional integral where the 
"time" dimension is compactified on a torus with radius T^ 1 [53]. Accordingly, non-zero 
temperature results in (anti-) periodic boundary conditions for (fermionic) bosonic fields 
in the Euclidean "time" direction with periodicity T~ x . As a consequence, the zeroth com- 
ponent of the Euclidean four-momentum becomes discrete, i.e. go = 2/7rT for bosons and 
go = (2/ + l)nT for fermions with integers 9 I. It is suggestive that if the relevant corre- 
lation length becomes much larger than the inverse temperature the compactified "time" 
dimension cannot be resolved anymore. This phenomenon is known as "dimensional reduc- 
tion" [54]. The effective three dimensional behavior of the theory can also be understood 
from the observation that for sufficiently high temperature the theory is dominated by 
classical statistical fluctuations. In particular, the critical exponents which describe the 
singular behavior of various quantities near a second order phase transition are those of 
the corresponding classical system 10 [55] (cf. section 5). An important particularity of the 
dimensionally reduced theory is the absence of fermionic degrees of freedom. This can be 
seen as a consequence of the fact that the "lightest" fermionic excitation is of the order 
ttT, whereas there exists a zero mode (I = 0) for a bosonic field. We have explicitly verified 
the phenomenon of dimensional reduction and the absence of fermions within the effective 
quark meson model at high temperature (cf. section 5). 

Near a second order (or a sufficiently weak first order) phase transition we therefore 
have to deal with a purely bosonic three dimensional theory. Standard high temperature 
perturbation theory can give a reliable description only if the relevant dimensionless cou- 
plings remain small near the transition. In three dimensions the (effective) scalar quartic 
coupling A has positive canonical mass dimension. If one denotes by tur the relevant mass 
at the critical temperature, standard perturbation theory corresponds to an expansion 
in the dimensionless effective coupling X/m R . For a given fixed coupling A a perturba- 
tive expansion will only converge if mj is sufficiently large, itlr ^> A. In consequence, 
standard perturbation theory is not applicable to second order {tur = 0) or weak first 
order {tur <C A) phase transitions due to the apparent infrared divergences in a perturba- 
tive treatment. Obviously in a suitable non-perturbative treatment a "running" effective 
renormalized coupling Xr has to be introduced. The running coupling approaches zero at 
a second order phase transition, with Xr/ ttir going to a constant as the correlation length 
m^ 1 tends to infinity. Typically no small effective coupling characterizes the interactions 
at the phase transition (cf. section 3). 

The three dimensional systems may be treated within the e-expansion [56] or by the 
use of perturbation series at fixed dimension [57]. Most of the lattice results come from 
the analysis of high temperature expansions [58] or from Monte Carlo simulations [59]. 
For the critical equation of state of the three dimensional 0(4) model a result from the 
e-expansion up to second order in e = 4 — d (d = 3) is available [60]. The reliability of this 

9 These are commonly called Matsubara frequencies. 

10 In the vicinity of the critical temperature of a second order phase transition the correlation length of 
the relevant bosonic field is £ ~ T~ X {\T — T c \/T c )~ 1 ' with a (positive) exponent u, and one always faces a 
high temperature situation for T T c , i.e. £ 3> T -1 . 
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expansion is not guaranteed a priori and results from a lattice Monte Carlo simulation [61] 
that has recently been obtained show substantial deviations from the e-expansion results 
(cf. section 5). The lattice simulations are quite demanding in the vicinity of the critical 
temperature of a second order (or weak first order) phase transition and require careful 
correction for finite size effects. As has been pointed out above our results agree quite well 
with the lattice study. 

The use of non-perturbative flow equations within the framework of the effective aver- 
age action is particularly suitable for dealing with theories which are plagued by infrared 
problems in perturbation theory. The employed formulation is both infrared and ultraviolet 
finite in arbitrary dimensions. In particular, the method is not restricted to an expansion 
in small couplings which has been crucial for the treatment of the effective quark meson 
model. 

Equation of state for second order and (weak) first order transitions 

The first two parts of this work are devoted to the study of scalar field theories near the 
critical temperature of second order and first order phase transitions. These investigations 
are performed directly in three dimensions. We obtain a detailed quantitative picture of 
the critical equation of state in the vicinity of the second order phase transition of the 
O(N) symmetric iV-component model [62, 63, 64]. Similar to the case N — 4, for N — 1 
our results have been confirmed by lattice simulations [65, 66] while there remains a consid- 
erable discrepancy to the e-expansion results 11 . We point out that our approximations are 
based on a derivative expansion for the effective average action that takes into account the 
most general non-derivative, i.e. potential, term consistent with the symmetries. It does 
not rely on any restriction of the potential to a polynomial form. It therefore allows to 
study the complete non-analytic behavior of the effective potential, or equivalently the free 
energy, near the critical temperature of the second order phase transition. For N = A our 
universal results in three dimensions are found to correspond to the results from the quark 
meson model near the critical temperature and the chiral limit. This correspondence is a 
manifestation of the phenomenon of dimensional reduction. The O(N) model has a wide 
range of different applications in statistical physics and high temperature quantum field 
theory. For N = A the model also describes the scalar sector of the electroweak standard 
model in the limit of vanishing gauge and Yukawa couplings. In condensed matter physics 
N = 3 corresponds to the well known Heisenberg model used to describe the ferromag- 
netic phase transition. There are other applications like the helium superfluid transition 
(N — 2), liquid- vapor transition (N = 1) or statistical properties of long polymer chains 
(N = 0) [67]. The critical properties of these systems are presented in section 3. 

The second part of this work provides a detailed investigation of first order phase 
transitions. In particular weak first order phase transitions are a much less thoroughly 
studied subject in the literature than second order transitions. This becomes especially 
desirable in the context of a possible weak first order chiral phase transition in QCD. We 
consider [68, 69] models with U (N) x U (N) symmetry with a scalar field in the (JV, N) 

11 For N — 1 the e-expansion for the scaling equation of state is available up to order e 3 [67]. 
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representation, described by an arbitrary complex NxN matrix. We compute the equation 
of state for N = 2 which can be used for an effective description of two flavor QCD in 
the case of a speculative "effective restoration" [31, 32] of the chiral anomaly at high 
temperature. The considered theory also has a relation to the electroweak phase transition 
in models with two Higgs doublets. In certain limiting cases the model describes a non- 
linear matrix model for unitary matrices or one for singular 2x2 matrices (cf. section 4). 

For a large part of the parameter space the U(2)xU (2) symmetric matrix model exhibits 
a weak first order phase transition. Though the corresponding classical or short distance 
action indicates a second order phase transition, the transition becomes first order once 
fluctuations are taken into account. The fluctuation induced first order phase transition is 
known in four dimensions as the Coleman- Weinberg phenomenon [70]. A reliable method 
for a description of this phenomenon in three dimensions may be of relevance for the QCD 
chiral phase transition. For the realistic case of three flavors it has been argued that the 
chiral transition can be mainly fluctuation induced [51, 52]. 

We put special emphasis on the investigation of (approximate) universal properties 
near weak first order phase transitions. Whereas for second order transitions the universal 
equation of state can be expressed as a function of only one scaling variable (cf. section 3), 
due to an additional mass scale at a first order transition there is a second dimensionless 
ratio. The equation of state or, equivalently, the effective potential or free energy therefore 
depends on two independent scaling variables in the universal region. We have succeeded 
to describe this situation [68] for the present matrix models and present the universal 
equation of state in section 4.7. 

Coarse graining and first order transitions 

Another challenging problem is the discussion of the dynamics of a first order phase 
transition [71] which usually relies on the study of a non-convex potential or free energy. 
The decay of metastable states is associated either with tunneling fluctuations through 
barriers in the potential [72], or, at non-zero temperature, with thermal fluctuations above 
them [73]. However, the effective potential [75], which seems at first sight a natural tool 
for such studies, is expected to be a convex quantity with no barrier. The resolution of this 
paradox lies in the realization that the effective potential is convex because the tunneling 
or thermal fluctuations are incorporated in it. These fluctuations are associated with low 
frequency modes, while the non-convex part of the potential is related to the classical 
potential and the integration of high frequency modes. A natural approach to the study of 
first order phase transitions separates the problem in two parts. First, the high frequency 
modes are integrated out, with the possible generation of new minima of the potential 
through radiative symmetry breaking [70]. Subsequently, the decay of metastable states 
is discussed with semiclassical techniques [72, 73], using the non-convex potential that has 
resulted from the first step. This leads one to the notion of the coarse grained free energy, 
which is fundamental in statistical physics. Every physical system has a characteristic 
length scale associated with it. The dynamics of smaller length scales is integrated out, 
and is incorporated in the parameters of the free energy one uses for the study of the 
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behavior at larger length scales. Here we note that by construction, i.e. the inclusion 
of fluctuations with characteristic momenta larger than a given infrared cutoff ~ k, the 
effective average action 1^ is the appropriate quantity for the study of physics at a scale k 
(cf. section 2). It therefore realizes the concept of a coarse grained free energy in the sense 
of L anger [71]. 

The use of the effective average action will allow us to address the old question of the va- 
lidity [68, 69] of the standard theory of spontaneous bubble nucleation proposed by Langer 
[71] . This requires first a meaningful definition of a coarse grained free energy with a coarse 
graining scale k and second the validity of a saddle point approximation for the treatment 
of fluctuations around the "critical bubble" (cf. section 4.8). Here only fluctuations with 
momenta smaller than k must be included. Both issues turn out to be closely related. The 
validity of the saddle point approximation typically requires small dimensionless couplings. 
On the other hand we observe for large effective couplings that the form of the relevant 
coarse grained effective potential Uu depends strongly on the coarse graining scale k. This 
means that the lowest order in the saddle point approximation (classical contribution) de- 
pends strongly on the details of the coarse graining procedure. Since the final results as 
nucleation rates etc. must be independent of the coarse graining prescription this is only 
compatible with a large contribution from the higher orders of the saddle point expansion. 
We will consider this issue in section 4.8 in a quantitative way. 

This work is organized as follows. We provide an introduction into the non-perturbative 
effective average action method in section 2. The two following sections concern our studies 
of second order and first order phase transitions in three dimensional scalar field theories. 
These parts are based on refs. [62, 63, 64] and [68, 69] respectively. Section 5 is devoted 
to the chiral phase transition for two flavor QCD within the effective quark meson model. 
This part is based on refs. [33, 74]. Each section is supposed to be rather self-contained 
with a short introduction and conclusions to allow for a selective approach to the presented 
material. Section 6 presents an outlook. The acknowledgments are found in section 7. The 
three appendices A-C contain a discussion and definitions of threshold functions used in 
preceding sections. In appendix D the connection of the average current quark mass with 
the effective source term of the quark meson model is established. This is used in section 5. 
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2 Non-perturbative flow equations 

2.1 Effective average action 

We will introduce here the effective average action [37] r*. with an infrared cutoff ~ k. The 
effective average action is based on the quantum field theoretical concept of the effective 
action [75] T, i.e. the generating functional of the 1PI Green functions. The field equations 
derived from T include all quantum effects. For a field theoretical description of thermal 
equilibrium this concept is easily generalized to a temperature dependent effective action 
which includes in addition the thermal fluctuations. In statistical physics T describes the 
free energy as a functional of some (space dependent) order parameter. In particular, for 
a constant order parameter, T yields the effective potential which encodes the equation of 
state. The effective average action T k is a simple generalization of the effective action, with 
the distinction that only fluctuations with momenta q 2 <: k 2 are included. In the language 
of statistical physics, is a coarse grained free energy with a coarse graining scale k. 
Lowering k results in a successive inclusion of fluctuations with momenta q 2 ^ k 2 and 
therefore permits to explore the theory on larger and larger length scales. We note that Tk 
is closely related to an effective action for averages of fields [37], where the average is taken 
over a volume of size ~ k~ d and which is similar in spirit to the block-spin action [76, 30] 
in lattice theories. In a theory with a physical ultraviolet cutoff A we can associate r a with 
the microscopic or classical action S since no fluctuations below A are effectively included. 
For example, in statistical systems on a lattice the scale A plays the role of the inverse 
lattice spacing. In the context of dimensional reduction in high temperature quantum field 
theory A appears as the ultraviolet cutoff scale for the effective three dimensional models 
(cf. section 5.5). By definition, the effective average action equals the effective action for 
k = 0, i.e. To = T, since the infrared cutoff is absent. Thus interpolates between the 
classical action S and the effective action r as k is lowered from A to zero. The ability 
to follow the evolution to k — * is equivalent to the ability to solve the theory. The 
dependence of the effective average action on the scale k is described by an exact flow 
equation which is presented in section 2.2. 

To be explicit we consider the path integral representation of the generating functional 
for the connected Green functions in d Euclidean dimensions with N real scalar fields Xa 
(a = 1 . . . N), classical action S and sources J a . We introduce a /c-dependent generating 
functional 

W k [J] = In j D X exp - A k S[ X ] + / d d xJ a {x) X a {x)^ , (2.1) 

with an additional infrared cutoff term which is quadratic in the fields and reads in mo- 
mentum space 12 

AkS[x] = \ I ^M^ai-^iq)- (2.2) 

12 In order to avoid a proliferation of symbols we distinguish functions or operators in position 
space and their Fourier transforms according to their arguments, with the convention Xa(q) = 
J d d x cxp(-iq^)xa{x) . 
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Without this term W k equals the usual generating functional for the connected Green 
functions. Here the infrared cutoff function R k is required to vanish for k — > and to 
diverge for k — > A and fixed q 2 . For A — > oo this can be achieved, for example, by the 
choice 

Zi.a 2 e- q2/k2 

Mi) = (2 - 3) 

Here Z k denotes an appropriate wave function renormalization constant which will be 
defined later. For fluctuations with small momenta q 2 <C k 2 the cutoff R k — Z k k 2 acts like 
an additional mass term. For q 2 3> k 2 the infrared cutoff vanishes such that the functional 
integration of the high momentum modes is not disturbed. The expectation value of x in 
the presence of A k S[x] and J reads 

™ - <*•<*» = w (2 ' 4) 

In terms of W k the effective average action is defined via a modified Legendre transform 

= -W k [J] + | d*xJ a (x)<f> a (x) - A k S[<j>]. (2.5) 

The subtraction of the infrared cutoff piece in (2.5) guarantees that the only difference 
between and T is the effective infrared cutoff in the fluctuations. For k > this 
allows the interpretation of T k as a reasonable coarse grained free energy (cf. sections 
2.2 and 4.8). Furthermore, one observes that does not need to be convex whereas 
the standard effective action Y = lim^o T k is a convex functional by its definition as a 
Legendre transform. To establish the property = S one may consider an alternative 
integral equation for T k . Using 

and eqs. (2.1), (2.5) one finds 

exp(-r fc [0]) = J D X exp (-S\x\ ~ ^ k S[ X - 0] + / ^^Ptx" " ^l) ■ ( 2 -7) 

For k > the term exp(— A k S[x — 4>]) constrains the functional integral and behaves 
~ 5[x — 0] for A; A thus leading to the property r fc — > S in this limit [38]. One will 
often not consider the limit A — > oo, but rather keep this scale fixed and associate it with 
a physical ultraviolet cutoff for the validity of the investigated field theory. In particular, 
near second order or weak first order phase transitions the universal properties of the theory 
become independent of the precise form of F\. 

The formulation of the effective average action can be easily generalized to fermionic 
degrees of freedom. In particular, it is possible to incorporate chiral fermions since a 
chirally invariant cutoff R k can be formulated [77]. Possible local gauge symmetries can 
be treated along similar lines [78, 17, 79, 80, 81, 45, 82, 83]. Here A k S may not be gauge 
invariant and the usual Ward identities receive corrections [81] which vanish for k — > 0. 
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2.2 Exact flow equation 

The dependence of the effective average action T k on the coarse graining scale k is described 
by an exact flow equation which can be derived in a straightforward way [38]. From eq. 
(2.5) and (2.4) one finds 

dt % = k-^ r fe | = -dt w k \j - dA k s[<t>] 



l -Tr{w^d t R k }. (2.8) 



The trace in momentum space reads Tr = J2 a I d d q/{2n) d and wjf* denotes the /c-dependent 
connected two-point function 



5J a {-q)5J b {q') 
= (Xa(q) Xb (-q)) ~ (Xa(q))( Xb (-q'))- (2.9) 

With the help of eq. (2.4) differentiated with respect to and (2.6) with respect to J one 
finds that eq. (2.8) can be rewritten as 



^(2)r - - 1-1 



r%>[<l>]+R k d t R k . (2.10) 



d t T k [<p] = ^TT 

The exact flow equation (2.10) describes the scale dependence of T k in terms of the inverse 
average propagator rj^ as given by the second functional derivative of T k with respect to 
the field components 

The trace in eq. (2.10) involves only one momentum integration as well as the summation 
over the internal indices. The additional cutoff function R k with a form like the one given 
in eq. (2.3) renders the momentum integration both infrared and ultraviolet finite. In par- 
ticular, the direct implementation of the additional mass-like term R k ~ Z k k 2 for q 2 <C k 2 
into the inverse average propagator makes the formulation suitable for dealing with theories 
which are plagued by infrared problems in perturbation theory. Up to exponentially small 
corrections the integration of the high momentum modes with q 2 ^> k 2 is not affected by 
the infrared cutoff. We note that the derivation of the exact flow equation does not depend 
on the particular choice of the cutoff function and it may sometimes be technically easier 
to use a different infrared cutoff function, as e.g. R k ~ k 2 . Ultraviolet finiteness, however, 
is related to a fast decay of d t R k for q 2 ^> k 2 . For a rapidly decaying cutoff function 
R k like the one given in eq. (2.3) only fluctuations with momenta in a small "window" 
around q 2 ~ k 2 effectively contribute to the trace on the r.h.s. of eq. (2.10) for given k. 
Accordingly, the large momentum fluctuations with q 2 ^ k 2 , which are not affected by the 
infrared cutoff, have already been integrated out. Since the momentum integrations for 
q 2 ^ k 2 are effectively cut off, T k realizes the concept of a coarse grained free energy in the 
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sense of Langer [71] (cf. section 4.8). Of course, the particular choice for the infrared cutoff 
function should have no effect on the physical results for k — > 0. Different choices of Rk 
correspond to different trajectories in the space of effective actions along which the unique 
infrared limit T is reached. Nevertheless, once approximations are applied not only the 
trajectory but also its end point may depend on the precise definition of the function R k . 
This dependence may be used to study the robustness of the approximation. 

The flow equation (2.10) closely resembles a one-loop equation: Replacing F k by the 
second functional derivative of the classical action, S^ 2 ), one obtains the corresponding 
one-loop result. The "renormalization group improvement" S'*- 2 -' — > rj^ turns the one- 
loop flow equation into an exact non-perturbative flow equation. It also turns the equation 
into a functional differential equation. Possible methods for its solution include standard 
perturbation theory in the case of a small coupling, the 1/A-expansion or the e-expansion. 
Particularly suitable for our purposes is the derivative expansion. We will allow for an arbi- 
trary field dependence of the non-derivative term, i.e. the effective potential. This enables 
us to study the non-analytic behavior of the equation of state near second order phase 
transitions or the simultaneous description of degenerate minima at first order transitions. 
With the help of the derivative expansion the functional differential equation (2.10) can 
be transformed into a coupled set of non-linear partial differential equations. If one does 
not want to resort to further approximations at this level, one needs appropriate tools for 
solving this type of partial differential equations. Analytical solutions can be found only in 
certain limiting cases and for most purposes numerical solutions seem the adequate tool. 
We have developed suitable numerical algorithms which will be addressed to in section 4.5. 

We emphasize that the flow equation (2.10) is closely connected to the Wilsonian (exact) 
renormalization group equation 13 [30, 39, 40, 41, 42, 43]. The latter describes how the 
Wilsonian effective action changes with the ultraviolet cutoff A. Polchinski's continuum 
version of the Wilsonian flow equation [42] can be transformed into eq. (2.10) by means 
of a Legendre transform and a suitable variable redefinition [85, 86]. In contrast to the 
effective average action the Wilsonian effective action does not generate the 1PI Green 
functions [87]. We also note that eq. (2.10) is compatible with all symmetries of the 
model. Extensions of the flow equations to gauge fields [78, 17, 79, 80, 81, 45, 82, 83] and 
fermions [77, 120] (cf. section 5) are available. The exact flow equation can be supplemented 
by an exact formalism for the introduction of composite field variables or, more generally, a 
change of variables [44]. We will exploit this crucial feature for a description of low-energy 
QCD in section 5. 



13 The exact flow equation may also be interpreted as a differential form of the Schwinger-Dyson equations 
[84]. 
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3 Equation of state and second order phase transi- 
tions 

3.1 Introduction 

The prototype for investigations concerning the restoration of a spontaneously broken sym- 
metry at high temperature is the iV-component scalar field theory with O(N) symmetry. 
For N = 4 it describes the scalar sector of the electroweak standard model in the limit of 
vanishing gauge and Yukawa couplings. It is also used as an effective model for the chiral 
phase transition in QCD in the limit of two quark flavors [31, 50, 26, 27]. In condensed 
matter physics N = 3 corresponds to the well known Heisenberg model used to describe 
the ferromagnetic phase transition. There are other applications like the helium super- 
fluid transition (N = 2), liquid- vapor transition (N — 1) or statistical properties of long 
polymer chains (N = 0) [67]. 

The equation of state (EOS) for a magnetic system is specified by the free energy den- 
sity (here denoted by U) as a function of arbitrary magnetization <fi and temperature T. 
All thermodynamic quantities can be derived from the function U((f>,T). For example, 
the response of the system to a homogeneous magnetic field H follows from dU/d<p = H. 
This permits the computation of <fi for arbitrary H and T. There is a close analogy to 
quantum field theory at non-vanishing temperature. Here U corresponds to the temper- 
ature dependent effective potential as a function of a scalar field 0. For instance, in the 
0(4) symmetric model for the chiral phase transition in two flavor QCD the meson field 
(f) has four components. In this picture, the average light quark mass rh is associated with 
the source H ~ rh and one is interested in the behavior during the phase transition (or 
crossover) for H ^ 0. The temperature and source dependent meson masses and zero 
momentum interactions are given by derivatives of U (cf. section 5). 

The applicability of the O(N) symmetric scalar model to a wide class of very different 
physical systems in the vicinity of the critical temperature T c at which the phase transition 
occurs is a manifestation of universality of critical phenomena. There exists a universal 
scaling form of the EOS in the vicinity of the second order phase transition which is not 
accessible to standard perturbation theory. The main reason is the effective three dimen- 
sional behavior of high temperature quantum field theory already discussed in section 1. In 
particular, the critical exponents which describe the singular behavior of various quantities 
near the phase transition are those of the corresponding classical statistical system (cf. 
section 5). 

The quantitative description of the scaling form of the EOS will be the main topic of 
this section [62, 63, 64]. The calculation of the effective potential U(<f),T) in the vicinity 
of the critical temperature of a second order phase transition is an old problem. One can 
prove through a general renormalization group analysis [30] the Widom scaling form [88] 
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of the EOS 14 

H = <f> s f((T-T c )/ ( f> 1 ^). (3.1) 

Only the limiting cases 0-^0 and <fi — > oo are quantitatively well described by critical 
exponents and amplitudes. The critical exponents (3 and S have been computed with high 
accuracy [67, 89, 90] but the scaling function / is more difficult to access. A particular 
difficulty for a perturbative computation in three dimensions arises from the singularities 
induced by massless Goldstone modes in the phase with spontaneous symmetry breaking 
for models with continuous symmetry (N ^ 1). The scaling function / has been computed 
within the e-expansion up to second order in e [60] (third order for N — 1 [67]). Recently, 
the scaling form of the EOS has also been studied for N = 1 and N = 4 using lattice 
Monte-Carlo simulations [66, 65, 61]. We will observe that the lattice results agree quite 
well with those obtained from the effective average action method, while there appears a 
considerable discrepancy to the e-expansion results. 

In the following we compute the scaling EOS for O(N) symmetric scalar TV-component 
models directly in three dimensions [62, 63, 64]. Section 3.2 presents a computation of 
the effective potential U = lim fc ^ Uk from a derivative expansion of the effective average 
action. The approximation takes into account the most general field dependence of the 
potential term. This will allow us to compute the non-analytic behavior of U in the vicinity 
of the second order phase transition. From U the universal scaling form of the equation of 
state is extracted in section 3.3. We conclude in section 3.4. 

3.2 Scale dependence of the effective average potential 

In this section we compute the effective average potential Uk(p) for an O(N) symmetric 
scalar field theory directly in three dimensions [62, 63, 64]. Here p = |0 a a and <p a 
denotes the iV-component real scalar field. For k — > one obtains the effective potential 
(Helmholtz free energy) U (p) = U(p). In the phase with spontaneous symmetry breaking 
the minimum of the potential occurs for k = at po 7^ 0. In the symmetric phase the 
minimum of Uk(p) ends at po — for k — 0. The two phases are separated by a scaling 
solution for which Uk/k 3 becomes independent of k once expressed in terms of a suitably 
rescaled field variable and the corresponding phase transition is of second order. 

Though the evolution equation (2.10) for the effective average action is exact, it remains 
a complicated functional differential equation. In practice, one has to find a truncation 
for Tk in order to obtain approximate solutions. An important feature of the exact flow 
equation is therefore its simple and intuitive form which helps to find a non-perturbative 
approximation scheme. The r.h.s. of eq. (2.10) expresses the scale dependence of T k in 
terms of the exact propagator. Known properties of the propagator can be used as a guide 
to find an appropriate truncation for the effective average action. For a scalar theory the 
propagator is a matrix characterized by mass terms and kinetic terms ~ Zq 2 . The mass 
matrix is given by the second derivative of the potential Uk with respect to the fields. In 

14 We frequently suppress in our notation an appropriate power of a suitable "microscopic" length scale 
A^ 1 which is used to render quantities dimensionless. 
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general Z can be a complicated function of the fields and momenta. We may exploit the 
fact that the function Z plays the role of a field and momentum dependent wave function 
renormalization factor. For second order phase transitions and approximately for weak 
first order phase transitions the behavior of Z is governed by the anomalous dimension 77. 
Typically for three dimensional scalar theories r\ is small. One therefore expects a weak 
dependence of Z on the fields and momenta (cf. also the discussion in section 6). 

Our truncation is the lowest order in a systematic derivative expansion of T k [37, 55, 
92, 93, 94, 62, 63, 95] 

T k = J d d x{u k {p) + X -Z k d^ a d^ a ). (3.2) 

We keep for the potential term the most general O(N) symmetric form U k (p), whereas 
the wave function renormalization is approximated by one /c-dependent parameter. Next 
order in the derivative expansion would be the generalization to a p-dependent wavefunc- 
tion renormalization Z k (p) plus a function Y k (p) accounting for a possible different index 
structure of the kinetic term for iV > 2 [37, 92, 95] (cf. section 6). Going further would 
require the consideration of terms with four derivatives and so on. The weak p-dependence 
of the wave function renormalization factor employed in the ansatz (3.2) has been estab- 
lished explicitly for the scaling solution [93, 95] and also recently for the equation of state 
(N = 1) [96]. ' 

If the ansatz (4.2) is inserted into the evolution equation for the effective average action 
(2.10) one can extract flow equations for the effective average potential U k (p) and for the 
wave function renormalization constant Z k (or equivalently the anomalous dimension rj) 15 . 
For a study of the behavior in the vicinity of the phase transition it is convenient to work 
with dimensionless renormalized fields 16 

p = Z k k 2 ~ d p , 
u k (p) = k- d U k (p). (3.3) 

With the truncation of eq. (3.2) the exact evolution equation for u' k = du k /dp [55, 92, 63] 
reduces then to the partial differential equation 

°-^= (-2 + r))u' k + (d-2 + r))pu'l 

-2v d (N - 1K/?K; 77) - 2v d (3u't + 2~pu' k ")l d {u' k + 2p<; v), (3.4) 

where t = ln(/c/A), with A the ultraviolet cutoff of the theory. The anomalous dimension 
i] is defined by 

d 

r) = -—\nZ k (3.5) 

and 

vt = ($) , (3.6) 



15 Details of the calculation may be found in section 4 for the discussion of scalar matrix models which 
covers the considerations involved here. 

16 We keep the number of dimensions d arbitrary and specialize only later to d = 3. 
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with v 3 = 1/87T 2 . The "threshold" functions l d (w; 77) result from the momentum integration 
on the r.h.s. of eq. (2.10), and read for n > 1, with y = q 2 /k 2 

l d n (w; rj) = -n f° dyy^ 1 ^ [y(l + r{y)) + w]~ {n+1) 
Jo ay 



\ I dyyir(y) [y(l + r(y)) + w]^ n+1) . (3.7) 
2 Jo 

Here r(y) depends on the choice of the momentum dependence of the infrared cutoff and 
we employ according to eq. (2.3) 



roo 



r (y) = Y~~^-y' (3 ' 8) 

This choice has the property lim 9 2^ R k /Z k k 2 = \im y ^ yr(y) = 1, whereas for q 2 3> k 2 
the effect of the infrared cutoff is exponentially suppressed. The "threshold" functions l d 
decrease rapidly for increased w and account for the decoupling of modes with mass term 
w ^> 1. We evaluate these functions numerically. Finally, the anomalous dimension is 
given in our truncation by 17 [55, 92, 63] 

r)(k) = ^j±K\ 2 m d 2)2 (2\K), (3.9) 
with k the location of the minimum of the potential and A the quartic coupling 



<(«) = A . (3.10) 



The function toJj 2 is given by 



m 



poo . 

22 {w) = / dyy~ 2 - 
Jo (1 



(1 + r) 2 {(1 + r)y + w] 

dr ( d \ 2 o / dr\ dr 

2y— + 2 y— r - 2y 2 [l + r + y—) — 



+ 



(1 + r)y (1 + r)y + w 



(3.11) 



We point out that the argument 2\k turns out generically to be of order one for the 
scaling solution. Therefore, k ~ A -1 and the mass effects are important, in contrast to 
perturbation theory where they are treated as small quantities ~ A. 

At a second order phase transition there is no mass scale present in the theory. In 
particular, one expects a scaling behavior of the rescaled effective average potential Uk(p)- 
This can be studied by following the trajectory describing the scale dependence of Uk(p) 

17 We neglect here for simplicity the implicit, linear 77-dependence of the function 2 - We have numer- 
ically verified this approximation to have only a minor effect on the value of rj. 
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Figure 1: The evolution of u' k (p) as k is lowered from A to zero for N = 1. The initial conditions (bare 
couplings) have been chosen such that the scaling solution is approached before the system evolves towards 
the symmetric phase with u' k (0) > 0. 

as k is lowered from A to zero. Near the phase transition the trajectory spends most of 
the "time" t in the vicinity of the /c-independent scaling solution of eq. (3.4) given by 
d t u'^(p) = 0. 18 Only at the end of the running the "near-critical" trajectories deviate 
from the scaling solution. For k — > they either end up in the symmetric phase with 
k = and positive constant mass term m 2 such that u' k (0) ~ m 2 /k 2 ; or they lead to a 
non-vanishing constant po indicating spontaneous symmetry breaking with k — > Z k 2 ~ d p . 
The equation of state involves the potential Uo(p) for temperatures away from the critical 
temperature. Its computation requires the solution for the running away from the critical 
trajectory which involves the full partial differential equation (3.4). We have developed 
two alternative algorithms for its numerical solution which will be discussed in section 4.5. 

In fig. 1 we present the results of the numerical integration of eq. (3.4) for d = 3 and 
N — 1. The function u' k (p) is plotted for various values of t — \n(k/A). The evolution 
starts at k — A (t — 0) where the average potential is equal to the classical potential 
(no effective integration of modes has been performed). We start with a quartic classical 
potential parameterized as 

u' A (p) = X A (p-K A ). (3.12) 

We arbitrarily choose Aa = 0.1 and fine tune k a so that a scaling solution is approached 
at later stages of the evolution. There is a critical value k ct ~ 6.396 x 10~ 2 for which 
the evolution leads to the scaling solution. For the results in fig. 1 a value k a slightly 
smaller than K cr is used. As k is lowered (and t turns negative), u' k (p) deviates from its 

18 Thc resulting ordinary differential equation for u*(p) has already been solved numerically for a some- 
what different choice of the infrared cutoff [93, 94, 95]. 
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initial linear shape. Subsequently it evolves towards a form which is independent of k and 
corresponds to the scaling solution dtu'^(p) = 0. It spends a long "time" t - which can 
be rendered arbitrarily long through appropriate fine tuning of k\ - in the vicinity of the 
scaling solution. During this "time", the minimum of the potential u' k {p) takes a fixed 
value re*, while the minimum of U k (p) evolves towards zero according to 

p (k) = kK*/Z k . (3.13) 

The longer u k (p) stays near the scaling solution, the smaller the resulting value of po(k) 
when the system deviates from it. As this value determines the mass scale for the renor- 
malized theory at k — 0, the scaling solution governs the behavior of the system very 
close to the phase transition, where the characteristic mass scale goes to zero. Another 
important property of the "near-critical" trajectories, which spend a long "time" t near 
the scaling solution, is that they become insensitive to the details of the classical theory 
which determine the initial conditions for the evolution. After u' k (p) has evolved away from 
its scaling form u'^(p), its shape is independent of the choice of Aa for the classical theory. 
This property gives rise to the universal behavior near second order phase transitions. For 
the solution depicted in fig. 1, Uk(p) evolves in such a way that its minimum runs to zero 
with u' k (0) subsequently increasing. Eventually the theory settles down in the symmetric 
phase with a positive constant renormalized mass term m 2 = k 2 u' k (0) as k — > 0. Another 
possibility is that the system ends up in the phase with spontaneous symmetry breaking. 
In this case k grows in such a way that po(k) approaches a constant value for k — > 0. 

The approach to the scaling solution and the deviation from it can also be seen in fig. 
2. The evolution of the running parameters n(t), X(t) starts with their initial classical 
values, leads to fixed point values «*, A* near the scaling solution, and finally ends up in 
the symmetric phase (k runs to zero). Similarly the anomalous dimension f](k), which is 
given by eq. (3.9), takes a fixed point value 77* when the scaling solution is approached. 
During this part of the evolution the wave function renormalization is given by 

Z k ~ AT* (3.14) 

according to eq. (3.5). When the parts of the evolution towards and away from the fixed 
point become negligible compared to the evolution near the fixed point - that is, very close 
to the phase transition - eq. (3.14) becomes a very good approximation for sufficiently low 
k. This indicates that 77* can be identified with the critical exponent 77. For the solution 
of fig. 2 (N = 1) we find = 4.07 x 1(T 2 , A* = 9.04 and 77* = 4.4 x 10" 2 . 

As we have already mentioned the details of the renormalized theory in the vicinity of 
the phase transition are independent of the classical coupling Aa- Also the initial form of the 
potential does not have to be of the quartic form of eq. (3.12) as long as the symmetries 
are respected. Moreover, the critical theory can be parameterized in terms of critical 
exponents [56], an example of which is the anomalous dimension 77. These exponents are 
universal quantities which depend only on the dimensionality of the system and its internal 
symmetries. For our three-dimensional theory they depend only on the value of N and 
can be easily extracted from our results. We concentrate here on the exponent u, which 
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parameterizes the behavior of the renormalized mass in the critical region. Other exponents 
are computed in sections 3.3 and 5.7 along with the critical equation of state. The other 
exponents are not independent quantities, but can be determined from r\ and v through 
universal scaling laws [56]. We define the exponent v through the renormalized mass term 
in the symmetric phase 

m 2 = 1-<£M. = k 2 u'(0) for k - 0. (3.15) 
Z k dp 

The behavior of m 2 in the critical region depends only on the distance from the phase 
transition, which can be expressed in terms of the difference of ka from the critical value 
k ct for which the renormalized theory has exactly m 2 = 0. The exponent v is determined 
from the relation 

m 2 ~ |5/tA| 2iy = \k\ — K cr \ 2u . (3.16) 

For a determination of v from our results we calculate m 2 for various values of near 
k ct . We subsequently plot ln(m 2 ) as a function of In \8ka\- This curve becomes linear for 
Sk\ — > and we obtain v from the constant slope. In the past the critical exponents of the 
O(N) symmetric theory were calculated from truncated versions of the partial differential 
equation (3.4) [55, 92]. The strategy was to turn eq. (3.4) into an infinite system of ordi- 
nary differential equations for the coefficients of a Taylor expansion around the "running" 
minimum of the potential. This infinite system was approximately solved by neglecting 
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p-derivatives of u k (p) higher than a given order. The apparent convergence of the proce- 
dure was checked by enlarging the level of truncation. We now have an alternative way of 
estimating the accuracy of this method. Our numerical solution of the partial differential 
equation (3.4) corresponds to an infinite level of truncation where all the higher derivatives 
are taken into account. In table 1 we present results obtained through the procedure of 
successive truncations and through our numerical solution for N = 3. We give the values 





ft* 


A* 


(3) 
Ui 


V 


V 


a 


6.57 x KT 2 


11.5 






0.745 


b 


8.01 x KT 2 


7.27 


52.8 




0.794 


c 


7.86 x 10~ 2 


6.64 


42.0 


3.6 x 10~ 2 


0.760 


d 


7.75 x 1(T 2 


6.94 


43.5 


3.8 x 10~ 2 


0.753 


e 


7.71 x 1(T 2 


7.03 


43.4 


3.8 x 10~ 2 


0.752 


f 


7.64 x 1(T 2 


7.07 


44.2 


3.8 x 10~ 2 


0.747 



Table 1: The minimum k of the potential u k (p), the derivatives A = u"(k), u k (k) for the 
scaling solution, and the critical exponents rj and v, in various approximations: (a)-(e) 
from refs. [55, 92] and (f) from the present section [63, 62]. N — 3. 

a) Truncation where only the evolution of k and A is considered and higher derivatives of 
the potential and the anomalous dimension are neglected. 

b) k, A, u^\k) are included. 

c) k, A, u k (k) are included and rj is approximated by eq. (3.9). 

d) with five parameters: k, A, u k 3 \n), u k 4 \n) and rj. 

e) as in d) and in addition u k 5 \n), u k % \n) are estimated. 

f) The partial differential equation (3.4) for u' k {p) is solved numerically and rj is approxi- 
mated by eq. (3.9). 

of k, A, u k 3 \n) for the scaling solution and the critical exponents rj, v. We observe how the 
results stabilize as more p-derivatives of u k (p) at p — k and the anomalous dimension are 
taken into account. The last line gives the results of our numerical solution of eq. (3.4). 
By comparing with the previous line we conclude that the inclusion of all the p-derivatives 
higher than u k 6 \n) and the term ~ r\ in the "threshold" function of eq. (3.7) generates an 
improvement of less than 1 % for the results. This is smaller than the error induced by 
the omission of the higher derivative terms in the average action, which typically gener- 
ates an uncertainty of the order of the anomalous dimension. In table 2 we compare our 
values for the critical exponents with results obtained from other methods 19 (such as the 

19 See also ref. [95] and references therein for a calculation of critical exponents using similar methods as 
in this work. 
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e-expansion, perturbation series at fixed dimension, lattice high temperature expansions 
and the 1/iV-expansion) . As expected r\ is rather poorly determined since it is the quan- 



N 


V 


V 





0.589 


0.5880(15) a 
0.5880(15) 6 
0.5878(6) c 


0.040 


0.027(4) a 
0.0320(25) 6 


1 


0.643 


0.6300(15) a 
0.6310(15) b 
0.6315(8) c 


0.044 


0.032(3) a 
0.0375(25) b 


2 


0.697 


0.6695(20) a 
0.671(5) b 
0.675(2) c 


0.042 


0.033(4) a 
0.040(3) 6 


3 


0.747 


0.705(3) a 
0.710(7) b 
0.716(2) c 


0.038 


0.033(4) a 
0.040(3) 6 


4 


0.787 


0.73(2) a * 
0.759(3) c 


0.034 


0.03(l) a * 


10 


0.904 


0.894(4) c 
0.877 d 


0.019 


0.025 d 


100 


0.990 


0.989^ 


0.002 


0.003 d 



Table 2: Critical exponents v and rj for various values of N. For comparison we list results 
obtained with other methods as summarized in [67], [97] and [89]: 

a) From perturbation series at fixed dimension to six-loop ( * seven-loop [98]) order. 

b) From the e-expansion at order e 5 . 

c) From lattice high temperature expansions. 

d) From the 1 /iV-expansion at order 1/N 2 . 

tity most seriously affected by the omission of the higher derivative terms in the average 
action. The exponent v is in agreement with the known results at the 1-5 % level, with a 
discrepancy roughly equal to the value of r\ for various N. 

In conclusion, the shape of the average potential is under good quantitative control for 
every scale k. This permits a quantitative understanding of the most important properties 
of the system at every length scale. We will exploit this in the following to extract the 
scaling form of the equation of state which is the main result here. 

3.3 Universal critical equation of state 

In this section we extract the Widom scaling form of the equation of state (EOS) from a 
solution [62, 64] of eqs. (3.4), (3.9) for the three dimensional O(N) model. Its asymptotic 
behavior yields the universal critical exponents and amplitude ratios. We also present fits 
for the scaling function for N — 1 and N = 3. (See in addition section 5.7 for iV = 4.) An 
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Figure 3: Logarithmic plot of / and df jdx for x > 0. 



alternative parameterization of the equation of state in terms of renormalized quantities is 
used in order to compute universal couplings. 

Eq. (3.1) establishes the scaling properties of the EOS. The external field H is related 
to the derivative of the effective potential U' = dU /dp by H a = U'ifia- The critical equation 
of state, relating the temperature, the external field and the order parameter, can then be 
written in the scaling form (0 = y/2p) 



TXT = /( x )> 



x = 



-5k a 



(3.17) 



with critical exponents 5 and (3. A measure of the distance from the phase transition is the 
difference 5k, a — k,\ — n cr . If is interpreted as a function of temperature, the deviation 
Ska is proportional to the deviation from the critical temperature, i.e. 5k\ = A(T)(T cr — T) 
with A(T cr ) > 0. For <fi — > oo our numerical solution for U' obeys U' ~ with high 
accuracy. The inferred value of 5 is displayed in table 3, and we have checked the scaling 
relation 5 — (5 — 77) / (1 + 77). The value of the critical exponent 77 is obtained from eq. (3.4) 
for the scaling solution. We have also verified explicitly that / depends only on the scaling 
variable x for the value of (3 given in table 3. In figs. 1 and 2 we plot log(/) and log(df/dx) 
as a function of log|a;| for iV = 1 and N = 3. Fig. 1 corresponds to the symmetric phase 
(x > 0), and fig. 2 to the phase with spontaneous symmetry breaking (x < 0). This is our 
main result of this section. 

One can easily extract the asymptotic behavior from the logarithmic plots and compare 
with known values of critical exponents and amplitudes. The curves become constant, both 
for x — > + and x — > CT with the same value, consistently with the regularity of f(x) at 
x — 0. For the universal function one obtains 



lim f(x) 

x->0 J V ' 



D, 



(3.18) 
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Figure 4: Logarithmic plot of / and df/dx for x < 0. 



and H = Dtp 6 on the critical isotherm. For x — > oo one observes that log(/) becomes a 
linear function of log(x) with constant slope 7. In this limit the universal function takes 
the form 

Jm/(s) = (C+)-V. (3.19) 

The amplitude C + and the critical exponent 7 characterize the behavior of the 'unrenor- 
malized' squared mass or inverse susceptibility 

m 2 = X- 1 = Jim (0 = (C+r^lV- 1 -^) • (3.20) 

We have verified the scaling relation 7//? = 5 — 1 that connects 7 with the exponents /3 and 
S appearing in the Widom scaling form (3.1). One observes that the zero- field magnetic 
susceptibility, or equivalently the inverse unrenormalized squared mass fnT 2 = x, is non- 
analytic for Ska — > in the symmetric phase: x — C + |^ k a| 7 - In this phase we find that 
the correlation length £ = (Zqx) 1 ^ 2 , which is equal to the inverse of the renormalized mass 
m R , behaves as £ = £ + ]5ka| -1/ with v — 7/(2 — r\). 

The spontaneously broken phase is characterized by a nonzero value 0o of the minimum 
of the effective potential U with H = (dU/d<f>)(<f>o) = 0. The appearance of spontaneous 
symmetry breaking below T c implies that f(x) has a zero x = —B~ l IP and one observes 
a singularity of the logarithmic plot in fig. 2. In particular, according to eq. (3.1) the 
minimum behaves as 0o — B{5kjC) 13 . Below the critical temperature, the longitudinal and 
transversal susceptibilities xl and xt are different for iV > 1 

^=w = * s '^ m - x p fix ^ x - = \% = ^ 1}(x) (3 - 21) 

(with /' = df/dx). This is related to the existence of massless Goldstone modes in the 
(N — 1) transverse directions, which causes the transversal susceptibility to diverge for 
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vanishing external field. Fluctuations of these modes induce the divergence of the zero- 
field longitudinal susceptibility. This can be concluded from the singularity of log(/') for 
N = 3 in fig. 2. The first x-derivative of the universal function vanishes as H — > 0, 
i.e. f{x = —B~ l IP) = for N > 1. For N — 1 there is a non-vanishing constant 
value for f(x = —B^ 1 ^) with a finite zero-field susceptibility x — C~(^ k a)~ 7 > where 
(C - ) -1 = B s ~ 1-1 ^ f\— B~ l IP)/(3. For a non-vanishing physical infrared cutoff k, the 
longitudinal susceptibility remains finite also for N > 1: Xl ~ (^Po) F° r N — 1 in 
the ordered phase, the correlation length behaves as £ = £~(5ka)~ u , and the renormalized 
minimum p 0R = Z p of the potential U scales as p i? — E{5k^) u . 

The amplitudes of singularities near the phase transition D, C ± , £ ± , B and E are given 
in table 3. They are not universal. There are two independent scales in the vicinity of 
the transition point which can be related to the deviation from the critical temperature 
and to the external source. All models in the same universality class can be related by a 
multiplicative rescaling of and 5ka or (T c — T). Accordingly there are only two inde- 
pendent amplitudes and exponents respectively. Ratios of amplitudes which are invariant 
under this rescaling are universal. We display the universal combinations C + /C~, 
R x = C + DB 5 -\ = (£+)/3/^ 1 /(<5+i) j b anc i in ta bl e 3. 
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0.747 


0.038 


7.45 


57.4 





~ 250 






C+ 


D 


B 


e 


E 


c+/c- 


r/r 


R x 


h 


£+E 


1 


0.0742 


15.88 


1.087 


0.257 


0.652 


4.29 


1.86 


1.61 


0.865 


0.168 
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0.0743 


8.02 


1.180 


0.263 


0.746 






1.11 


0.845 


0.196 



Table 3: Parameters for the equation of state. 



The asymptotic behavior observed for the universal function can be used in order to 
obtain a semi-analytical expression for fix). We find that the following fit reproduces the 
numerical values for both / and df/dx with a 1% accuracy (apart from the immediate 
vicinity of the zero of / for iV = 3): 



f fit (x) = D(l + B 1//3 x y(l + 0a;) A (l + cx) 



7— a— A 



(3.22) 



with c = (C+ J D J B«//?e A )- 1 /(7-a-A) an d a = 1 ( a = 2) for N = 1 (N > 1). The fitting 
parameters are chosen as = 0.569 (1.312) and A = 0.180 (-0.595) for TV = 1 (3). 

There is an alternative parameterization of the equation of state in terms of renormal- 
ized quantities. In the symmetric phase (<5ka < 0) we consider the dimensionless quantity 



F(s) 



m 



C + x-ifix), 



R 



mji z 



(3.23) 
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with p R = Z p and U R n ^ = Z n U < " n \ The derivatives of F at s = yield the universal 
couplings 

^ EM^^ ^(0) = ^(0) S ^, (3.24) 
as m,R rriR ds z 

and similarly for higher derivatives. We confirm that the potential U is well approximated 
by a 6 (or p 3 ) polynomial [55, 92, 66]. For N = 1, our value for vr in the symmetric 
phase agrees with the results of a numerical simulation [66] within the expected accuracy 
(related to the size of the anomalous dimension 77). However, Xr/uir seems to be too small. 
This may be connected to the more crude approximation for the anomalous dimension in 
the symmetric phase, as compared to the ordered one. For a recent detailed comparison 
of the dimensionless renormalized couplings in the symmetric phase with results obtained 
from other methods we refer the reader to ref. [99]. Small field expansions of the effec- 
tive potential have also been used to reconstruct the equation of state (N — 1) within 
summed perturbation series at fixed dimension or from high temperature expansions [100]. 
A comparison including the results obtained here can be found in ref. [100]. 
In the ordered phase (5k a > 0) we consider the ratio 

G(s) = ^= l -B->E-*{-x)-y{x), s=^ = B-\-x)^. (3.25) 

PoR 1 POR 

The values for the universal couplings 

dG = Uj^A s , = K(pm) = Or, (3.26) 

as p R p R ds z 

as well as Xr/uir and u R are given in table 3. One observes that for iV > 1 the renormalized 
quartic coupling Xr vanishes in the ordered phase. The coupling u R takes a finite value and 
higher derivatives d n G/ds n (l) (for n < 3) diverge. This is due to the presence of massless 
fluctuations. 

For comparison, figure 5 is taken from ref. [65] which presents results from a lattice 
Monte Carlo investigation of the equation of state for iV = 1 (curve V), along with our 
results (curve W) and those from the e-expansion to order e°, . . . , e 3 (curves — 3). The 
plot shows very good agreement of the lattice study with our results and the fit of eq. (3.22) 
in the phase with spontaneous symmetry breaking, whereas it deviates substantially from 
the results of the e-expansion. We also have computed the equation of state for iV = 4. We 
present these results in section 5.7 where we obtain the universal critical equation of state 
within a study of an effective quark meson model at non-zero temperature. Comparison 
of our iV = 4 results with a recent lattice Monte-Carlo simulation [61] shows agreement for 
the scaling function within a few percent in the symmetric and the spontaneously broken 
phase (cf. section 5.7). 
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In summary, our numerical solution of eq. (3.3) gives a very detailed picture of the critical 
equation of state of the three dimensional O(N) model. The numerical uncertainties are es- 
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Figure 5: The plot shows the "Ising equation of state" (N = 1) computed by different methods. The 
external source j = dU /d<j) is plotted as a function of the magnetization <f>. In contrast to the parameters 
used in the text here all curves are normalized by the conditions dU /d<j)(l) — and d 2 U / d((f) 2 (1) = 1. 
The curves — 3 are from the e-expansion in the parametric representation of the equation of state [67] , 
in orders e°, . . . , e 3 , respectively. The curve W represents our result and the curve V denotes the lattice 
result of ref. [65] where this plot is taken from. 



timated by comparison of results obtained through two independent integration algorithms 
[63, 68] (cf. section 4.5). They are small, typically less than 0.3% for critical exponents and 
(1 - 3)% for amplitudes. The scaling relations between the critical exponents are fulfilled 
within a deviation of 2 x 10~ 4 . The dominant quantitative error stems from the truncation 
of the exact flow equation and is related to the size of the anomalous dimension rj ~ 4% 
(cf. also the discussion in section 6). 

We have obtained a precise description of the universal equation of state in the vicinity 
of the critical temperature of a second order phase transition. We will see in the next 
section that universal properties can also be observed for weak first order phase transitions. 
In particular, we will compute the universal equation of state for weak first order phase 
transitions in scalar matrix models in section 4.7. 
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4 Equation of state and first order phase transitions 
4.1 Introduction 

Matrix models are extensively discussed in statistical physics. Beyond the O(N) symmetric 
Heisenberg models ("vector models"), which we have discussed in the previous section, they 
correspond to the simplest scalar field theories. There is a wide set of different applications 
as the metal insulator transition [101] or liquid crystals [102] or strings and random surfaces 
[103]. The universal behavior of these models in the vicinity of a second order or weak 
first order phase transition is determined by the symmetries and the field content of the 
corresponding field theories. We will consider here [68] models with U(N) x U (N) symmetry 
with a scalar field in the (N, N) representation, described by an arbitrary complex N x N 
matrix if. We do not impose nonlinear constraints for ip a priori but rather use a "classical" 
potential. This enforces nonlinear constraints in certain limiting cases. Among those, our 
model describes a nonlinear matrix model for unitary matrices or one for singular 2x2 
matrices. The universal critical behavior does not depend on the details of the classical 
potential and there is no difference between the linear and nonlinear models in the vicinity 
of the limiting cases. We concentrate here on three dimensions, relevant for statistical 
physics and critical phenomena in high temperature field theory. 

The cases N = 2, 3 have a relation to high temperature strong interaction physics. 
At vanishing temperature the four dimensional models can be used for a description of 
the pseudoscalar and scalar mesons for iV quark flavors (cf. section 5). For N = 3 the 
antihermitean part of ip describes here the pions, kaons, i] and rj' whereas the hermitean 
part accounts for the nonet of scalar ++ mesons. 20 For nonzero temperature T the effects of 
fluctuations with momenta p 2 ^ (ttT) 2 are described by the corresponding three dimensional 
models. These models account for the long distance physics and obtain after integrating 
out the short distance fluctuations by virtue of dimensional reduction which will be treated 
in detail in section 5.5. In particular, the three dimensional models embody the essential 
dynamics in the immediate vicinity of a second order or weak first order chiral phase 
transition [31, 50, 26, 27]. The four dimensional models at nonvanishing temperature have 
also been used for investigations of the temperature dependence of meson masses [105, 52]. 
The simple model investigated in this section is not yet realistic - it neglects the effect 
of the axial anomaly which reduces the chiral flavor symmetry to SU(N) x SU(N). For 
simplicity we will concentrate here on N = 2, but our methods can be generalized to N = 3 
and the inclusion of the axial anomaly. 

The case N = 2 also has a relation to the electroweak phase transition in models with 
two Higgs doublets. Our model corresponds here to the critical behavior in a special class 
of left-right symmetric theories in the limit where the gauge couplings are neglected. Even 
though vanishing gauge couplings are not a good approximation for typical realistic models 
one would like to understand this limiting case reliably. 

For the present matrix models one wants to know if the phase transition becomes 

20 See ref. [104] for a recent phenomenological analysis. 
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second order in certain regions of parameter space. In the context of flow equations this is 
equivalent to the question if the system of running couplings admits a fixed point which is 
infrared stable (except one relevant direction corresponding to T— T c ). We will employ here 
the average action method introduced in section 2. We find that the phase transition for 
the investigated matrix models with N = 2 and symmetry breaking pattern U(2)xU (2) — > 
U(2) is always (fluctuation induced) first order, except for a boundary case with enhanced 
0(8) symmetry. For a large part of parameter space the transition is weak and one finds 
large renormalized dimensionless couplings near the critical temperature. If the running 
of the couplings towards approximate fixed points (there are no exact fixed points) is 
sufficiently fast the large distance physics looses memory of the details of the short distance 
or classical action. In this case the physics near the phase transition is described by an 
universal equation of state. 

Besides the possible practical applications it is a theoretical challenge to find the uni- 
versal equation of state for weak first order phase transitions. Whereas for second order 
transitions the universal equation of state can be expressed as a function of only one scal- 
ing variable (cf. section 3.3), the additional mass scale at a first order transition induces a 
second dimensionless ratio. The equation of state or, equivalently, the effective potential 
or free energy therefore depends on two different scaling variables in the universal region. 
We have succeeded to describe this situation for the present matrix models and give the 
universal equation of state in section 4.7. 

An important aspect in the study of first order phase transitions is the possible cre- 
ation of an out of equilibrium situation, as has been discussed in section 1. An old and 
challenging problem in this context concerns the question of the validity of Langer's the- 
ory [71] of spontaneous bubble nucleation. In Langer's approach the decay of metastable 
states is described in terms of a steady state solution of the associated Fokker-Planck 
equation. The nucleation rate is obtained from the integrated probability current flowing 
across the saddle point ( "critical bubble" ) of a suitably defined coarse grained free energy. 
The treatment requires first a meaningful definition of a coarse grained free energy with 
a coarse graining scale k. Second, the validity of a saddle point approximation for the 
treatment of fluctuations around the critical bubble has to be ensured. Here only fluctu- 
ations with momenta smaller than k must be included. We will see that the two issues 
are closely related. The validity of the saddle point approximation typically requires small 
dimensionless couplings. On the other hand we observe for large couplings that the form 
of the coarse grained effective potential £4 depends strongly on k even for scales where 
the location of the minima of Uk is essentially independent of k. This means that the 
lowest order in the saddle point approximation (classical contribution) depends strongly 
on the details of the coarse graining procedure. Since the final results as nucleation rates 
etc. must be independent of the coarse graining prescription this is only compatible with a 
large contribution from the higher orders of the saddle point expansion. Section 4.8 deals 
with this issue in a quantitative way. 

The text is organized as follows. In section 4.2 we define the U{2) x U(2) symmetric 
matrix model and we establish the connection to a matrix model for unitary matrices and 
to one for singular complex 2x2 matrices. There we also give an interpretation of the 
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model as the coupled system of two 5?7(2)-doublets for the weak interaction Higgs sector. 
Section 4.3 is devoted to an overview over the phase structure and the coarse grained 
effective potential Uk for the three dimensional theory. The evolution equation for Uk and 
its scaling form is computed in section 4.4. A method for its numerical solution is discussed 
in section 4.5 which also contains the flow equation for the wave function renormalization 
constant or the anomalous dimension rj. A detailed account on the renormalization group 
flow is presented in section 4.6. We compute the universal form of the equation of state for 
weak first order phase transitions in section 4.7 and we extract critical exponents and the 
corresponding index relations. The dependence of the coarse grained effective potential on 
the coarse graining scale is studied in detail in section 4.8. The quantitative analysis is 
applied to Langer's theory of bubble formation. Section 4.9 contains the conclusions. 



4.2 Scalar matrix model 

We consider a,U(2)xU (2) symmetric effective action for a scalar field tp which transforms in 
the (2, 2) representation with respect to the subgroup 577(2) x 577(2). Here tp is represented 
by a complex 2x2 matrix and the transformations are 

tp -> UtpV^ , 

tp\ vtp-^W 

where U and V are unitary 2x2 matrices corresponding to the two distinct U (2) factors. 

We classify the invariants for the construction of the effective average action by the 
number of derivatives. The lowest order in a systematic derivative expansion [92, 93] of T k 
is given by 

Y k = J d d x {U k (tp, tp^) + Zkd,tp* ab d^tp ab } (a, b=l,2). (4.2) 

The term with no derivatives defines the scalar potential Uk which is an arbitrary function 
of traces of powers of tp^tp. The most general U(2) x U (2) symmetric scalar potential can 
be expressed as a function of only two independent invariants, 

p = tr(y?V) 

r = 2tr(^V--p) =2tr(y?V) - P 2 ■ 

Here we have used for later convenience the traceless matrix tp^tp — |p to construct the 
second invariant. Higher invariants, tr (tp^tp — |p) for n > 2, can be expressed as functions 
of p and t [46]. 

For the derivative part we consider a standard kinetic term with a scale dependent 
wave function renormalization constant Z k - The first correction to the kinetic term would 
include field dependent wave function renormalizations Z k (p, r) plus functions not specified 
in eq. (4.2) which account for a different index structure of invariants with two derivatives. 
These wave function renormalizations may be defined at zero momentum. The next level 
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involves invariants with four derivatives and so on. We define Z k at the minimum p , tq of 
U k and at vanishing momenta q 2 , 

Z k = Z k (p = p ,r = T -q 2 = 0). (4.4) 

The factor Z k appearing in the definition of the infrared cutoff R k in eq. (2.3) is identified 
with (4.4). The /^-dependence of this function is given by the anomalous dimension 

V (k) = ~j^Z k . (4.5) 

If the ansatz (4.2) is inserted into the flow equation for the effective average action 
(2.10) one obtains flow equations for the effective average potential U k (p,r) and for the 
wave function renormalization constant Z k (or equivalently the anomalous dimension rj). 
This is done in sections 4.4 and 4.5. These flow equations have to be integrated starting 
from some short distance scale A and one has to specify Ua and Z A as initial conditions. At 
the scale A, where I\ can be taken as the classical or short distance action, no integration 
of fluctuations has been performed. The short distance potential is taken to be a quartic 
potential which is parametrized by two quartic couplings Aia, A 2 a and a mass term. We 
start in the spontaneously broken regime where the minimum of the potential occurs at a 
nonvanishing field value and there is a negative mass term at the origin of the potential 

04 >o), 

U A (p, t) = -p\p + -A 1A p 2 + -A 2A r (4.6) 

and Z\ — 1. The potential is bounded from below provided Aia > and A 2 a > _ 2Aia- For 
A 2 a > one observes the potential minimum for the configuration ip a b = ip5 a b corresponding 
to the spontaneous symmetry breaking down to the diagonal U (2) subgroup of U(2)xU (2). 
For negative A 2 a the potential is minimized by the configuration = (p5 a \5 a b which 
corresponds to the symmetry breaking pattern U(2) x U(2) — > U(l) x (/(l) x U(l). In 
the special case A 2 a = the theory exhibits an enhanced 0(8) symmetry. This constitutes 
the boundary between two phases with different symmetry breaking patterns. 

The limits of infinite couplings correspond to nonlinear constraints in the matrix model. 
For Aia — > oo with fixed ratio p\/Xi\ one finds the constraint tr = fi\/Xi\- By a 

convenient choice of Z\ (rescaling of tp) this can be brought to the form tr(</?'</?) = 2. On 
the other hand, the limit A 2 a — > +oo enforces the constraint ip^<p = | tr (ip^ if). Combining 
the limits Aia - ► oo, A 2 a — > oo the constraint reads (p^ip = 1 and we deal with a matrix 
model for unitary matrices. (These considerations generalize to arbitrary N.) Another 
interesting limit obtains for Aia = — |A 2 a + Aa, Aa > if A 2 a — > — oo. In this case the 
nonlinear constraint reads (tr v^V) 2 = tr(y?V) 2 which implies for N = 2 that det (p = 0. 
This is a matrix model for singular complex 2x2 matrices. 

One can also interpret our model as the coupled system of two SU (2)-doublets for the 
weak interaction Higgs sector. This is simply done by decomposing the matrix ip a f, into 
two two-component complex fundamental representations of one of the SU (2) subgroups, 
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<£ab — > (fiib, V 9 2b- The present model corresponds to a particular left-right symmetric model 
with interactions specified by 

p = (p\(pi + ip\ip 2 2 

t = (ip\ Vl - + 4 (ip\(p 2 ) (y4v?i) • 

We observe that for a typical weak interaction symmetry breaking pattern the expectation 
values of <pi and <p 2 should be aligned in the same direction or one of them should vanish. In 
the present model this corresponds to the choice A 2A < 0. The phase structure of a related 
model without the term ~ (ip\ip2)(<f2fi) has been investigated previously [108] and shows 
second or first order transitions 21 . Combining these results with the outcome of this work 
leads already to a detailed qualitative overview over the phase pattern in a more general 
setting with three independent couplings for the quartic invariants (vtvi + v4^2) 2 , (vtvi ~~ 
^4<^ 2 ) 2 and (<pl<p 2 ) (ylsVi)- We also note that the special case A 2 a = 2Aia corresponds to 
two Heisenberg models interacting only by a term sensitive to the alignment between tpi 
and </?2, i.e. a quartic interaction of the form ((/?{(/?i) 2 + ((/9 2 y?2) 2 + 2(</?j<p 2 )(</4 < / ; 'i)- 

The model is now completely specified and it remains to extract the flow equations 
for Uk and Z k . Before this is carried out in sections 4.4 and 4.5 we present an overview 
over the phase structure in the next section. These results are obtained from a numerical 
solution of the evolution equations. 

4.3 Phase structure 

In this section we consider the U{2) x U (2) symmetric model in three space dimensions. 
The aim is to give an overview of our results concerning the phase structure and the 
effective average potential. We concentrate here on the spontaneous symmetry breaking 
with a residual U(2) symmetry group. This symmetry breaking can be observed for a 
configuration which is proportional to the identity and with (4.3) one finds r = 0. In this 
case we shall use an expansion of Uk(p, r) around r = keeping only the linear term in r. 
This amounts to assuming 

^(p,r = 0) = for n>2. (4.8) 

We will motivate this truncation in section 4.5 where we present a more detailed analysis. 
We make no expansion of t4(p, r) in terms of p. This allows the description of a first order 
phase transition where a second local minimum of Uk(p) = Ut(p,r = 0) appears. The 
p-dependence also gives information about the equation of state of the system. 

For the considered symmetry breaking pattern the short distance potential U\ given in 
eq. (4.6) is parametrized by positive quartic couplings, 

Aia, A 2A > (4.9) 



21 First order phase transitions and coarse graining have also been discussed in a multi-scalar model with 
Zi symmetry [109]. 
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and the location of its minimum is given by 

Poa = £aAia- (4-10) 

To study the phase structure of the model we integrate the flow equation for the effective 
average potential U k (cf. sect. 4.4, 4.5) for a variety of initial conditions poa, Aia and A 2 a- 
In particular, for general Aia, A 2 a > we are able to find a critical value poA = Poc for 
which the system exhibits a first order phase transition. In this case the evolution of U k 
leads at some scale fc 2 < A to the appearance of a second local minimum at the origin of 
the effective average potential and both minima become degenerate in the limit k — > 0. 
If po(k) > denotes the /c-dependent outer minimum of the potential (U' k (po) = 0, where 
the prime on denotes the derivative with respect to p at fixed k) at a first order phase 
transition one has 

lim(C4(0) - U k (po)) = 0. (4.11) 

fc^O 

A measure of the distance from the phase transition is the difference 5k,a = (poa — Poc)/A. 
If p\ and therefore poA is interpreted as a function of temperature, the deviation is 
proportional to the deviation from the critical temperature T c , i.e. 5k\ = A(T)(T C — T) 
with A(T C ) > 0. 

We consider in the following the effective average potential Uk for a nonzero scale k. This 
allows to observe the nonconvex part of the potential (cf. sect. 4.8). As an example we show 
in fig. 6 the effective average potential U k=kf for Aia = Aia/A = 0.1 and A 2 a = A 2 a/A = 2 as 
a function of the renormalized field ipn = (p^/2) 1 / 2 with p R = Z k=kf p. The scale kf is some 
characteristic scale below which the location of the minimum po{k) becomes essentially 
independent of k. Its precise definition is given below. We have normalized U kf and ip R to 
powers of the renormalized minimum <po«(^/) = (por(^/)/2) 1 ^ 2 with pm{kf) = Z kf p (kf). 
The potential is shown for various values of deviations from the critical temperature or 
Sk\. For the given examples 5k\ = —0.03, —0.015 the minimum at the origin becomes 
the absolute minimum and the system is in the symmetric (disordered) phase. Here <po R 
denotes the minimum in the metastable ordered phase. In contrast, for 5k\ = 0.04, 0.1 the 
absolute minimum is located at <Pr/<Por = 1 which characterizes the spontaneously broken 
phase. For large enough Ska the local minimum at the origin vanishes. For 5k a = the 
two distinct minima are degenerate in height 22 . As a consequence the order parameter 
makes a discontinuous jump at the phase transition which characterizes the transition 
to be first order. It is instructive to consider some characteristic values of the effective 
average potential. In fig. 7 we consider for Aia = 0.1, A 2 a = 2 the value of the renormalized 
minimum p 0R (kf) and the radial mass term as a function of — 5k a or temperature. In the 
spontaneously broken phase the renormalized radial mass squared is given by (cf. section 
4.4) 

m R (k f ) = 2Z kf 1 p U' i ; f (po), (4.12) 

22 We note that the critical temperature is determined by condition (4.11) in the limit k — » 0. Nevertheless 
for the employed nonvanishing scale k — kf the minima of Uk become almost degenerate at the critical 
temperature. 
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Figure 6: The effective average potential Uk=k t as a function of the renormalized field ip^. The potential 
is shown for various values of 5n A ~ T c — T. The parameters for the short distance potential U A are (1) 
<5k a = -0.03, (2) 5n A = -0.015, (3) 5n A = 0, (4) 5n A = 0.04, (5) 8n A = 0.1 and A iA = 0.1, A 2A = 2. 

in the symmetric phase the renormalized mass term reads 

ml R (k f ) = Z h 1 U' kf (0). (4.13) 

At the critical temperature (Ska = 0) one observes the discontinuity Ap 0R = poR.(kf) and 
the jump in the mass term Am R = m R (kf) — moji(kf) = m c R — itlq R . (Here the index 
"c" denotes 5k a = 0). The ratio Ap 0R /A is a rough measure for the "strength" of the 
first order transition. For Apofi/A 1 the phase transition is weak in the sense that 
typical masses are small compared to A. In consequence, the long-wavelength fluctuations 
play a dominant role and the system exhibits universal behavior, i.e. it becomes largely 
independent of the details at the short distance scale A -1 . We will discuss the universal 
behavior in more detail below. 

In order to characterize the strength of the phase transition for arbitrary positive values 
of Aia and A 2 a we consider lines of constant Ap 0R /A in the Aia, A 2 a plane. In fig. 8 this is 
done for the logarithms of these quantities. For fixed A 2 a one observes that the discontinuity 
at the phase transition weakens with increased Aia- O n the other hand for given Aia one 
finds a larger jump in the order parameter for increased A 2 a- This is true up to a saturation 
point where Ap 0R / A becomes independent of A 2 a- In the plot this can be observed from 
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-0.14 -0.12 -0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 

-Ska 

Figure 7: The minimum pofl and the mass term m R in units of the momentum scale A as a function of 
— 5k\ or temperature (Aia = 0.1, A 2 a = 2,k = kj). For Ska = one observes the jump in the renormalized 
order parameter Ap _R and mass Ato«. 

the vertical part of the line of constant h^Apoij/A). This phenomenon occurs for arbitrary 
nonvanishing Ap 0R /A in the strong A2A coupling limit and is discussed in section 4.6. 

In the following we give a detailed quantitative description of the first order phase 
transitions and a separation in weak and strong transitions. We consider some charac- 
teristic quantities for the effective average potential in dependence on the short distance 
parameters Aia an d A2A f° r Ska = 0. We consider the discontinuity in the renormalized 
order parameter Ap 0R and the inverse correlation lengths (mass terms) m c R and m,Q R in 
the ordered and the disordered phase respectively. Fig. 9 shows the logarithm of Ap 0R in 
units of A as a function of the logarithm of the initial coupling A 2 a- We have connected 
the calculated values obtained for various fixed Aia = 0.1, 2 and Aia — 4 by straight lines. 
The values are listed in table 4. For A 2 a/Aia ^ 1 the curves show constant positive slope. 
In this range Apo R follows a power law behavior 

Ap 0R = R(\ 2A ) 9 , 6 = 1.93. (4.14) 

The critical exponent 9 is obtained from the slope of the curve in fig. 9 for A 2 a/Aia <C 1. 
The exponent is universal and, therefore, does not depend on the specific value for Aia- On 
the other hand, the amplitude R grows with decreasing Aia- For vanishing A 2 a the order 
parameter changes continuously at the transition point and one observes a second order 
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Figure 8: Lines of constant jump of the renormalized order parameter A/?o,r/A at the phase transition 
in the ln(AiA), ln(A2A) plane. The curves correspond to (1) \n(Ap R/ A) = —4.0, (2) ln(Ap it/A) = —4.4, 
(3) ln(A/WA) = -10.2, (4) ln(Ap ij/A) = -14.3. 



phase transition as expected for the 0(8) symmetric vector model. As A2A/A1A becomes 
larger than one the curves deviate substantially from the linear behavior. The deviation 
depends on the specific choice of the short distance potential. For A 2A /Ai A 3> 1 the curves 
flatten. In this range Ap 0R becomes insensitive to a variation of the quartic coupling A 2 a- 

In addition to the jump in the order parameter we present the mass terms m c R and 
which we normalize to Apo R . In fig- 10 these ratios are plotted versus the logarithm of the 
ratio of the initial quartic couplings A 2 a/Aia- Again values obtained for fixed Aia = 0.1, 2 
and Aia = 4 are connected by straight lines. The universal range is set by the condition 
m c R / Ap 0R ~ const (equivalently for rriQ R / Ap 0R ) . The universal ratios are m c R / Ap 0R = 1.69 
and rriQ R / Ap 0R = 1.26 as can be seen from table 4. For the given curves universality holds 
approximately for A 2 a/Aia^1/2 and becomes "exact" in the limit A 2 a/Aia — > 0. In this 
range we obtain 

m R = S(\ 2A ) e , m c 0R = S(\ 2A ) e . (4.15) 

The universal features of the system are not restricted to the weak coupling region of 
A 2A - This is demonstrated in fig. 10 for values up to A 2A — 2. The ratios m c R /Ap 0R and 
m,Q R /Apo R deviate from the universal values as A 2 a/Aia is increased. For fixed A 2 a a larger 
Aia results in a weaker transition concerning Ap 0R /A. The ratio m c R / Ap 0R increases with 
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Figure 9: The logarithm of the discontinuity of the renormalized order parameter A/?o.r/A as a function 
of ln(A2A). Data points for fixed (1) Aia = 0.1, (2) Aia = 2, (3) Aia = 4 are connected by straight lines. 

Aia for small fixed A 2 a whereas in the asymptotic region, A 2 a/Aia 1, one observes from 
fig. 10 that this tendency is reversed and m c R / Ap 0R , m,Q R /Ap 0R start to decrease at about 
A 1A ~2. 

In summary, the above results show that though the short distance potential U\ indi- 
cates a second order phase transition, the transition becomes first order once fluctuations 
are taken into account. This fluctuation induced first order phase transition is known in 
four dimensions as the Coleman- Weinberg phenomenon [70] . Previous studies of the three 
dimensional U(2) x U(2) symmetric model using the e-expansion [110, 31] already indi- 
cated that the phase transition should be a fluctuation induced first order transition. The 
validity of the e-expansion for weak first order transitions is, however, not clear a priori 
since the expansion parameter is not small - there are known cases where it fails to predict 
correctly the order of the transition [111]. The question of the order of the phase tran- 
sition has been addressed also in lattice studies [112] and in high-temperature expansion 
[113]. All studies are consistent with the first order nature of the transition and with the 
absence of nonperturbative infrared stable fixed points. It is, however, notoriously difficult 
to distinguish by these methods between weak first order and second order transitions. 
Our method gives here a clear and unambiguous answer and allows a detailed quantitative 
description of the phase transition. The universal form of the equation of state for weak 
first order phase transitions is presented in section 4.7. 
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Table 4: The discontinuity in the renormalized order parameter Apo R and the critical inverse correlation 
lengths m c R and m,Q R in the ordered and the disordered phase respectively. For small A 2 a/Aia the ratios 
m R j Apo_R and vn^ R l Apofi become universal. 



In the following we specify the scale kf for which we have given the effective average 
potential U^. We observe that [/& depends strongly on the infrared cutoff k as long as k is 
larger than the scale /C2 where the second minimum of the potential appears. Below /c 2 the 
two minima start to become almost degenerate for T near T c and the running of po(k) stops 
rather soon. The nonvanishing value of /c 2 induces a physical infrared cutoff and represents 
a characteristic scale for the first order phase transition. We stop the integration of the 
flow equation for the effective average potential at a scale kf < k2 which is determined 
in terms of the curvature (mass term) at the top of the potential barrier that separates 
the two local minima of £4 at the origin and at po(k). The top of the potential barrier at 
ps(/c) is determined by 

U' k (p B ) = (4.16) 
for < psik) < po(k) and for the renormalized mass term at ps(/c) one obtains 



m% R (k) = 2Z^p B Ul{p B ) < 0. 



(4.17) 
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Figure 10: The inverse correlation lengths m c R and m,Q R in the ordered and the disordered phase respec- 
tively. They are normalized to Apofl and given as a function of ui^a/Aia)- Data points for fixed (1) 
Aia = 0.1, (2) Aia = 2, (3) Aia = 4 are connected by straight lines. 



We fix our final value for the running by 

kf-\m% R {k f )\ 

For this choice the coarse grained effective potential U^ f essentially includes all fluctua- 
tions with momenta larger than the mass \m,B ,r\ at the top of the potential barrier. It is a 
nonconvex function which is the appropriate quantity for the study of physical processes 
such as tunneling or inflation. The nonconvex part of Uk is considered in detail in section 
4.8. There we will also discuss that the appropriate choice for a coarse graining scale k is 
often far from obvious. 



(4.18) 



4.4 Scale dependence of the effective average potential 

In this section we derive the flow equation for the effective average potential Uk- We point 
out a related investigation of the four dimensional SU(N) x SU (N) symmetric linear sigma 
model coupled to fermions that has been studied previously within the framework of the 
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effective average action [46]. There the flow equations for a polynomial approximation of 
Uk can be found for arbitrary dimension d. To study the equation of state of the three 
dimensional theory we keep here the most general form for U k . 

We define U k by evaluating the average action for a constant field with T k = VtU k 
where Vt denotes the volume. To evaluate the r.h.s. of (2.10) with the ansatz (4.2) we 
expand Y k around a constant background configuration. With the help of the U{2) x U (2) 
transformations the matrix field tp can be turned into a standard diagonal form with 
real nonnegative eigenvalues. Without loss of generality the evolution equation for the 
effective potential can therefore be obtained by calculating the trace in (2.10) for small field 
fluctuations Xab around a constant background configuration which is real and diagonal, 

<Pab = VaSab , V*a = <Pa- (4.19) 

We separate the fluctuation field into its real and imaginary part, Xab = -^(XRab+iXiab) and 
perform the second functional derivatives of T k with respect to the eight real components. 
For the constant configuration (4.19) it turns out that has a block diagonal form 
because mixed derivatives with respect to real and imaginary parts of the field vanish. The 
remaining submatrices S 2 Fk/5xR^X C R an d 5 2 T k / Sxt^xf 1 can be diagonalized in order to 
find the inverse of rj^ + R k under the trace occuring in eq. (2.10). Here the momentum 

(2) 

independent part of rjj. defines the mass matrix by the second functional derivatives of 
Uk- The eight eigenvalues of the mass matrix are 

(Mt? = U' k + 2(p±(p 2 -ry/ 2 )d T U k , 
{Mff = U' k ±2r^d T U k 

corresponding to second derivatives with respect to xi an d 

{Mff = (Mt) 2 , 

{Mff = U' k + pU£ + 2pd T U k + ATd T U' k + ApTd 2 T U k 

± {r {U' k ' + 4d T U k + Apd T U' k + 4rd 2 T U k ) 2 (421) 

+ (P 2 " t) (U> k > - 2d T U k - ^rd 2 U k f} 1/2 

corresponding to second derivatives with respect to xr- Here the eigenvalues are expressed 
in terms of the invariants p and r using 

V\ = \(P + ^I% p\ = \{p-r^ 2 ) (4.22) 

and we adopt the convention that a prime on U k (p, r) denotes the derivative with respect 
to p at fixed r and k and d™U k = d n U k /(dr) n . 

The flow equation for the effective average potential is simply expressed in terms of the 



(4.20) 
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mass eigenvalues 



1 f d^q 

m Uk ^ T) = 2lwYdi Rk{q) 

f 2 



+ 57^ I TJI^Z ^2 + 57^ I 7~h7+7Z 3vi ( 4 - 2 3) 



\P k (q) + (M+(p,T)y P k (q) + (Mr(p,r)y P k (q) + (M+(p,r)) 2 
1 1 1 

+ ^ ; ; . „, , . 7TTX7 + 



P fc (g) + (M 2 -(p,r)) 2 P fc (g) + (M+(p,r)) 2 P fc (g) + (M 4 -(p, r)) 2 J ' 

On the right hand side of the evolution equation appears the (massless) inverse average 
propagator 

P k (q) = Z k q 2 + R k (q) = - _^ 2/fc2 (4.24) 

which incorporates the infrared cutoff function R k given by eq. (2.3). The only approxima- 
tion so far is due to the derivative expansion (4.2) of T k which enters into the flow equation 
(4.23) through the form of P k . The mass eigenvalues (4.20) and (4.21) appearing in the 
above flow equation are exact since we have kept for the potential the most general form 
U k (p,r). 



Spontaneous symmetry breaking and mass spectra 

In the following we consider spontaneous symmetry breaking patterns and the corre- 
sponding mass spectra for a few special cases. For the origin at ip a b = all eigenvalues equal 
U k (0,0). If the origin is the absolute minimum of the potential we are in the symmetric 
regime where all excitations have mass squared U' k (0,0). 

Spontaneous symmetry breaking to the diagonal U(2) subgroup of Z7" (2) x U(2) can 
be observed for a field configuration which is proportional to the identity matrix, i.e. 
(fab = tfSab- The invariants (4.3) take on values p = 2ip 2 and r = 0. The relevant 
information for this symmetry breaking pattern is contained in U k (p) = U k (p,r = 0). In 
case of spontaneous symmetry breaking there is a nonvanishing value for the minimum 
Po of the potential. With U' k (po) = one finds the expected four massless Goldstone 
bosons with (Mf) 2 = (M^) 2 = (M^) 2 = 0. In addition there are three massive scalars 
in the adjoint representation of the unbroken diagonal SU(2) with mass squared (M^) 2 = 
(M 3 + ) 2 = (M^) 2 = Ap d T U k and one singlet with mass squared (M 4 + ) 2 = 2p U k . The 
situation corresponds to chiral symmetry breaking in two flavor QCD in absence of quark 
masses and the chiral anomaly. The Goldstone modes are the pseudoscalar pions and the 
i] (or T)'), the scalar triplet has the quantum numbers of a and the singlet is the so-called 
cr-field. 

Another interesting case is the spontaneous symmetry breaking down to a residual 
U(l) x f/(l) x U(l) subgroup of U(2) x U{2) which can be observed for the configuration 
¥>ab — f^ai^ab (p = ^ \ t = ip A = p 2 ) . Corresponding to the number of broken generators 
one observes the five massless Goldstone bosons (Mj 4 ) 2 = (M^) 2 = ) 2 = for the 
minimum of the potential at U' k + 2p d T U k = 0. In addition there are two scalars with mass 
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squared (M 2 ") 2 = (M 4 ~) 2 = U' k - 2p d T U k and one with (M 4 + ) 2 = U' k + 2p U' k ' + 6p d T U k + 
8p 2 d T U k + 8p 3 d 2 U k . 

We finally point out the special case where the potential is independent of the second 
invariant r. In this case there is an enhanced 0(8) symmetry instead of U(2)xU (2). With 
<9™£4 = and U' k (p ) = one observes the expected seven massless Goldstone bosons and 
one massive mode with mass squared 2p U k . 



Scaling form of the flow equation 

For the 0(8) symmetric model in the limit A 2 a = one expects a region of the parameter 
space which is characterized by renormalized masses much smaller than the ultraviolet 
cutoff or inverse microscopic length scale of the theory. In particular, in the absence of 
a mass scale one expects a scaling behavior of the effective average potential U k . The 
behavior of U k at or near a second order phase transition is most conveniently studied 
using the scaling form of the evolution equation. This form is also appropriate for an 
investigation that has to deal with weak first order phase transitions as encountered in the 
present model for A 2 a > 0. The remaining part of this section is devoted to the derivation 
of the scaling form (4.33) of the flow equation (4.23). 

In the present form of eq. (4.23) the r.h.s. shows an explicit dependence on the scale k 
once the momentum integration is performed. By a proper choice of variables we cast the 
evolution equation into a form where the scale no longer appears explicitly. We introduce a 
dimensionless potential u k = k~ d U k and express it in terms of dimensionless renormalized 
fields 

p = Z k k 2 ~ d p, 

(4.25) 

f = Z 2 k 4 ~ 2d r. 1 ' 

The derivatives of u k are given by 

<9?4 m) (P, f) = Z k - 2n - m kV n+m -V d - 4n - 2m d?uj : m) (p, r) . (4.26) 

(Note that u k m ^ denotes m derivatives with respect to p at fixed f and k, while U k m ^ denotes 
m derivatives with respect to p at fixed r and k). With 

d 

;Uk(P,T)\p,f = -du k (p,T) + (d-2 + r])pu' k (p,T) + (2d-A + 27])rdrU k (p,T) 



dt 

+k' d ^-U k (p(p),r(f))^ T 



(4.27) 



dt 

one obtains from (4.23) the evolution equation for the dimensionless potential. Here the 

anomalous dimension r\ arises from the t-derivative acting on Z k and is given by eq. (4.5). 
It is convenient to introduce dimensionless integrals by 

1 1 fk'-wP^ - [*"> - *)] < 4 ' 28 > 

where 

v-, 1 = 2 d+1 ^T ( i) . (4.29) 
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The explicit form of l d and l d reads 



tfM = -/ ° O ^ +1 ^b(l + r( 2/ )) + c]- 1 , 



1 ;-oo d 

*oM = o/ rfyy2 r ( y )[ y (i + r ( y )) + a; ]- 

2 jo 



(4.30) 



with the dimensionless infrared cutoff function 



r(y) = (4.31) 

The behavior of these functions is discussed in section 4.5 where we also consider their 
derivatives with respect to uj. Using the notation 

l*(u>;r,) = l d (u) - (4.32) 

one obtains the flow equation for the dimensionless potential 

d 

— u k (p, r) = -du k (p, r) + (d-2 + rj)pu' k (p, f) + (2d - 4 + 2r])rd f u k (p, r) 

+4^ ((m+(p, f )) 2 ; 77) + 4^ d ((m7(p, f )) 2 ; 17) + 2v d l d ((m+(p, f )) 2 ; 77) (4-33) 
+2^/ d (K (p, f )) 2 ; 77) + 2v d l d ((m+(p, f )) 2 ; 77) + 2^ ((m7 (^r)) 2 ;^) 

where the dimensionless mass terms are related to (4.21) according to 

(mf(p, r)) 2 = ^ — >-. (4.34) 

Eq. (4.33) is the scaling form of the flow equation we are looking for. For a f-independent 
potential it reduces to the evolution equation for the 0(8) symmetric model [92, 93]. The 
potential u k at a second order phase transition is given by a /c-independent (scaling) solu- 
tion duk/dt = [92, 93]. For this solution all the /c-dependent functions on the r.h.s. of 
eq. (4.33) become independent of k. For a weak first order phase transition these functions 
will show a weak /c-dependence for k larger than the inherent mass scale of the system (cf. 
section 4.6). There is no particular advantage of the scaling form of the flow equation for 
strong first order phase transitions. 



4.5 Solving the flow equation 

Eq. (4.33) describes the scale dependence of the effective average potential u k by a nonlinear 
partial differential equation for the three variables t, p and f . A particular difficulty for its 
analytical study is that the integral Iq(u; rf) (cf. eq. (4.30)) can be done analytically only for 
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certain limits of the arguments. The complicated form of the equation therefore suggests 
a numerical solution. We will use a method that relies on a simultaneous expansion of 
the potential around a number of field values pi,fj, i = 1, j = 1, . . . , /' for given 
numbers /, /'. The expansions around different points are matched to obtain the general 
field dependence of the potential. As a consequence we cast the partial differential equation 
(4.33) into a system of ordinary differential equations. The method is developed in [63] for 
a computation of the critical equation of state for O(N) symmetric Heisenberg models [62] 
(cf. section 3). The generalization to the present model is described below. This approach 
has some favorable aspects. The main advantage is that it allows a very efficient integration 
of the differential equations using a standard Runge-Kutta algorithm without the occurance 
of numerical instabilities 23 . The coupled set of ordinary differential equations describes the 
flow of couplings defined as derivatives of the potential at given points, e.g. at the minimum 
of the potential. These couplings often allow direct physical interpretation and some of 
their properties can be read off from the analytic structure of their flow equations. We will 
exploit this fact to explain the results we obtain from the numerical solution in section 4.6. 

We concentrate in the following on spontaneous symmetry breaking with a residual U (2) 
symmetry group. As we have already pointed out in section 4.2 this symmetry breaking 
can be observed for a configuration which is proportional to the identity and we have f = 0. 
In this case the eigenvalues (4.20) and (4.21) of the mass matrix with (4.34) are given by 



and on the r.h.s. of the partial differential equation (4.33) for Uk{p) = Uk(p, f — 0) only the 
functions u' k (p), u k (p) and dfU k (p) appear. At fixed p = pi the /c-dependence of u k is then 
determined by the couplings u' k (pi), u k (pi) and dfU k (pi). We determine these couplings 
through the use of flow equations which are obtained by taking the derivative in eq. (4.33) 
with respect to p and f evaluated at p — pi, f — 0. These flow equations for u' k , u k 
and dfUk involve also higher derivatives of the potential as u' k , u k 4 \ dfu' k and dfU k . The 
procedure will be to evaluate the flow equations for u' k and u k at different points pi, f — 
for i — 1, . . . , I and to estimate the higher p-derivatives appearing on the right hand side of 
the flow equations by imposing matching conditions. The same procedure can be applied 
to dfUk and in order to obtain an equivalent matching we also consider the flow equation 
for dfu' k . One could proceed in a similar way for dfuk which appears on the right hand 
side of the evolution equation of dfU k and so on. However, since we are interested in the 

23 Thcse instabilities arise through the inevitable small round-off errors at every integration step. The 
standard discretised version of eq. (4.33), using finite difference quotients for the derivatives, becomes 
strongly oscillating after a few integration steps. Already for the most simple linear partial differential 
equations one finds that the resulting error may grow exponentially with the number of steps (see e.g. rcf. 
[114] for an analytically solved example). A possible finite difference scheme uses "backward" difference 
quotients as we have exemplified in ref. [63]. 





u' k + 4pd f U k , 



(4.35) 
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p-dependence of the potential at r = we shall use a truncated expansion 24 in r with 

dy k (p, f = 0) = for n > 2. (4.36) 

In three space dimensions the neglected (p-dependent) operators have negative canonical 
mass dimension. We make no expansion in terms of p since the general p-dependence 
allows a description of a first order phase transition where a second local minimum of 
Uk{p) appears. The approximation (4.36) only affects the flow equations for dfU k and dfu' k 
(cf. eqs. (4.43) and (4.44)). The form of the flow equations for u' k and u k is not affected 
by the truncation (cf. eqs. (4.41) and (4.42)). From u' k we obtain the effective average 
potential u k by simple integration. We have tested the sensitivity of our results for u' k 
to a change in d f u k by neglecting the p-dependence of the f-derivative. We observed no 
qualitative change of the results. We expect that the main truncation error is due to the 
derivative expansion (4.2) for the effective average action. 
To simplify notation we introduce 

e (p) = u' k {p,f = 0), 
X 1 (p) = <(p,f = 0), (4.37) 
A 2 (p) = 4dfU k (p,T = Q). 

Higher derivatives are denoted by primes on the p-dependent quartic "couplings", i.e. 
X[ = u' k) X' 2 = dfu' k etc. It is convenient to introduce functions 77) that can be related 
to Zq(cj; 77) (4.30) by differentiation with respect to the mass argument: 

l*(u>;r,) = lf(uj) - riif(u) 

= -j^lifar,), (4.38) 
1 d 

= ~-Q^ l n( UJ > r l) fOT U ^ 1 - 

The explicit form of I* and I* is given in eq. (C.l) in the appendix. For uo > — 1 they 
are positive, monotonically decreasing functions of uo. In leading order 1^ and 1% vanish 
~ ^-( n+1 ) for arguments uj ^> 1. They introduce a "threshold" behavior that accounts 
for the decoupling of modes with mass squared larger than the infrared cutoff Z k k 2 . For 
vanishing argument they are of order unity. As w — > —1 the functions 1%, 1% exhibit a 
pole for d < 2(n + 1). The pole structure is discussed in the appendix. We also use two- 
parameter functions ^ i n2 (^i, ^2', v) [46]. For n\ — n<i — 1 their relation to the functions 
Z^(a>; rf) can be expressed as 

lf 1 (uj 1 ,uj2;i]) = hf^r)) - l((uj 2 ;r)j\ for W!^U2, 



0J 2 — OJi 

lf 1 (uj,uj;r)) = li(u;rj) 



(4.39) 



24 In principle one could also consider points with f 7^ in the neighborhood of f = and use the 
additional information to estimate the higher f-derivatives as it is done for the higher p-derivatives. 
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and 



C+i,n 2 (^i>^2;?7) = 



1 d 



CnMi^v), t fl2 (wi,w 2 ;i) =li atni (u 2 ,u 1 ;ri). (4.40) 



With the help of these functions the scale dependence of e is described by 
de 



dt 



(-2 + r])e + (d - 2 + r^pAi - 6v d (Ai + A 2 + pA' 2 )Z?(e + pA 2 ; 77) 



-2u d (3Ai + 2pA;)Z?(e + 2pA i; 77) - 8u d AiZ?(e; 77) 



(4.41) 



and for Ai one finds 
9Ai 



at 



= (d-4 + 277)A 1 + (d-2 + 77)/3A' 1 



(Ax + A 2 + p\' 2 ) 2 l d 2 (e + pA 2 ; 77) - (A; + 2A' 2 + p%)l*(e + pA 2 ; 77) 
+2v d [(3Ai + 2/5A / 1 ) 2 Z 2 z (e + 2pA i; 77) - (5A' X + 2p)![)l d {e + 2pA i; 77) 
+87;, [(A 1 ) 2 Z 2 d ( e ;77)-A / 1 Zf( e ;77y. 

Similarly the scale dependence of A 2 is given by 
<9A 2 



(4.42) 



dt 



(d - 4 + 27/)A 2 + (d - 2 + r])p\' 2 - 4v d (\ 2 yi d :1 (e + pA 2 , e; 7/) 



+2t; d [3(A 2 ) 2 + 12AiA 2 + 8pA' 2 (Ai + A 2 ) + 4p 2 (A' 2 ) 
-Uv d X' 2 l d (e + pA 2 ; 77) - 27; d (5A' 2 + 2p\' 2 ')l d (e + 2pX 1 ; 77) 
+2t;, [(A 2 ) 2 Z^(e;77)-4A 2 Zf(e;77) 



Z( 1 (e + pA 2 ,e + 2pA 1 ;77) 



(4.43) 



and for A 2 it reads 
d\' 2 



dt 



(2d - 6 + 377) A' 2 + (d - 2 + rj)p\' 2 ' 



+47; d (A 2 )^ 



(Ai + A 2 + p\' 2 )lUe + pA 2 , e; 77) + AiZf 2 (e + pA 2 , e; 77) 



-27;„ 



3A 2 (4Ai + A 2 ) + 4pA' 2 (2Ai + 2A 2 + pA'2) 



(Ai + A 2 + pA' 2 ) 



I 2 d i(e + pA 2 , e + 2pAi; 77) + (3A X + 2pX[)l d 2 (e + pA 2 , e + 2pA x ; 77) 



A' 2 (7A 2 + 10Ai) + 6A 2 A; + 4p( A 2 (Ai + A 2 + pA' 2 ) + A' 2 (2A 2 + A'^ 



+4v d 

lf A (e + pA 2 , e + 2pAi; 77) - 8udA 2 A' 2 Zf + pA 2 , e; 77) 
+2v d (3X 1 + 2pA' 1 )(5A , 2 + 2p\' 2 r )l d (e + 2pA i; 77) + 14u d A' 2 (Ai + A 2 + pA' 2 ) 
l d 2 {e + pA 2 ; 77) - 2v d {lX 2 ' + 2pA' 2 ")Z? (e + 2pA l5 v) ~ Uv A X^lf(e + pA 2 ; 77) 
-4t;,(A 2 ) 2 A 1 Zf (e; 77) + 4u d A' 2 (2Ai + A 2 )Z 2 d (e; 77) - 8u d A£Z?(e; 77). 



(4.44) 
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We evaluate the above flow equations at different points pi for i = and use a 

set of matching conditions that has been proposed in ref. [63]. The generalization of 
these conditions to the present model is obtained by considering fourth order polynomial 
expansions of e(p) and A 2 (p) around some arbitrary point p iy 



(4.45) 



(e)i(p) = e i + X 1>i (p-p i ) + ^ i (p-p i ) 2 + ^X'l i (p-p i f, 

with €i = e(pi), A 2 ,j = A 2 (pi) etc. Using similar expressions for 

(AiMp) = §z(eW) , (A' 2 ),(p) = §z(X2)i(p) (4.46) 

the matching is done by imposing continuity at half distance between neighboring expan- 
sion points, 

( £ ),(^±l) =(e ), +1 (^±l) , (Al)( = ( A ,), +1 (^±1 

(A 2 ), (^P 1 ) = (A 2 ),+, , (A' 2 ) i f^±l N )=(A' 2 ), + /« + «« 



2 

for i = 1, — 1 and 



(4.47) 



= (AiW f^±i) . (AS), (&±£±0 = (A?) j+1 ^ + ^« 



2 

(4.48) 

for the initial and end points, j = I and j — I — 1. Together these 4(Z — 1) conditions 
(4.47) for all I — 1 > 2 intermediate points and the four conditions (4.48) make up two 
independent algebraic systems of each 21 equations. From the first set of equations one 
obtains a unique solution for A' 1;i and A" j. The second set is identical in structure and A^, 
A 2 'j can be obtained from the solutions for A' 1;i and A" f with the substitutions €j — > A 2j -, 
Ai,j — > A' 2 j, A' X j — > A 2 j, A'/j — > A 2 j for j = 1,...,/. With the help of these algebraic 
solutions we eliminate A'^p), X'l(p), A 2 (p) and A 2 '(p) in the flow equations (4.41) - (4.44) 
for all / points pV 25 Therefore, equations (4.41) - (4.44) are turned into a closed system of 
Al ordinary differential equations for the unknowns e(p~j), Ai(p"i), A 2 (p~j) and A' 2 (pj). 

25 The algebraic solutions and A", (equivalently A^ and A^'J do incorporate information from the 
whole range of points pj with j = 1, . . . , I. It is a feature of the matching conditions (4.47), (4.48) that 
the contributions from points Pj^i to A^ i , A'/j rapidly decrease with increasing \i — j\. For equal spacings 
between neighboring expansion points contributions from points pj with j > i + 1 (j < i — I) arc typically 
suppressed by a factor ^,10~\ l ~^ +1 as compared to the contribution from the nearest neighbor point p i+i 
(Pi-i)- As a consequence solutions A^ i , X'{ i for inner points with 1 < i € 1 become independent from 
boundary points. We observe approximate translational invariance for inner point solutions, i.e. A^ j±n 
and A' 1 ' i±n are approximately obtained from the solutions \' 1A and X'{ A with the substitutions ej — > £j±„, 
Ai,j — > Xi.j± n for 1 < j < I if i and i ± n are sufficiently far away from the boundaries. The decoupling 
from distant points and the translational invariance for inner points can be used to obtain approximate 
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If there is a minimum of the potential at nonvanishing k = fio we use expansion points 
that are proportional to the minimum, i.e. pi = —k with i — 1, . . . , I and fixed integer n. 
The condition e(/t) = can be used to obtain the scale dependence of n{k): 



dn \i-i de \ 
¥ = - [Al(K)] di Vp - 

= -(d-2 + V )K + 6v d (l + A ^+^) ) if {kX2{k) . v) (449) 



+ 2v d ^3 + If (2«Ai(«); 77) + (0; 77). 

To make contact with /5-functions for the couplings at the potential minimum k we point 
out the relation 

dAg(«) gAg +A (m+l) f > 

Similar relations hold for e(pi), Xi(pi) etc., e.g. 

de(pi) de i — 1 dn 

~&r = 9t'^ + ^r Ai(Pi) dt' (451) 

dXi(pi) d\i. [~}-\t(~ \ dK 
dt ~ dt Vp=h+ n XApl) df 

We integrate the 4/ — 1 differential equations (4.41) - (4.44) for the couplings e(pi), Xi(pi), 
X 2 (pi) and X' 2 (pi) (with d/dt replaced by d/dt according to (4.51)) and the one for k (4.49) 
with a fifth-order Runge-Kutta algorithm using the embedded fourth-order method for 
precision control. The general p-dependence is recovered by patching the simultaneous 
expansions around different points at half distance between neighboring expansion points. 
The polynomial patching improves with decreasing distance which is used to check the 
stability of numerical results. In order to check for systematic errors the algorithm has 
been verified by comparison with a conventional "backward" finite difference method [63] . 

It remains to compute the anomalous dimension 77 defined in (4.5) which describes the 
scale dependence of the wave function renormalization Z^. We consider a space dependent 
distortion of the constant background field configuration (4.19) of the form 

Vab(x) = Va 5 ah + [5 V e~ iQx + 5<p*e lQx ] S ab . (4.52) 

Insertion of the above configuration into the parametrization (4.2) of Tk yields 

Z k = Z t ( fi ,r,Q> = 0) = l-L-Jjm^^^^l^, (4.53) 



expressions which become useful if a large number of expansion points is considered. We use the exact 
algebraic solution for I = 10 points. For I > 10 we apply the approximate translational invariance to 
generate from X' 15 , X' 16 additional solutions A' 15+2i , e+2i for i = 1, . . . , (I — 10)/2 with / even and 
equivalently for A" 5 , A" 6 . With X[ ; _ 3 , . . . ,X' lt and A" ; _ 3 , . . . ,X'( t from the calculation with 10 points 
one obtains the desired generalization. Wc have used runs with different choices of I in order to check the 
stability of the numerical results. 



4.6 Renormalization group flow 
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To obtain the flow equation of the wave function renormalization one expands the effective 
average action around a configuration of the form (4.52) and evaluates the r.h.s. of eq. 
(2.10). This computation has been done in ref. [46] for a "Goldstone" configuration with 



J ab 



S a i$b2 - S a2 5 bl (4.54) 



and (p a Sab = ¥&ab corresponding to a symmetry breaking pattern with residual U (2) sym- 
metry. The result of ref. [46] can be easily generalized to arbitrary fixed field values of p 
and we find 

r](k) = 4^p[4(A 1 ) 2 mJ 2 (e,e + 2pA 1 ^) + (A 2 ) 2 m^ 2 (e,e + pA 2 ;r ? )]. (4.55) 

The definition of the threshold function 

m^wi, u 2 \ rj) = m 22 (ui, uo 2 ) - rim 2 2 (ui,u 2 ) (4.56) 

can be found in appendix A. For vanishing arguments the functions m 22 and rh 22 are 
of order unity. They are symmetric with respect to their arguments and in leading order 
m 2 2 {Q, uj) ~ ^2,2(0? <jS) ~ uj~ 2 for uj >> 1. According to eq. (4.4) we use p — k to define the 
uniform wave function renormalization 



Z k = Z k {K). (4.57) 

We point out that according to our truncation of the effective average action with eq. (4.55) 
the anomalous dimension r\ is exactly zero at p = 0. This is an artefact of the truncation 
and we expect the symmetric phase to be more affected by truncation errors than the 
spontaneously broken phase. We typically observe small values for r)(k) = — d(\nZ k )/dt 
(of the order of a few per cent). The smallness of i] is crucial for our approximation 
of a uniform wave function renormalization to give quantitatively reliable results for the 
equation of state. For the universal equation of state given in sect. 4.7 one has rj = 0.022 
as given by the corresponding index of the 0(8) symmetric "vector" model. 



4.6 Renormalization group flow 

To understand the detailed picture of the phase structure presented in section 4.3 we will 
consider the flow of some characteristic quantities for the effective average potential as the 
infrared cutoff k is lowered. We will always consider in this section the trajectories for the 
critical "temperature" , i.e. 5k a = 0, and we follow the flow for different values of the short 
distance parameters Aia and A2A. The discussion for sufficiently small 5k a is analogous. 
In particular, we compare the renormalization group flow of these quantities for a weak 
and a strong first order phase transition. In some limiting cases their behavior can be 
studied analytically. For the discussion we will frequently consider the flow equations for 
the quartic "couplings" Ai(p), A 2 (p) eqs. (4.42), (4.43) and for the minimum k eq. (4.49). 
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In fig. 11, 12 we follow the flow of the dimensionless renormalized minimum k and the 
radial mass term m 2 = 2k\i(k) in comparison to their dimensionful counterparts poi? — kn 
and m 2 R = k 2 fh 2 in units of the momentum scale A. We also consider the dimensionless 
renormalized mass term fh\ = kX 2 (k) corresponding to the curvature of the potential in the 
direction of the second invariant f. The height of the potential barrier Us(k) = k 3 Uk(f>B) 
with u' k {pB) = 0, < pB < k, and the height of the outer minimum U (k) = k 3 Uk(n) 
is also displayed and will be discussed in section 4.8. Fig. 11 shows these quantities as 
a function of t = \n(k/A) for Aia = 2, A 2 a = 0.1. One observes that the flow can be 
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Fi gure 11: Scale dependence of the dimensionless renormalized masses to 2 , to 2 ., minimum k and dimen- 



sionful counterparts m 2 R — k 2 rh 2 , p R = kn in units of A. We also show f7s(fc) and Uo(k), the value of 
the potential at the top of the potential barrier and at the minimum /9 ij, respectively. The short distance 
parameters are Aia = 2, A 2 a = 0.1 and Ska = 0. 



separated into two parts. The first part ranging from t — to t ~ —6 is characterized by 
k ~ const and small fh 2 . It is instructive to consider what happens in the case fh\ = 0. 
In this case A 2 = and the flow is governed by the Wilson-Fisher fixed point of the 0(8) 
symmetric theory. At the corresponding second order phase transition the evolution of 
leads to the scaling solution of (4.33) which obtains for du^/dt = 0. As a consequence 
becomes a /c-independent function that takes on constant (fixed point) values [92, 93]. In 
particular, the minimum k of the potential takes on its fixed point value n(k) = k*. The 
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fixed point is not attractive in the U(2) x U(2) symmetric theory and A2A is an additional 
relevant parameter for the system. For small A2 the evolution is governed by an anomalous 
dimension d\ijdt = A\ 2 with A < 0, leading to the increasing fh\ as k is lowered. 

The system exhibits scaling behavior only for sufficiently small A 2 . As m| increases the 
quartic coupling Ai and therefore the radial mass term m 2 is driven to smaller values as 
can be observed from fig. 11. For nonvanishing A2 the corresponding qualitative change 
in the flow equation (4.42) for Ai is the occurance of a term ~ A|. It allows to drive Ai 
to negative values in a certain range of p < k and, therefore, to create a potential barrier 
inducing a first order phase transition. We observe from the plot that at t& — 9.5 a second 
minimum arises (Ub 7^ 0). The corresponding value of k = Ae* = k 2 sets a characteristic 
scale for the first order phase transition. Below this scale the dimensionless, renormalized 
quantities approximately scale according to their canonical dimension. The dimensionful 
quantities like poi? or m 2 R show only a weak scale dependence in this range. In contrast to 
the above example of a weak first order phase transition with characteristic renormalized 
masses much smaller than A fig. 12 shows the flow of the corresponding quantities for a 
strong first order phase transition. The short distance parameters employed are Aia = 0.1, 
A2A = 2. Here the range with k ~ const is almost absent and one observes no approximate 
scaling behavior. 




-2 -1.5 -1 -0.5 



t = \n(k/A) 

Figure 12: Same as fig. 11 for Aia = 0.1 and A 2 a = 2. 
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With the above examples it becomes easy to understand the phase structure presented 
in section 4.3. For the given curves of figs. 9 and 10 we distinguished between the range 
^2a/Aia <C 1 and A2A/A1A 3> 1 to denote the weak and the strong first order region. For 
^2a/Aia <C 1 the initial renormalization group flow is dominated by the Wilson-Fisher 
fixed point of the 0(8) symmetric theory. In this range the irrelevant couplings are driven 
close to the fixed point for some "time" \t\ = —ln(k/A), loosing their memory on the 
initial conditions given by the short distance potential u\. As a consequence we are able 
to observe universal behavior as is demonstrated in fig. 10. 

To discuss the case A 2 a/Aia > 1 we consider the flow equations for the couplings at 
the minimum k ^ of the potential given by (4.49) and (4.50) with (4.42), (4.43). In the 
limit of an infinite mass term m| = k\2{k) — > oo the /^-functions for \\{k) and k become 
independent from \2(k) due to the threshold functions, with ^(^2) ~ (/tA2)~ (n+1 - ) for 
large k\2(k). As a consequence (5^ and f3 K equal the /^-functions for an 0(5) symmetric 
model. We argue in the following that in this large coupling limit fluctuations of massless 
Goldstone bosons lead to an attractive fixed point for A 2 («). We take the flow equation 
(4.50), (4.43) for A2(/c) keeping only terms with positive canonical mass dimension for a 
qualitative discussion. (This amounts to the approximation \^\k) = A 2 n ^(ft) = for 
n > 1.) To be explicit, one may consider the case for given Aia = 2. The critical cutoff 
value for the potential minimum is ka — 0.2 for A 2 a 1. For kA 2 (k) 1 and taking 
i] ~ the /^-function for A 2 (/t) is to a good approximation given by (d = 3) 

dA 2( K ) _ X , 0,. 1 \ ^2,3 



dt 



-\ 2 (K)+2v 3 (\ 2 (K)yi*(0). (4.58) 



The second term on the r.h.s. of eq. (4.58) is due to massless Goldstone modes which 
give the dominant contribution in the considered range. The solution of (4.58) implies an 
attractive fixed point for A 2 (/t) with a value 

Starting from A2A one finds for the "time" \t\ necessary to reach a given A 2 (/t) > A2*(k) 

2A 



A 2 («) ( 1 - 



This converges to a finite value for A2A - ► 00. The further evolution therefore becomes 
insensitive to the initial value for A2A in the large coupling limit. The flow of and 
k is not affected by the initial running of A 2 (k) and quantities like Ap QR /A or m R /Ap 0R 
become independent of A 2 a if the coupling is sufficiently large. This qualitative discussion 
is confirmed by the numerical solution of the full set of equations presented in figs. 9 and 
10. For the fixed point value we obtain A2*(k) = 38.02. We point out that an analogous 
discussion for the large coupling region of Aia cannot be made. This can be seen by 
considering the mass term at the origin of the short distance potential (4.6) given by 
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u A (0, 0) = — k, a Xi A . Due to the pole of l^(u,rf) at uo = — 1 for n > 1/2 (cf. appendix C) 
one obtains the constraint 

« a Aia<1 . (4.61) 

In the limit Aia — > oo the mass term 2/t A Ai A at the minimum k of the potential at the 
critical temperature therefore remains finite. 

4.7 Universal equation of state for weak first order phase tran- 
sitions 

We presented in section 4.3 some characteristic quantities for the effective average potential 
which become universal at the phase transition for a sufficiently small quartic coupling 
^2A = A 2 a/A of the short distance potential Ua (4.6). The aim of this section is to 
generalize this observation and to find a universal scaling form of the equation of state for 
weak first order phase transitions. The equation of state relates the derivative of the free 
energy U = lim^o £4 to an external source, dU/dp = j. Here the derivative has to be 
evaluated in the outer convex region of the potential. For instance, for the meson model of 
strong interactions the source j is proportional to the average quark mass [31, 6] and the 
equation of state permits to study the quark mass dependence of properties of the chiral 
phase transition. We will compute the equation of state for a nonzero coarse graining 
scale k. It therefore contains information for quantities like the "classical" bubble surface 
tension in the context of Langer's theory of bubble formation which will be discussed in 
section 4.8. 

In three dimensions the U(2) x U(2) symmetric model exhibits a second order phase 
transition in the limit of a vanishing quartic coupling X 2 \ due to an enhanced 0(8) symme- 
try. In this case there is no scale present in the theory at the critical temperature. In the 
vicinity of the critical temperature (small \Sk\\ ~ \T C — T\) and for small enough A 2 a one 
therefore expects a scaling behavior of the effective average potential Uk and accordingly a 
universal scaling form of the equation of state. At the second order phase transition in the 
0(8) symmetric model there are only two independent scales that can be related to the 
deviation from the critical temperature and to the external source or tp. As a consequence 
the properly rescaled potential U / p\ or U / p( 5+1 )/ 2 (with the usual critical exponent 5) can 
only depend on one dimensionless ratio. A possible set of variables to represent the two 
independent scales are the renormalized minimum of the potential </?oii = (por/2) 1 ^ 2 (or the 
renormalized mass for the symmetric phase) and the renormalized field pn = (p^/2) 1 / 2 . 
The rescaled potential will then only depend on the scaling variable z = Pr/por [62]. An- 
other possible choice is the Widom scaling variable x = —5k a / p 1 ^ [88]. In the U(2)xU (2) 
symmetric theory A 2 a is an additional relevant parameter which renders the phase transi- 
tion first order and introduces a new scale, e.g. the nonvanishing value for the jump in the 
renormalized order parameter A<p R = (^-Por/2) 1 ^ 2 at the critical temperature or 5k a = 0. 
In the universal range we therefore observe three independent scales and the scaling form 
of the equation of state will depend on two dimensionless ratios. The rescaled potential 
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U I ifQ R can then be written as a universal function G 




(4.63) 



(4.62) 



The relation (4.62) is the scaling form of the equation of state we are looking for. At a 
second order phase transition the variable v vanishes identically and G(z, 0) describes the 
scaling equation of state for the model with 0(8) symmetry [62]. The variable v accounts 
for the additional scale present at the first order phase transition. We note that z = 1 
corresponds to a vanishing source and G(l,v) describes the temperature dependence of 
the free energy for j — 0. In this case v — 1 denotes the critical temperature T c whereas 
for T < T c one has v < 1. Accordingly v > 1 obtains for T > T c and (pon describes here 
the local minimum corresponding to the metastable ordered phase. The function G(z, 1) 
accounts for the dependence of the free energy on j for T = T c . 

We consider the scaling form (4.62) of the equation of state for a nonzero coarse graining 

1/2 

scale k. The renormalized field is given by ^r = Z k ' p. We pointed out in section 4.3 that 
there is a characteristic scale k 2 for the first order phase transition where the second local 
minimum of the effective average potential appears. For weak first order phase transitions 
one finds poR ~ &2- To observe the scaling form of the equation of state the infrared cutoff 
k has to run below k 2 with k <C k 2 . For the scale kf defined in eq. (4.18) we observe 
universal behavior to high accuracy (cf. fig. 10 and the corresponding universal ratios in 
table 4 for small A 2 a/Ai A ). The result for the universal function U kf /(fQ R = G kf (z,v) is 
presented in fig. 13. For v — 1 one has </?oi?(^/) = A<pon(kf) which denotes the critical 
temperature. Accordingly v > 1 denotes temperatures above and v < 1 temperatures 
below the critical temperature. One observes that Gk f (z, 1) shows two almost degenerate 
minima. (They become exactly degenerate in the limit k — > 0). For the given examples 
v = 1.18, 1.07 the minimum at the origin becomes the absolute minimum and the system 
is in the symmetric phase. In contrast, for v = 0.90, 0.74 the absolute minimum is located 
at z — 1 which characterizes the spontaneously broken phase. For small enough v the local 
minimum at the origin vanishes. 

We explicitly verified that the universal function G kf depends only on the scaling vari- 
ables z and v by choosing various values for 5k^ and for the quartic couplings of the short 
distance potential, Aia and A 2 a- In section 4.3 we observed that the model shows univer- 
sal behavior for a certain range of the parameter space. For given Aia and small enough 
A 2 a one always observes universal behavior. For Aia = 0.1, 2 and 4 it is demonstrated 
that (approximate) universality holds for A 2 a/Aia^1/2 (cf. fig. 10). For Aia around 2 one 
observes from figs. 9, 10 and table 4 that the system is to a good accuracy described by 
its universal properties for even larger values of A 2 a- The corresponding phase transitions 
cannot be considered as particularly weak first order. The universal function G kf therefore 
accounts for a quite large range of the parameter space. 
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Figure 13: Universal shape of the coarse grained potential (k = kf) as a function of the scaling variable 
z = pr/<Por = (pr/por) 1/2 for different values of v = A^ok/^ok = (Ap ii/ 'por) 1/2 ■ The employed values 
for v are (1) v = 1.18, (2) v = 1.07, (3) v = 1, (4) v = 0.90, (5) v — 0.74. For vanishing sources one has 
z = 1. In this case v = 1 denotes the critical temperature T c . Similarly v > 1 corresponds to T > T c with 
</?ofl denoting the minimum in the metastable ordered phase. 

We emphasize that the universal form of the effective potential given in fig. 13 depends 
on the scale kf where the integration of the flow equations is stopped (cf. eq. (4.18)). 
A different prescription for kf will, in general, lead to a different form of the effective 
potential. We may interpret this as a scheme dependence describing the effect of different 
coarse graining procedures. This is fundamentally different from nonuniversal corrections 
since Gk f is completely independent of details of the short distance or classical action 
and in this sense universal. A more quantitative discussion of this scheme dependence 
will be presented in the next section. We note that fluctuations on scales k < kf do not 
influence substantially the location of the minima of the coarse grained potential and the 
form of Uk(<fn) for ip R > (p 0R remains almost /c-independent. Here dU^Jdip = j(kf) with 
j(k f ) ~ ]im k ^ j(k) = j. 26 

We have established the scaling equation of state using the renormalized variables 
(pon and A<£>or- Alternatively the scaling equation of state can be presented in terms of the 
short distance parameters 5k,\, A2A and the unrenormalized field </?. For the 0(8) symmetric 

26 The role of massless Goldstone boson fluctuations for the universal form of the effective average 
potential in the limit k — > has been discussed previously for the 0(8) symmetric model [62]. 
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model (A2A = 0) this is known as the Widom scaling form of the equation of state [88] 
and in this case the relation between p R , po R and ip, <5ka is solely determined by critical 
exponents and amplitudes [62]. For the U{2) x U(2) symmetric model the dependence 
of <p R , <p 0R and Ap 0R on the parameters 5k\, A 2 a and the unrenormalized field ip can be 
expressed in terms of scaling functions. In general, only for limiting cases as 5k a = or 
A 2 a = the relation is determined by critical exponents and amplitudes. We will consider 
here these limits. 

We consider the renormalized minimum p 0R in two limits which are denoted by A<p 0R = 
<A).r(^ k a — 0) and <Pq R = Por(^2A = 0). The behavior of Aip 0R is described in terms of the 
exponent 9 according to eq. (4.14), 

Ap 0R ~ (A 2 a) 9/2 , 9 = 1.93. (4.64) 

The dependence of the minimum p^ R of the 0(8) symmetric potential on the temperature 
is characterized by the critical exponent u, 

~ (<S«a) v/2 , v = 0.882. (4.65) 

The exponent v for the 0(8) symmetric model is determined analogously to 9 as described 
in section 9 27 . We can also introduce a critical exponent ( for the jump of the unrenor- 
malized order parameter 

Ap ~ (A 2A ) f , C = 0.988 . (4.66) 

With 

pi ~ {5n K f, (3 = 0.451 (4.67) 
it is related to 9 and v by the additional index relation 

^ = ^ = 1.95. (4.68) 

We have verified this numerically. For the case 5k\ = \ 2 \ = one obtains 

j ~ <p s . (4.69) 

The exponent 5 is related to the anomalous dimension r\ via the usual index relation 
5 = (5 — 77) /(l + rj). From the scaling solution of eq. (4.55) we obtain 77 = 0.0224. 
With the help of the above relations one immediately verifies that for A 2 a = 

z~(-x)- /3 , v = (4.70) 

and for 5k\ = 

z~y~ c , v = l. (4.71) 
Here we have used the Widom scaling variable x and the new scaling variable y given by 



27 



For the 0(8) symmetric model (A2A = 0) we consider the minimum ip® R at k — 0. 
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4.8 Coarse graining and Langer's theory of bubble nucleation 

The coarse grained effective potential Uk results from the integration of fluctuations with 
momenta larger than ~ k. It is a nonconvex function whereas the standard effective poten- 
tial U = lim^o Uk has to be convex by its definition as a Legendre transform (cf. section 
2.2). The difference between U and Uk is due to fluctuations with characteristic length 
scales larger than the inverse coarse graining scale ~ k~ l . The role of these fluctuations 
for the approach to convexity has been established explicitly [106, 68, 69]. The study of 
first order phase transitions usually relies on the nonconvex "part" of the potential. As 
an illustration we consider the change of the system from the high temperature to the 
low temperature phase by bubble nucleation as described by Langer's theory [71]. On the 
one hand, the approach relies on the definition of a suitable coarse grained free energy Y k 
with a coarse graining scale k and, on the other hand, a saddle point approximation for 
the treatment of fluctuations around the "critical bubble" is employed. The problem is 
therefore separated into two parts: One part concerns the treatment of fluctuations with 
momenta q 2 ^ k 2 which are included in the coarse grained free energy. The second part 
deals with an estimate of fluctuations around the bubble for which only fluctuations with 
momenta smaller than k must be considered. These issues will be discussed in this section 
in a quantitative way. In particular, we will give a criterion for the validity of Langer's 
formula. 

One may consider a system that starts at some high temperature T > T c and investigate 
what happens as T is lowered as a function of time as for example during the evolution of 
the early universe [10]. For large enough temperature the origin of the potential (</? = 0) 
is the only minimum and the system is therefore originally in the symmetric phase. As 
T approaches T c a second local minimum develops at tpo > 0. This becomes the absolute 
minimum below T c . Nevertheless, the potential barrier prevents the system to change 
smoothly to the ordered phase. For a short while where T is in the vicinity of T c but below 
T c the system remains therefore in a state with higher energy density as compared to the 
state corresponding to the absolute minimum away from the origin. This is the so-called 
"false vacuum" in high energy physics or the metastable state in statistical physics. Such 
a state is unstable with respect to fluctuations which penetrate or cross the barrier. The 
picture is familiar from the condensation of vapor. The false vacuum corresponds to the 
supercooled vapor phase and the true vacuum to the fluid phase. Bubbles of the true 
vacuum (droplets) occur through thermal or quantum fluctuations 28 . If a bubble is large 
enough so that the decrease in volume energy exceeds the surface energy it will grow. The 
phase transition is completed once the whole space is filled with the true vacuum. On 
the other hand, small bubbles shrink due to the surface tension. The critical bubble is 
just large enough that it does not shrink. To be explicit we consider a spherical bubble 
where the bubble wall with "thickness" A is thin as compared to the bubble radius R, i.e. 
A <C R. In leading order the coarse grained free energy for such a bubble configuration 

28 In the real world the condensation of vapor is triggered by impurities but this is not the issue here. 
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can be decomposed in a volume and a surface term [72, 73] 

rf = -^R^ + 4rrR 2 a k . (4.73) 
In the thin wall approximation one obtains for the surface tension a k in our conventions 

a k = 2 J d^2Z k U k {v). (4.74) 
o 

For the difference in the free energy density e one has 

e = [/(0) -U(<p ) = -]imU (k). (4.75) 

We include in e fluctuations with arbitrarily small momenta. In contrast, the long wave- 
length contributions to the true surface tension 29 are effectively cut off by the characteristic 
length scale of the bubble surface. We include for the computation of the "classical" surface 
tension o k only fluctuations with q 2 ^k 2 . The modes with q 2 ^.k 2 contribute to the "fluc- 
tuation determinant" A k (cf. eq. (4.76)). The determination of a suitable coarse graining 
scale k for the computation of a k is discussed below. 

The critical bubble maximizes T k ^ with respect to the radius. It minimizes the coarse 
grained free energy with respect to other deformations because the spherical form is en- 
ergetically favorable. The critical bubble therefore represents a saddle point in the space 
of possible "bubble configurations" . In Langer's theory of bubble formation one considers 
a saddle point expansion around the critical bubble. There is exactly one negative mode 
that corresponds to the shrinking or growth of the bubble and there are infinitely many 
positive modes (there are also translational zero modes) 30 . The bubble nucleation rate f , 
which describes the probability per unit volume per unit time for the transition to the new 
vacuum, can be written in the form [72, 73, 71] 

f = A k exp(-rf [<p®]) = A k exp {-'-ffj (4.76) 

where rj^ is evaluated for ip^ corresponding to the critical bubble and approximated by 
(4.73). The exponential term with the coarse grained free energy rj^ denotes the lowest 
order or classical contribution. The prefactor A k contains several factors that depend on the 
details of the system under investigation. In particular, A k accounts for the contribution to 
the free energy from the fluctuations with momenta smaller than k. It depends on k through 
the effective ultraviolet cutoff for these fluctuations which is present since fluctuations with 

29 The true surface tension is a "measurable quantity" . It is independent of k and all fluctuations must 
be included. It therefore differs, in general, from a k which includes only part of the fluctuations. 

30 Langer's theory is not restricted to the thin wall approximation which is considered here for simplicity. 
In particular, the property of the critical bubble to represent a saddle point with exactly one negative 
mode is independent from the thin wall approximation. 
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momenta larger than k are already included in rjj. \ip\ ] 31 . Langer's formula for bubble 
nucleation amounts essentially to a perturbative one loop estimate of A k . 

For a determination of a useful choice of k it is convenient to place the discussion in 
a more general context which does not rely on the thin wall approximation or a saddle 
point approximation. What one is finally interested in is the free energy T[{p b ] for bubble 
configurations of a given shape. The "true critical bubble" ip b corresponds to a saddle 
point in the space of "bubble configurations" which are characterized by boundary condi- 
tions connecting the false and the true vacuum. The nucleation rate is then proportional 
exp— r[ipl\. The coarse graining can be seen as a convenient strategy to evaluate r[ip b ] 
by separating contributions from different momentum scales. We propose to choose the 
coarse graining scale k somewhat above but in the vicinity of the inverse thickness A -1 
of the bubble wall. We will argue below that in this case the corrections to the effective 
surface tension from the prefactor A k should be best accessible. 

In fact, we can write f = B exp — r where B contains dynamical factors and T[ipl] 
does not include contributions from fluctuations of the negative mode and the translational 
modes present for the critical bubble. The prefactor in eq. (4.76) can then be written as 

A k = Bexp~(Sk + Vk) (4.77) 

where 

S k = T[(p c b ] - T k [p c b ] , rj k = r k [<pt] - — -f . (4.78) 

The term r\ k includes the difference between the true critical bubble and the configuration 
used to estimate a k as well as a correction term to e to be discussed below. We first 
concentrate on 5 k which describes the difference between the free energy and the coarse 
grained free energy for the critical bubble. As mentioned above this is due to fluctuations 
with momenta q 2 ^ k 2 and incorporates the dominant /c-dependence of A k . Since the bubble 
provides for inherent effective infrared cutoff scales ~ R^ 1 or A -1 the contribution 5 k is 
both infrared and ultraviolet finite. The larger k, the more fluctuations are included in 5 k 
and from this point of view one wants to take k as low as possible. On the other hand, 
k should not be taken smaller than A -1 if the approximation used for a computation of 
T k relies on almost constant field configurations rather than real bubbles, as is usually the 
case. Only for k sufficiently large compared to A -1 the difference between an evaluation 
of the potential and kinetic terms in T k for almost constant field configurations (e.g. by 
a derivative expansion) rather than for bubbles remains small. In this way the technique 
of course graining combines a relatively simple treatment of the modes with q 2 k 2 for 
which the detailed properties of the bubble are irrelevant with an estimate of fluctuations 
around the bubble for which the short distance physics (q 2 ^ k 2 ) needs not to be considered 
anymore. It is clear that k is only a technical construct and for physical quantities the 
/c-dependence of 8 k and rj^ must cancel. More precisely, this concerns the sum 5 k + rj k + 
167rcr|/3e 2 . For thin wall bubbles the most important contribution to i] k is easily identified: 

31 The effective average action [37] afso provides the formai tool how the ultraviolet cutoff <~ k is imple- 
mented in the remaining functional integral for large length scale fluctuations. 
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By our definition of e we have included contributions from fluctuations with length scales 
kR. They should not be present in the effective action for a bubble with finite radius. 
Therefore rj k contains a correction term (167r/3)<rf (e~ 2 (i?) — e~ 2 ) which replaces effectively 
e by e(R) in eq. (4.76). We can evaluate e(R) in terms of the coarse grained free energy at 
a scale k R 

e(R) ~ U kR (0) - U kR ( Vo ) = -U (k R ) , k R = i . (4.79) 

For A <C R one should not confound k R with the coarse graining scale k since one has the 
inequality 

k R < £ fc . (4.80) 

Only for A ~ i? the clear separation between k R and disappears. We note that at the 
critical temperature one has R — > oo and therefore e(-R) = e . 

Since <7fc enters the nucleation rate (4.76) exponentially even small changes with k will 
have large effects. If one finds a strong dependence of a k on the coarse graining scale k 
this is only compatible with a large contribution from the higher orders of the saddle point 
expansion. The /c-dependence of a k therefore gives direct information about the validity 
of Langer's formula. We find a strong scale dependence of a k if the phase transition 
is characterized by large effective dimensionless couplings. A weak scale dependence is 
observed for small effective couplings. This gives a very consistent picture: The validity 
of the saddle point approximation typically requires small dimensionless couplings. In 
this case also the details of the coarse graining are not of crucial importance within an 
appropriate range of k. The remaining part of this section provides a detailed quantitative 
discussion. 

We consider the dependence of the effective average potential U k and the surface tension 
o k on the coarse graining scale k near the critical temperature T c for three examples in 
detail. They are distinguished by different choices for the quartic couplings Aia and A 2 a of 
the short distance potential U\ given by eq. (4.6). The choice Aia/A = O.f and A 2 a/A = 2 
corresponds to a strong first order phase transition with renormalized masses not much 
smaller than the cutoff scale A. The renormalized couplings will turn out small enough 
such that the notion of a coarse grained potential U k and a surface tension a k can be 
used without detailed information on the coarse graining scale within a certain range of 
k. In contrast we give two examples where the dependence of U k and a k on the coarse 
graining scale becomes of crucial importance. The choice Aia/A = 2 and A 2 a/A = 0.1 
leads to a weak first order phase transition with small renormalized masses and the system 
shows universal behavior (cf. sect. 4.3). For Aia/A = 4 and A 2 a/A = 70 one observes a 
relatively strong first order phase transition. Nevertheless for both examples the coarse 
grained potential and the surface tension show a similarly high sensitivity on the scale k. 

In figs. 14, 15 the scale dependence of U k is shown for a fixed temperature in the vicinity 
of T = T c or 5k a = 0. We plot U k in units of A 3 as a function of (f R /A 1/2 = (p R /2A) 1 ^ 2 = 
(Z k p/2h) 1 / 2 for different values of k. Each curve differs in k = Ae* by At = 1/18 and the 
first curve to the left with a negative curvature at the origin corresponds to t ~ —0.40 
for fig. 14 (t ~ -9.3 for fig. 15). For > tk - 0.40 (0 > tk, - 9.3) there is only 
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Figure 14: The coarse grained potential in dependence on the coarse graining scale k for fixed (almost 
critical) temperature. The example corresponds to Aia/A = 0.1, A 2 a/A = 2. 

one minimum of the potential away from the origin which lies below the plotted range in 
fig. 14 (15). Lowering k results in a successive inclusion of fluctuations on larger length 
scales and one observes the appearance of a second local minimum at the origin of the 
potential. The barrier between both minima increases and reaches its maximum value. As 
k is lowered further the potential barrier starts to decrease whereas the location of the 
minima becomes almost ^-independent and degenerate in height. In this region also the 
outer convex part of the potential shows an almost stable profile due to the decoupling 
of the massive modes. We stop the integration at the scale k — kf given by eq. (4.18). 
The coarse grained effective potential Uk with k = kf essentially includes all fluctuations 
with squared momenta larger than the scale Imf^l given by the curvature at the top of 
the potential barrier (cf. eq. (4.17)). Successive inclusion of fluctuations with momenta 
smaller than kf would result in a further decrease of the potential barrier. The flattening 
of the barrier is induced by the pole structure of the threshold functions l\{uS) appearing 
in the flow equations for the couplings (e.g. l\(u>) appearing in eq. (4.42) exhibits a pole at 
uo = — 1 according to ~ (uo + l) -1 / 2 for uo near —1 (cf. appendix C)). The argument 

uo corresponds to dimensionless mass terms as U' k (pR)/k 2 or {U k (pn) + 2p R U k (p R ))/k 2 . In 
the nonconvex region of the potential the curvature is negative. Since the pole cannot be 
crossed, negative U' k or U k must go to zero with k 2 and as a consequence the potential 
barrier flattens. For the scalar "vector" model the approach to convexity in the limit k — > 
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Figure 15: The coarse grained potential for Aia/A = 2, A 2 a/A = 0.1. 



has been studied analytically previously [106]. Here we are interested in the potential for 
a nonzero coarse graining scale k. 

The most significant difference between fig. 14 and fig. 15 is the /c-dependence of the 
potential barrier in a region where the minima become degenerate and almost independent 
of k. Fig. 14 shows a barrier with a weak scale dependence in the range of k where the 
location of the minima stabilizes. In contrast, in fig. 15 one observes a barrier with a strong 
scale dependence in this region. Figs. 11 and 12 (cf. sect. 4.6) show the relative height 
U (k) = U k ((p ) — U k (0) between the two local minima and the height of the potential 
barrier = £4(</?b) — £4(0) ((9C4/9</?) — 0, < ips < V^o) as a function of 

t = \n(k/A). Accordingly one observes that for Aia/A = 0.1, A 2 a/A = 2 the top of the 
potential barrier shows a weak scale dependence in a region of k where U (k) is small 
whereas for Aia/A = 2, A 2 a/A = 0.1 the top of the potential barrier depends strongly on 
the coarse graining scale in this region. The surface tension a k is displayed in fig. 16 and 
shows a corresponding behavior. Here we consider a k also for the short distance parameters 
Aia/A = 4, A 2 a/A = 70. In fig. 16 the surface tension o k is normalized to o" max ((x max /A 2 = 
1.67 x 10~ 2 (8.41 x 10- n )(1.01 x 10~ 3 )_ for A 1A /A = 0.1(2)(4), A 2A /A = 2(0.1)(70)) and 
given as a function of \n(k/kf). 32 For Aia/A = 0.1, A 2 a/A = 2 the curve exhibits a small 

32 The integration according to eq. (4.74) is performed between the two zeros ip = and ip = <p' £,tp of 

u k . 
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curvature around its maximum and a max ~ cr fe=fe/ . One observes for the second and the 
third example a comparably large curvature around a max and o k=kf becomes considerably 
smaller than the maximum. One may consider what happens if this scale dependence is 



<7 n 




In (k/kf) 

Figure 16: The normalized surface tension CTfe/(T max as a function of ln(k/kf). The short distance 
parameters are (1) Aia/A = 0.1, A 2A /A = 2, (2) Aia/A = 2, A 2A /A = 0.1, (3) Aia/A = 4, A 2A /A = 70. 

not taken into account correctly. If one takes the difference in a k from its maximum value 
to o"fe=fc / as a rough measure for its uncertainty this is about 30% for the last example. 
Entering exponentially in eq. (4.76) this would lead to tremendous errors. 

We have pointed out above that the coarse graining scale k should not be taken smaller 
than the inverse bubble wall thickness A -1 . This ensures that the detailed properties of 
the bubble are irrelevant for the computation of the "classical" surface tension a k . We 
estimate the bubble wall thickness A by 



A(fc) = 2Z k 



PoO) 



(4.81) 



where we have taken the gradient energy as half the total surface energy and approximated 
the mean field gradient at the bubble wall by tpo/A. For the given examples we observe 
A^ 1 (kf) ~ kf/4. We choose k ~ kf as the coarse graining scale. 

In order to quantify the differences between the three examples we have displayed some 
characteristic quantities in table 5. The renormalized couplings 



Air = U£ f (p 0R ) , \ 2 r = 4:d T U kf (p 0R ) 



(4.82) 
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A2A 


Air 


Mr 


m< R 


m c R 


m 2R 


k f 


"a" 


~Y 


VTl R 


m c R 


m c 9 R 


T 


A 


A 


0.1 


2 


0.228 


8.26 


0.235 


1.55 x 10- 1 


6.62 x 10- 1 


1.011 x 10- 1 


2 


0.1 


0.845 


15.0 


0.335 


2.04 x 10~ 5 


6.10 x 10" 5 


1.145 x 10~ 5 


4 


70 


0.980 


16.8 


0.341 


6.96 x 10~ 2 


2.04 x 10- 1 


3.781 x 10~ 2 



Table 5: The effective dimensionless couplings XiR/m c R and \2r/iti c r . The couplings and the mass terms 
m R , m\ R are evaluated at the scale kf and 5k\ = 0. 



are normalized to the mass term 

m c R = (2 Po rAir) 1/2 . (4.83) 

In addition we give the mass term 

mm = (PorA 2 r) 1/2 (4.84) 

corresponding to the curvature of the potential in the direction of the second invariant 
t. In comparison with figs. 14-16 one observes that the smaller the effective couplings 
the weaker the scale dependence of Uk and a^. In particular, a reasonably weak scale 
dependence of Uk and Ok requires 



Air _ 1 m c R 



< 1 . (4.85) 



This establishes a quantitative criterion for the range where Langer's theory can be used 
without paying too much attention to the precise definition of the coarse graining. Com- 
parison with fig. 10 (cf. sect. 9) shows that this condition is not realized for the range 
of couplings leading to universal behavior and for large Aia/A. The only region where 
the saddle point expansion is expected to converge reasonably well is for small Aia/A 
and large A2A/A. The second and the third example given in fig. 16, which exhibit a 
strong ^-dependence, show similar values for the effective couplings. More precisely, for 
the relatively strong phase transition with slightly larger effective couplings one observes 
an increased scale dependence as compared to the weak phase transition. Therefore, the 
strength of the transition is not the primary criterion for the convergence of the saddle 
point expansion. In table 5 also the renormalized masses in units of A which indicate the 
strength of the phase transition and kf/A are presented. 

In addition to the dependence on k, the coarse grained free energy depends also on 
the precise shape of the infrared cutoff function Rk{q) (cf. eq. 4.24). Analytical studies in 
the Abelian Higgs model indicate [107] that this scheme dependence is rather weak for the 
effective potential and the surface tension. 
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In summary, we have shown that the coarse grained free energy cannot be defined with- 
out detailed information on the coarse graining scale k unless the effective dimensionless 
couplings are small. Only for small couplings we observe a weak /c-dependence of the sur- 
face tension in a range where the location of the minima of the potential remains almost 
fixed. There is a close relation between the dependence of the coarse grained free energy on 
the coarse graining scale and the reliability of the saddle point approximation in Langer's 
theory of bubble nucleation. For a strong /c-dependence of U k a small variation in the coarse 
graining scale can induce large changes in the predicted nucleation rate in lowest order in 
a saddle point approximation. In this case the /c-dependence of the prefactor A k has also 
to be computed. For strong dimensionless couplings a realistic estimate of the nucleation 
rate therefore needs the capability to compute InAk with the same accuracy as 16nal/3e 2 
and a check of the cancelation of the /c-dependence in the combined expression (4.76). Our 
observation that the details of the coarse graining prescription become less important in 
the case of small dimensionless couplings is consistent with the fact that typically small 
couplings are needed for a reliable saddle point approximation for A k . For the electroweak 
high temperature phase transition a small ^-dependence of a k is found for a small mass 
Mh of the Higgs scalar whereas for Mh near the W-boson mass the picture resembles our 
fig. 15 [115]. This corresponds to the observation [116, 117, 118] that the saddle point 
approximation around the critical bubble converges well only for a small enough mass of 
the Higgs scalar. 

4.9 Conclusions 

We have presented in this chapter a detailed investigation of the phase transition in three 
dimensional models for complex 2x2 matrices. They are characterized by two quartic 
couplings Aia and A2A. In the limit Aia — > 00, A2A — > 00 this also covers the model of 
unitary matrices. The picture arising from this study is unambiguous: 

(1) One observes two symmetry breaking patterns for A2A > and A2A < respectively. 
The case A2A = denotes the boundary between the two phases with different symmetry 
breaking patterns. In this special case the theory exhibits an enhanced 0(8) symmetry. The 
phase transition is always first order for the investigated symmetry breaking U(2) xU(2) — > 
U(2) (A 2 a > 0). For A 2 a = the 0(8) symmetric Heisenberg model is recovered and one 
finds a second order phase transition. 

(2) The strength of the phase transition depends on the size of the classical quartic 
couplings Aia/A and A 2 a/A. They describe the short distance or classical action at a 
momentum scale A. The strength of the transition can be parametrized by m c R /A with 
rrf R a characteristic inverse correlation length at the critical temperature. For fixed A2A 
the strength of the transition decreases with increasing Aia- This is analogous to the 
Coleman- Weinberg effect in four dimensions. 

(3) For a wide range of classical couplings the critical behavior near the phase tran- 
sition is universal. This means that it becomes largely independent of the details of the 
classical action once everything is expressed in terms of the relevant renormalized parame- 
ters. In particular, characteristic ratios like m^/Ap _R (critical inverse correlation length in 
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the ordered phase over discontinuity in the order parameter) or m^/Apo/? (same for the 
disordered phase) are not influenced by the addition of new terms in the classical action 
as far as the symmetries are respected. 

(4) The range of short distance parameters Aia, A 2 a for which the phase transition 
exhibits universal behavior is not only determined by the strength of the phase transition 
as measured by m c R /A. For a given Aia/ A and small enough A 2 a/A one always observes 
universal behavior. In the range of small Aia/ A the essential criterion for universal behavior 
is given by the size of A 2 a/Aiaj with approximate universality for A 2 a < Aia- For strong 
couplings universality extends to larger A 2 a/Aia and occurs for much larger m c R /A (cf. 
table 4). 

(5) We have investigated how various characteristic quantities like the discontinuity in 
the order parameter Apo or the corresponding renormalized quantity Apoi? or critical corre- 
lation lengths depend on the classical parameters. In particular, at the critical temperature 
one finds universal critical exponents for not too large A 2 a, 

ApoR ~ (A 2A ) e , = 1.93, / x 

4.86 

Ap ~ (A 2A ) 2C , C = 0.988. 

These exponents are related by a scaling relation to the critical correlation length and 
order parameter exponents v and (3 of the 0(8) symmetric Heisenberg model according to 
9/C, — u//3 — 1.95 [y = 0.882, j3 = 0.451 in our calculation for A 2 a = 0). Small values 
of A 2 a can be associated with a perturbation of the 0(8) symmetric model and 6,( are 
related to the corresponding crossover exponents. On the other hand, Ap R (Ap ) becomes 
independent of A 2 a in the infinite coupling limit. 

(6) We have computed the universal equation of state. The equation of state relates 
the derivative of the free energy U to an external source, dU/dip = j. From there one 
can extract universal ratios e.g. for the jump in the order parameter (Apof>/m c R = 0.592) 
or for the ratios of critical correlation lengths in the disordered (symmetric) and ordered 
(spontaneously broken) phase {rriQ R /m R = 0.746). It specifies critical couplings (Xm/m R = 
0.845, A 2 ij/m^ = 15.0). The universal behavior of the potential for large field arguments 
p ^> po is U ~ p 3 /( 1+r 7) provided pn is sufficiently small as compared to A. Here the critical 
exponent rj which characterizes the dependence of the potential on the unrenormalized 
field p is found to be rj = 0.022. For large p the universal equation of state equals the one 
for the 0(8) symmetric Heisenberg model and rj specifies the anomalous dimension or the 
critical exponent S — (5 — 77) / (1 + rj). The equation of state is computed for a nonzero 
coarse graining scale k. It therefore contains information for quantities like the "classical" 
bubble surface tension in the context of Langer's theory of bubble formation. 

(7) We have investigated the dependence of the coarse grained effective potential Uk(p) 
and the "classical" surface tension on the coarse graining scale k with special emphasis 
on the question of the validity of Langer's theory of bubble formation. We find a strong 
scale dependence of [/& and aj, if the phase transition is characterized by large dimensionless 
couplings. A weak scale dependence is observed for small effective couplings. There is a 
close relation between the dependence of the coarse grained free energy on the coarse 
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graining scale and the reliability of the saddle point approximation in Langer's theory of 
bubble nucleation. A strong /c-dependence of is only compatible with a large contribution 
from the higher orders of the expansion. We obtain a very consistent picture: The validity 
of the saddle point approximation typically requires small dimensionless couplings. In 
this case also the details of the coarse graining are not of crucial importance within an 
appropriate range of k. The quantitative criterion for the validity of Langer's formula is 
in our case \i R /m c R <C 1. 

(8) Our method is not restricted to the study of the universal behavior. We can compute 
the effective potential for arbitrary values of the initial parameters and have done this for 
particular examples. 

The uncertainties of our results induced by the numerical integration of the flow equa- 
tions are well under control and small. They are negligible compared to the expected error 
induced by our truncation. For a significant improvement of our treatment one would have 
to include higher order terms in the derivative expansion employed for the effective average 
action. For weak first order or second order phase transitions we expect the error to be 
related to the anomalous dimension rj = 0.022. For the special case of the enhanced 0(8) 
symmetry one can compare e.g. with known values for critical exponents obtained by other 
methods [67, 89, 90]. A comparison of our results for the critical exponents (3 and v with 
the results of the most sophisticated calculations show agreement within a few per cent. 
The anomalous dimension is also well determined, even though it is most affected by our 
truncation. 

Finally, we should mention that our approach can be extended to systems with reduced 
SU (N) x SU (N) symmetry. They obtain by adding to the classical potential a term 
involving the invariant £ = detip + det</?L (Note that £ is not invariant with respect to 
U(N) x U(N)). This will give an even richer pattern of phase transitions and permits a 
close contact to realistic meson models in QCD where the axial anomaly is incorporated. 
Finally one can extend the three dimensional treatment to a four dimensional study of field 
theories at nonvanishing temperature. How this can be used to approach the chiral phase 
transition in QCD is presented in the next section. 
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5 The equation of state for two flavor QCD 
5.1 Introduction 

We have pointed out in section 1 that strong interactions in thermal equilibrium at high 
temperature T differ in important aspects from the well tested vacuum or zero tempera- 
ture properties. A phase transition at some critical temperature T c or a relatively sharp 
crossover may separate the high and low temperature physics [25]. Many experimental 
activities at heavy ion colliders [119] search for signs of such a transition. It was realized 
early that the transition should be closely related to a qualitative change in the chiral con- 
densate according to the general observation that spontaneous symmetry breaking tends to 
be absent in a high temperature situation. A series of stimulating contributions [31, 26, 27] 
pointed out that for sufficiently small up and down quark masses, m u and m d , and for a 
sufficiently large mass of the strange quark, m s , the chiral transition is expected to belong 
to the universality class of the 0(4) Heisenberg model. It was suggested [26, 27] that a 
large correlation length may be responsible for important fluctuations or lead to a disori- 
ented chiral condensate [28] . The question how small m u and ma would have to be in order 
to see a large correlation length near T c and if this scenario could be realized for realistic 
values of the current quark masses remained, however, unanswered. The reason was the 
missing link between the universal behavior near T c and zero current quark mass on one 
hand and the known physical properties at T = for realistic quark masses on the other 
hand. 

It is the purpose of the present section to provide this link [33]. We present here the 
equation of state for two flavor QCD within an effective quark meson model. The equation 
of state expresses the chiral condensate as a function of temperature and the average 

current quark mass rh = (m u + m d )j c l. This connects explicitly the universal critical 
behavior for T — * T c and rh — > with the temperature dependence for a realistic value 
m p h ys - Since our discussion covers the whole temperature range < T&1.7T C we can 
fix 

^phys such that the (zero temperature) pion mass is m n — 135 MeV. The condensate 



plays here the role of an order parameter. Its precise definition will be given in 



dependence of in the chiral limit rh — 0. Here the lower curve is the full result 

for arbitrary T whereas the upper curve corresponds to the universal scaling form of the 
equation of state for the 0(4) Heisenberg model. We see perfect agreement of both curves 
for T sufficiently close to T c = 100.7 MeV. This demonstrates the capability of our method 
to cover the critical behavior and, in particular, to reproduce the critical exponents of the 
0(4)-model. We have determined (cf. section 5.7) the universal critical equation of state 
as well as the non-universal amplitudes. This provides the full functional dependence of 
(ipip^J {T,rh) for small T — T c and rh. The curves (b), (c) and (d) are for non-vanishing 
values of the average current quark mass rh. Curve (c) corresponds to m p h ys or, equivalently, 
m^T = 0) = 135 MeV. One observes a crossover in the range T = (1.2 — 1.5)T C . The 
0(4) universal equation of state (upper curve) gives a reasonable approximation in this 




Curve (a) gives the temperature 
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Figure 17: The plot shows the chiral condensate (tpip) as a function of temperature T. Lines (a), (6), (c), 
(d) correspond at zero temperature to = 0, 45 MeV, 135 MeV, 230 MeV, respectively. For each pair of 
curves the lower one represents the full T-dependence of (ipi>) whereas the upper one shows for comparison 
the universal scaling form of the equation of state for the 0(4) Heisenberg model. The critical temperature 
for zero quark mass is T c = 100.7 MeV. The chiral condensate is normalized at a scale fc$ ~ 620 MeV. 

temperature range. The transition turns out to be much less dramatic than for rh — 0. 
We have also plotted in curve (b) the results for comparably small quark masses ~ 1 MeV, 
i.e. rh = m p h ys /10, for which the T = value of m n equals 45 MeV. The crossover is 
considerably sharper but a substantial deviation from the chiral limit remains even for 
such small values of rh. In order to facilitate comparison with lattice simulations which are 
typically performed for larger values of m n we also present results for m n (T = 0) = 230 MeV 
in curve (d). One may define a "pseudocritical temperature" T pc associated to the smooth 
crossover as the inflection point of (vpip^ (T) as usually done in lattice simulations. Our 

results for this definition of T pc are denoted by T^) and are presented in table 6 for the 
four different values of rh or, equivalently, m n (T = 0). The value for the pseudocritical 
temperature for m n = 230 MeV compares well with the lattice results for two flavor QCD 
(cf. section 5.7). One should mention, though, that a determination of T pc according to 
this definition is subject to sizeable numerical uncertainties for large pion masses as the 
curve in figure 17 is almost linear around the inflection point for quite a large temperature 
range. A problematic point in lattice simulations is the extrapolation to realistic values 
of m n or even to the chiral limit. Our results may serve here as an analytic guide. The 
overall picture shows the approximate validity of the 0(4) scaling behavior over a large 
temperature interval in the vicinity of and above T c once the (non-universal) amplitudes 
are properly computed. 

A second important result of our investigations is the temperature dependence of the 
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MeV 





45 


135 


230 


T (i) 

J-pc 

MeV 


100.7 


~ 110 


~ 130 


~ 150 


T (2) 
- 1 pc 

MeV 


100.7 


113 


128 





Table 6: The table shows the critical and "pseudocritical" temperatures for various values of the zero 
temperature pion mass. Here T$ is defined as the inflection point of (ipip) (T) whereas T$ is the location 
of the maximum of the sigma correlation length (see section 5.6). 



space-like pion correlation length m^ 1 (T). (We will often call m n (T) the temperature 
dependent pion mass since it coincides with the physical pion mass for T = 0.) The plot 
for m^T) in figure 18 again shows the second order phase transition in the chiral limit 
rh = 0. For T <T C the pions are massless Goldstone bosons whereas for T > T c they form 
with the sigma particle a degenerate vector of 0(4) with mass increasing as a function 
of temperature. For m = the behavior for small positive T — T c is characterized by 
the critical exponent u, i.e. m w (T) = (^y 1 T C ((T - T c ) /T c f and we obtain v = 0.787, 
£ + = 0.270. For rh > we find that m 7r (T) remains almost constant for T ^ T c with only 
a very slight dip for T near T c /2. For T > T c the correlation length decreases rapidly and 
for T ^> T c the precise value of m becomes irrelevant. We see that the universal critical 
behavior near T c is quite smoothly connected to T = 0. The full functional dependence 
of m^T, m) allows us to compute the overall size of the pion correlation length near the 
critical temperature and we find m^TpJ ~ 1. 7771^(0) for the realistic value m p hy S - This 
correlation length is even smaller than the vacuum (T = 0) one and gives no indication for 
strong fluctuations of pions with long wavelength. It would be interesting to see if a decrease 
of the pion correlation length at and above T c is experimentally observable. It should be 
emphasized, however, that a tricritical behavior with a massless excitation remains possible 
for three flavors. This would not be characterized by the universal behavior of the 0(4)- 
model. We also point out that the present investigation for the two flavor case does not 
take into account a speculative "effective restoration" of the axial £7a(1) symmetry at high 
temperature [31, 32]. We will comment on these issues in section 5.8. 

Our method is based on the effective average action T k [37] which has been introduced 
in section 2. Since in most cases the flow equation (2.10) for T k can not be solved exactly 
the capacity to devise useful truncations in a non-perturbative context becomes crucial. 
This requires first of all an identification of the degrees of freedom which are most relevant 
for a given problem. In the present work we concentrate on the chiral aspects of QCD 33 . 
Spontaneous chiral symmetry breaking occurs through the expectation value of a (complex) 
scalar field <3> afe which transforms as (N, N) under the chiral flavor group SUl(N) x SUr(N) 

33 For a study of chiral symmetry breaking in QED using related exact renormalization group techniques 
see ref. [121]. 



5.1 Introduction 



75 




100 - 



I 1 1 1 1 L 1 1 1 1 

20 40 60 80 100 120 140 160 180 

Tj MeV 

Figure 18: The plot shows as a function of temperature T for three different values of the average 
light current quark mass rh. The solid line corresponds to the realistic value m = m p hys whereas the dotted 
line represents the situation without explicit chiral symmetry breaking, i.e., rh — 0. The intermediate, 
dashed line assumes rh = m phys /10. 

with N the number of light quark flavors. More precisely, the expectation value 

($ ab } = a 5 ab (5.1) 

induces for a ^ a spontaneous breaking of the chiral group to a vector-like subgroup, 
SU L (N) x SU R (N) — ► SU L+R (N) = SU V (N). In addition, non-vanishing current quark 
masses m u ,mrf,m s break the chiral group explicitly and also lift the SUy(N) degeneracy 
of the spectrum if they are unequal. The physical degrees of freedom contained in the 
field <E> a & are pseudoscalar and scalar mesons which can be understood as quark-antiquark 
bound states. It is obvious that any analytical description of the chiral transition has to 
include at least part of these (pseudo-) scalar fields as the most relevant degrees of freedom. 

In the present work we use for k smaller than a "compositeness scale" ~ 600 MeV a 
description in terms of <£> a b and quark degrees of freedom. The quarks acquire a constituent 
quark mass M q through the chiral condensate <r which forms in our picture for k x ss — 
400 MeV. This effective quark meson model can be obtained from QCD by "integrating 
out" the gluon degrees of freedom and introducing fields for composite operators [44, 45]. 
This will be explained in more detail in the first part of section 5.2. In this picture the 
scale k$ is associated to the scale at which the formation of mesonic bound states can be 
observed in the flow of the effective (momentum dependent) four-quark interaction. We 
will restrict our discussion here to two flavor QCD with equal quark masses m u = rrid = rh. 
Since in this case the scalar triplet ao and the pseudoscalar singlet (associated with the 
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rf) have typical masses around 34 1 GeV we will neglect them for k < This reduces 
the scalar degrees of freedom of our effective model to a four component vector of 0(4), 
consisting of the three pions and the "sigma resonance" . 

We imagine that all other degrees of freedom besides the quarks ip and the scalars <3> are 
integrated out. This is reflected in the precise form of the effective average action T k<j? [tp, $] 
at the scale k$ which serves as an initial value for the solution of the flow equation. The flow 
of I\ [tp, <3>] for k < k<$, is then entirely due to the quark and meson fluctuations which are 
not yet included in T k<s> [ifj, $]. Obviously, the initial value T k<s> may be a quite complicated 
functional of ip and $ containing, in particular, important non-local behavior. We will 
nevertheless use a rather simple truncation in terms of standard kinetic terms and a most 
general form of the scalar potential U k) i.e. 35 



Here T k is invariant under the chiral flavor symmetry SUl(2) x SUr(2) and the only explicit 
symmetry breaking arises through the source term j ~ in. We will consider the flow of 
the most general form of U k consistent with the symmetries (without any restriction to 
a polynomial form as typically used in a perturbative context). On the other hand, our 
approximations for the kinetic terms are rather crude and parameterized by only two 
running wave function renormalization constants, Zq> >k and Z^. The same holds for the 
effective Yukawa coupling h k . The main approximations in this work concern 

(i) the simple form of the derivative terms and the Yukawa coupling, in particular, the 

neglect of higher derivative terms (and terms with two derivatives and higher powers 
of $). This is partly motivated by the observation that at the scale k$ and for small 
temperatures the possible strong non-localities related to confinement affect most 
likely only the quarks in a momentum range q 2 ^ (300 MeV) 2 . Details of the quark 
propagator and interactions in this momentum range are not very important in our 
context (see section 5.2). 

(ii) the neglect of interactions involving more than two quark fields. This is motivated 

by the fact that the dominant multi-quark interactions are already incorporated in 
the mesonic description. Six-quark interactions beyond those contained effectively 

34 More precisely, because of the anomalous Ua(1) breaking in QCD these mesons are significantly heavier 
than the remaining degrees of freedom in the range of scales k where the dynamics of the model is strongly 
influenced by mesonic fluctuations. The situation becomes more involved if the model is considered at 
high temperature which is discussed in section 5.8. 

35 0ur Euclidean conventions (h k is real) are specified in refs. [46, 77]. 
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in Uk could be related to baryons and play probably only a minor role for the meson 
physics considered here. 

We will choose a normalization of ip, $ such that Z^^ = = 1. We therefore need 
as initial values at the scale kq> the scalar wave function renormalization Z<p^ and the 
shape of the potential Uk 9 - We will make here the important assumption that Z$^ is 
small at the compositeness scale (similarly to what is usually assumed in Nambu-Jona- 
Lasinio-like models). This results in a large value of the renormalized Yukawa coupling 
hk = Z^ 2 Z^ k hk- A large value of is phenomenologically suggested by the comparably 
large value of the constituent quark mass M q . The latter is related to the value of the 
Yukawa coupling for k — > and the pion decay constant f w = 92.4 MeV by M q = hf w /2 
(with h = hk=o), and M q ~ 300 MeV implies h 2 /An ~ 3.4. For increasing k the value of 
the Yukawa coupling grows rapidly for k ^ M q . Our assumption of a large initial value 
for hkf, is therefore equivalent to the assumption that the truncation (5.2) can be used 
up to the vicinity of the Landau pole of hj,. The existence of a strong Yukawa coupling 
enhances the predictive power of our approach considerably. It implies a fast approach of 
the running couplings to partial infrared fixed points [46]. In consequence, the detailed form 
of Uk<s, becomes unimportant, except for the value of one relevant parameter corresponding 
to the scalar mass term m\ . In this work we fix m\ such that /„. = 92.4 MeV for 
m n = 135 MeV. The possibility of such a choice is highly non-trivial since f w can actually 
be predicted [46] in our setting within a relatively narrow range. The value/,,- = 92.4 MeV 
(for m n = 135 MeV) sets our unit of mass for two flavor QCD which is, of course, not 
directly accessible by observation. Besides m| (or /„.) the other input parameter used in 
this work is the constituent quark mass M q which determines the scale fc$ at which 
becomes very large. We consider a range 300 MeV ^ M q ^ 350 MeV and find a rather 
weak dependence of our results on the precise value of M q . We also observe that the limit 
hk^ — ► oo can be considered as the lowest order of a systematic expansion in which is 
obviously highly non-perturbative. 

A generalization of our method to the realistic case of three light flavors is possible 
and work in this direction is in progress. For the time being we expect that many features 
found for N = 2 will carry over to the realistic case, especially the critical behavior for 
T — > T c and m — > (for fixed m s ^ 0). Nevertheless, some quantities like (vpip^ (T = 0), 
the difference between f n for realistic quark masses and m = or the mass of the sigma 
resonance at T = may be modified. This will also affect the non-universal amplitudes 
in the critical equation of state and, in particular, the value of T c . In the picture of 
the two flavor quark meson model these changes occur through an effective temperature 
dependence of the initial values of couplings at the scale k$. This effect, which is due 
to the temperature dependence of effects from fluctuations not considered in the present 
work is discussed briefly in section 5.8. It remains perfectly conceivable that this additional 
temperature dependence may result in a first order phase transition or a tricritical behavior 
for realistic values of m for the three flavor case. Details will depend on the strange quark 
mass. We observe, however, that the temperature dependence in the limit rh — > involves 
for T < T c only information from the running of couplings in the range k ^ 300 MeV. 
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(The running for k ^ 3T effectively drops out in the comparison between the thermal 
equilibrium results and those for T = 0.) In this range of temperatures our model should 
be quite reliable. 

Finally, we mention that we have concentrated here only on the ^-dependent part of 
the effective action which is related to chiral symmetry breaking. The ^-independent part 
of the free energy also depends on T and only part of this temperature dependence is 
induced by the scalar and quark fluctuations considered in the present work. Most likely, 
the gluon degrees of freedom cannot be neglected for this purpose. This is the reason why 
we do not give results for "overall quantities" like energy density or pressure as a function 
of T. 

The text is organized as follows: In section 5.2 we review the linear quark meson model 
at vanishing temperature. We begin with an overview of the different scales appearing in 
strong interaction physics. Subsequently, the flow equations for the linear quark meson 
model are introduced and their approximate partial fixed point behavior is discussed in 
detail leading to a "prediction" of the chiral condensate (^pip^j. In section 5.5 the exact 
renormalization group formulation of field theories in thermal equilibrium is given. It is 
demonstrated how mass threshold functions in the flow equations smoothly decouple all 
massive Matsubara modes as the temperature increases, therefore leading to a "dimensional 
reduction" of the model. Section 5.6 contains our results for the linear quark meson model 
at non-vanishing temperature. Here we discuss the T-dependences of the parameters and 
physical observables of the linear quark meson model in detail for a temperature range 
< T ^ 170 MeV including the (pseudo) critical temperature T c of the chiral transition. 
The critical behavior of the model near T c and rh = 0, where rh denotes the light average 
current quark mass, is carefully analyzed in section 5.7. There we present the universal 
scaling form of the equation of state including a fit for the corresponding scaling function. 
Also the universal critical exponents and amplitude ratios are given there. The effects of 
additional degrees of freedom of strong interaction physics not included in the linear 0(4) 
symmetric quark meson model are addressed in section 5.8. Here we also comment on 
differences between the linear quark meson model and chiral perturbation theory. Some 
technical details concerning the quark mass term and the definition of threshold functions 
at vanishing and non-vanishing temperature are presented in three appendices. 

5.2 The quark meson model 

Before discussing the finite temperature behavior of strong interaction physics we will 
review some of its zero temperature features. This will be done within the framework of a 
linear quark meson model as an effective description for QCD for scales below the mesonic 
compositeness scale of approximately ~ 600 MeV. Relating this model to QCD in a 
semi-quantitative way in subsection 5.3 will allow us to gain some information on the 
initial value for the effective average action at the compositeness scale k$. We emphasize, 
however, that the quantitative aspects of the derivation of the effective quark meson model 
from QCD will not be relevant for our practical calculations in the mesonic sector. This 
is related to the "infrared stability" for large Yukawa coupling h k ^ as discussed in the 
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introduction and which will be made quantitative in subsection 5.4. 
5.3 A short (scale) history of QCD 

For an evaluation of the trace on the right hand side of the flow equation (2.10) only a 
small momentum range q 2 ~ k 2 contributes substantially. One therefore only needs to take 
into account those fluctuations which are important in this momentum interval. Here we 
are interested in the description of chiral symmetry breaking. The relevant fluctuations in 
relation to this phenomenon may change with the scale k and we begin by summarizing 
the qualitatively different scale intervals which appear for meson physics in QCD. Some of 
this will be explained in more detail in the remainder of this section whereas other aspects 
are well known. Details of this discussion may also be found in refs. [44, 46, 122]. We will 
distinguish five qualitatively different ranges of scales: 

1. At sufficiently high momentum scales, say, 

k>k p ~ 1.5 GeV 

the relevant degrees of freedom of strong interactions are quarks and gluons and their 
dynamics is well described by perturbative QCD. 

2. For decreasing momentum scales in the range 

fc$ ~ 600 MeV £ k £ k p ~ 1.5 GeV 

the dynamical degrees of freedom are still quarks and gluons. Yet, as k is lowered part 
of their dynamics becomes dominated by effective non-local four quark interactions 
which cannot be fully accessed perturbatively. 

3. At still lower scales this situation changes dramatically. Quarks and gluons are 
supplemented by mesonic bound states as additional degrees of freedom which are 
formed at a scale k$ ~ 600 MeV. We emphasize that k^ is well separated from 
Aqcd — 200 MeV where confinement sets in and from the constituent masses of the 
quarks M q ~ (300 — 350) MeV. This implies that below the compositeness scale kq> 
there exists a hybrid description in term of quarks and mesons! It is important to note 
that for scales not too much smaller than kq> chiral symmetry remains unbroken. This 
situation holds down to a scale k x sB — 400 MeV at which the scalar meson potential 
develops a non-trivial minimum thus breaking chiral symmetry spontaneously. The 
meson dynamics within the range 

k xSB ^ 400 MeV £ k £ k$ ~ 600 MeV 

is dominated by light current quarks with a strong Yukawa coupling h 2 k / (47r) a s (k) 
to mesons. We will thus assume that the leading gluon effects are included below k$ 
already in the formation of mesons. Near k x sB also fluctuations of the light scalar 
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mesons become important as their initially large renormalized mass approaches zero. 
Other hadronic bound states like vector mesons or baryons should have masses larger 
than those of the lightest scalar mesons, in particular near k x sB, and give therefore 
only subleading contributions to the dynamics. This leads us to a simple linear model 
of quarks and scalar mesons as an effective description of QCD for scales below k$. 

4. As one evolves to scales below k x sB the Yukawa coupling decreases whereas a s in- 
creases. Of course, getting closer to Aq CD it is no longer justified to neglect in the 
quark sector the QCD effects which go beyond the dynamics of the effective quark 
meson model in our truncation (5.2). On the other hand, the final IR value of 
the Yukawa coupling h is fixed by the typical values of constituent quark masses 
M q ~ 300 MeV to be h 2 /(Air) ~ 3.4. One may therefore speculate that the domina- 
tion of the Yukawa interaction persists even for the interval 

M q ~ 300 MeV < k < k xSB ~ 400 MeV 

below which the quarks decouple from the evolution of the mesonic degrees of freedom 
altogether. Of course, details of the gluonic interactions are expected to be crucial 
for an understanding of quark and gluon confinement. Strong interaction effects may 
dramatically change the momentum dependence of the quark propagator for k and 
q 2 around Aq C d- Yet, there is no coupling of the gluons to the color neutral mesons. 
As long as one is only interested in the dynamics of the mesons one is led to expect 
that confinement effects are quantitatively not too important. 

5. Because of the effective decoupling of the quarks and therefore of the whole colored 
sector the details of confinement have only little influence on the mesonic dynamics 
for scales 

k< M q ~ 300 MeV . 

Here quarks and gluons disappear effectively from the spectrum and one is left with 
the pions. They are the only particles whose propagation is not suppressed by a large 
mass. For scales below the pion mass the flow of the couplings stops. 

Let us now try to understand these different ranges of scales in more detail. We may 
start at k p — 1.5 GeV where we assume that all gluonic degrees of freedom have been inte- 
grated out while we have kept an effective infrared cutoff ~ k p in the quark propagators. 
Details of this procedure were outlined in ref. [45]. This results in a non-trivial momentum 
dependence of the quark propagator and effective non-local four and higher quark interac- 
tions. Because of the infrared cutoff the resulting effective action for the quarks resembles 
closely the one for heavy quarks (at least for Euclidean momenta). The dominant effect 
is the appearance of an effective quark potential similar to the one for the charm quark 
which describes the effective four quark interactions. 

We next have to remove the infrared cutoff for the quarks, k — > 0. This task can be 
carried out by means of the exact flow equation for quarks only, starting at k p with an 
initial value r fcp [■?/>] as obtained after integrating out the gluons. A first investigation in 
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Figure 19: The left diagram represents the one gluon exchange t-channel contribution to the four quark 
vertex at the scale k p ~ 1.5 GeV. It is assumed here that the gluon propagator is modified such that it 
accounts for the linearly rising term in the heavy quark potential. The right diagram displays the scalar 
meson s-channel exchange found at the compositeness scale fc$ ~ 600 MeV. 

this direction [44] used a truncation with a chirally invariant four quark interaction whose 
most general momentum dependence was retained. A crucial point is, of course, the initial 
value for this momentum dependence at k p . The ansatz used in ref. [44] is obtained by 
Fierz transforming the heavy quark potential and keeping, for simplicity, only the scalar 
meson channel while neglecting the p-meson and pomeron channels which are also present. 
The effective heavy quark potential was approximated there by a one gluon exchange term 
~ a s (k p ) supplemented by a linearly rising string tension term. This ansatz corresponds 
to the four quark interaction generated by the flavor neutral t-channel one gluon exchange 
depicted in figure 19 with an appropriately modified gluon propagator and quark gluon 
vertex in order to account for the linearly rising part of the potential. 

The evolution equation for the four quark interaction can be derived from the fermionic 
version of eq. (2.10). It is by far not clear that the evolution of the effective four quark 
vertex will lead at lower scales to a momentum dependence representing the (s-channel) 
exchange of colorless mesonic bound states. Yet, at the compositeness scale 



one finds [44] an approximate Bethe-Salpeter factorization of the four quark amplitude with 
precisely this property. This situation is described by the right Feynman diagram in figure 
19. In particular, it was possible to extract the amputated Bethe-Salpeter wave function as 
well as the mesonic bound state propagator displaying a pole-like structure in the s-channel 
if it is continued to negative s = (pi + P2) 2 - I n the limit where the momentum dependence 



k<t, ~ 600 MeV 



(5.3) 
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of the Bethe-Salpeter wave function and the bound state propagator is neglected the 
effective action resembles 36 the Nambu-Jona-Lasinio model [123, 124]. It is therefore 
not surprising that our description of the dynamics for k < k$ will parallel certain aspects 
of the investigations of this model, even though we are not bound to the approximations 
used typically in such studies (large-iV c expansion, perturbative renormalization group, 
etc.). 

It is clear that for scales k ^ /c$ a description of strong interaction physics in terms of 
quark fields alone would be rather inefficient. Finding physically reasonable truncations 
of the effective average action should be much easier once composite fields for the mesons 
are introduced. The exact renormalization group equation can indeed be supplemented by 
an exact formalism for the introduction of composite field variables or, more generally, a 
change of variables [44]. In the context of QCD this amounts to the replacement of the 
dominant part of the four quark interactions by scalar meson fields with Yukawa couplings 
to the quarks. In turn, this substitutes the effective quark action at the scale k$ by the 
effective quark meson action given in eq. (5.2) in the introduction 37 . The term in the 
effective potential which is quadratic in $, = ra£ tr<&t<I> + . . ., turns out to be 
positive as a consequence of the attractiveness of the four quark interaction inducing it. 
Its value was found for the simple truncations used in ref. [44] to be ~ 120 MeV. The 
higher order terms in cannot be determined in the four quark approximation since 
they correspond to terms involving six or more quark fields. (Their values will not be 
needed for our quantitative investigations as is discussed in subsection 5.4). The initial 
value of the (bare) Yukawa coupling corresponds to the amputated Bethe-Salpeter wave 
function. Neglecting its momentum dependence it can be normalized to = 1. Moreover, 
the quark wave function renormalization Z^ ;k is normalized to one at the scale k§ for 
convenience. One may add that we have refrained here for simplicity from considering 
four quark operators with vector and pseudo-vector spin structure. Their inclusion is 
straightforward and would lead to vector and pseudo-vector mesons in the effective action. 

In view of the possible large truncation errors made in ref. [44] we will take (5.3) and 
the above value of m kii only as order of magnitude estimates. Furthermore, we will assume, 
as motivated in the introduction and usually done in large-iV c computations within the 
NJL model, that 

Z* M « 1 • (5-4) 

As a consequence, the initial value of the renormalized Yukawa coupling — Z^J^ Z^^ 
is much larger than one and we will be able to exploit the infrared stable features of 
the flow equations. As a typical coupling we take = 100 in order to simulate the limit 
hk iy — > oo. The effective potential £4(<&) must be invariant under the chiral SUl(N) x 
SU R (N) flavor symmetry. In fact, the axial anomaly of QCD breaks the Abelian C/a(1) 

36 0ur solution of the flow equation for Tk with Z&^z = (see below) may be considered as a solution 
of the NJL model with a particular form of the ultraviolet cutoff dictated by the shape of Rk(q 2 ) as given 
in cq. (2.3). 

37 We note that no double counting problem arises in this procedure. 
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symmetry. The resulting Ua(1) violating multi-quark interactions 38 lead to corresponding 
Ua(1) violating terms in Uk{&). Accordingly, the most general effective potential U\. is a 
function of the N + 1 independent C and P conserving SUl{N) x SUr(N) invariants 

p = tr , 

n ~ tr(V$--^p) , i = 2,...,JV, (5.5) 
£ = det $ + det $ f . 

For a given initial form of £4 all quantities in our truncation of Tk (5.2) are now fixed and 
we may follow the flow of T k to k — > 0. In this context it is important that the formalism 
for composite fields [44] also induces an infrared cutoff in the meson propagator. The flow 
equations are therefore exactly of the form (2.10), with quarks and mesons treated on an 
equal footing. At the compositeness scale the quadratic term of = Tr $t$ + ... is 
positive and the minimum of therefore occurs for $ = 0. Spontaneous chiral symmetry 
breaking is described by a non-vanishing expectation value ($) in absence of quark masses. 
This follows from the change of the shape of the effective potential as k flows from k$> 
to zero. The large renormalized Yukawa coupling rapidly drives the scalar mass term to 
negative values and leads to a potential minimum away from the origin at some scale 
^xSB < k§ such that finally (<&) = a ^ for k — > [44, 46]. This concludes our overview 
of the general features of chiral symmetry breaking in the context of flow equations for 
QCD. 

We will concentrate in this work on the two flavor case (N = 2) and comment on the 
effects of including the strange quark in section 5.8. Furthermore we will neglect isospin 
violation and therefore consider a singlet source term j proportional to the average light 
current quark mass m = \{m u + mj. Due to the C/A(l)-anomaly there is a mass split 
for the mesons described by $. The scalar triplet (do) and the pseudoscalar singlet (rj') 
receive a large mass whereas the pseudoscalar triplet (it) and the scalar singlet (a) remains 
light. From the measured values m v /, m ao ~ 1 GeV it is evident that a decoupling of these 
mesons is presumably a very realistic limit 39 . It can be achieved in a chirally invariant 
way and leads to the well known 0(4) symmetric Gell-Mann-Levy linear sigma model [5] 
which is, however, coupled to quarks now. This is the two flavor linear quark meson model 
which we will study in the remainder of this work. For this model the effective potential 
Uk is a function of p only. 

It remains to determine the source j as a function of the average current quark mass 
m. This is carried out in appendix D and we obtain in our normalization with Z^ >k<6 = 1, 
hk<s, = 1, 

j = 1m\jn . (5.6) 

38 A first attempt for the computation of the anomaly term in the fermionic effective average action can 
be found in ref. [125]. 

39 In thermal equilibrium at high temperature this decoupling is not obvious. We will comment on this 
point in section 5.8. 
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It is remarkable that higher order terms do not influence the relation between j and m. 
Only the quadratic term to| enters which is in our scenario the only relevant coupling. 
This feature is an important ingredient for the predictive power of the model as far as the 
absolute size of the current quark mass is concerned. 

The quantities which are directly connected to chiral symmetry breaking depend on 
the /c-dependent expectation value ($) fc = <7 0j fc as given by 

^(P = 2< fe ) = -^. (5.7) 
op 2a 0ik 

In terms of the renormalized expectation value 

°o,fc = Z^la 0jk (5.8) 
we obtain the following expressions for phenomenological observables from (5.2) for 40 d — 4 

fn,k = 2<7 ,£; , 

fa) k = - 2m l* [z^ k /2 (?o,k - rh] , 



v -i/2 ^L m _ 7-1/2 3 

0,k 



Oo,fc 2cr, 



(5.9) 



™2 



2 ry-l/2 m k s , m a ■> 2 

m a,k = Z <3>,k — + 4\k(To, k ■ 

&0,k 

Here we have defined the dimensionless, renormalized couplings 

\ -7-2 9 2 U k _ 2 

" " *> k W (P ~ °' kh (5.10) 
h k = Z^ k /2 Z^ k h k . 

We will mainly be interested in the "physical values" of the quantities (5.9) in the limit 
k — > where the infrared cutoff is removed, i.e. f n = fn,k=o, m 1 — m 1 k =oi e ^ c - We point 
out that the formalism of composite fields provides the link [44] to the chiral condensate 
(ipi/j') since the expectation value a is related to the expectation value of a composite 
quark-ant iquark operator. 



5.4 Flow equations and infrared stability 

At first sight, a reliable computation of seems a very difficult task. Without a 

truncation T k is described by an infinite number of parameters (couplings, wave function 

40 We note that the expressions (5.9) obey the well known Gell-Mann-Oakes-Renner relation m\f^ = 
-2m{#) + 0{rh 2 ) [126]. 
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renormalizations, etc.) as can be seen if is expanded in powers of fields and deriva- 
tives. For instance, the sigma mass is obtained as a zero of the exact inverse propagator, 
linifc^olfc (?) | *=<$), which formally receives contributions from terms in Tk with arbitrar- 
ily high powers of derivatives and the expectation value a . Realistic non-perturbative 
truncations of T k which reduce the problem to a manageable size are crucial. We will 
follow here a twofold strategy: 

• Physical observables like meson masses, decay constants, etc., can be expanded in 
powers of (current) quark masses in a similar way as in chiral perturbation theory [3]. 
To a given finite order of this expansion only a finite number of terms of a simulta- 
neous expansion of T k in powers of derivatives and <£> are required if the expansion 
point is chosen properly. Details of this procedure and some results can be found 
in [104, 6, 7]. 

• Because of an approximate partial IR fixed point behavior of the flow equations in 
the symmetric regime, i.e. for k x sB < k < the values of many parameters of 
Tk for k — > will be almost independent of their initial values at the compositeness 
scale k$. For large enough only a few relevant parameters need to be computed 
accurately from QCD. They can alternatively be determined from phenomenology. 
Because of the present lack of an explicit QCD computation we will pursue the latter 
approach. 

In combination, these two points open the possibility for a perhaps unexpected degree of 
predictive power within the linear quark meson model. We wish to stress, however, that 
a perturbative treatment of the model at hand, e.g., using perturbative RG techniques, 
cannot be expected to yield reliable results. The renormalized Yukawa coupling is very 
large at the scale k$. Even the IR value of h k is still relatively big 

2M a 

h k=0 = — ± ~ 6.5 (5.11) 

and hk increases with k. The dynamics of the linear quark meson model is therefore clearly 
non-perturbative for all scales k < k$. 

We will now turn to the flow equations for the linear quark meson model. We first note 
that the flow equations for Tk and r*,— | J d 4 xtr + &j) are identical. The source term 
therefore does not need to be considered explicitly and only appears in the condition (5.7) 
for (<&). It is convenient to work with dimensionless and renormalized variables therefore 
eliminating all explicit /c-dependence. With 

u(t,p) = k- d U k (p) , p = Z^ k k 2 - d p (5.12) 

and using (5.2) as a first truncation of the effective average action one obtains the flow 
equation (t = ha(k/k&)) 

d 

—u = —du + (d — 2 + r]<s>)pu' 

dt f i , (5-13) 

+ 2v d |3/ V; ^) + l%(u' + 2pu"- 77*) - 2i +l N c l { Q F)d {-ph 2 - ^)| . 
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Here v^ 1 = 2 d+1 n d / 2 T(d/2) and primes denote derivatives with respect to p. The number 
of quark colors is denoted as N c . We will always use in the following iV c = 3. Eq. (5.13) is a 
partial differential equation for the effective potential u(t, p) which has to be supplemented 
by the flow equation for the Yukawa coupling and expressions for the anomalous dimensions 
77$, 77^,. The symbols l d , l^ d denote bosonic and fermionic mass threshold functions, 
respectively, which are defined in appendix A. They describe the decoupling of massive 
modes and provide an important non-perturbative ingredient. For instance, the bosonic 
threshold functions 



n + 6 nfi ^_ lu2n _ d f d d q 1 dR k 

(2n) d Z^ <k dt [P(q 2 ) + k 2 w] 



involve the inverse average propagator P(q 2 ) = q 2 + Z^\R k (q 2 ) where the infrared cutoff is 
manifest. These functions decrease ~ w ~( n+1 ) for Since typically w = M 2 /k 2 with 

M a mass of the model, the main effect of the threshold functions is to cut off fluctuations 
of particles with masses M 2 ^> k 2 . Once the scale k is changed below a certain mass 
threshold, the corresponding particle no longer contributes to the evolution and decouples 
smoothly. 

The dimensionless renormalized expectation value k = 2k 2 ~ d Z$, tk d\ k , with a 0; k the 
/c-dependent VEV of $, may be computed for each k directly from the condition (5.7) 

u'{t,K) = -j=k- d -¥zzf = e g . (5.15) 



Note that k = in the symmetric regime for vanishing source term. Equation (5.15) allows 
us to follow the flow of k according to 

j t K= 7^rx[^- d - 2]e '-4 u ' M } (516) 

with A = u"(t, k). We define the Yukawa coupling for p = k and its flow equation reads [46] 



^-h 2 = (rf-4 + 2r7^ + ?7$)/i 2 -2?j d /i 4 |3/^ B)d (^ 2K ,e s ;?7 V ,,?7 ( i,) 



(5.17) 
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Similarly, the scalar and quark anomalous dimensions are infered from 



77$ = 




$,fc — 



4j 4/tA 2 m2 2 (e s , e g + 2Xn] 77$) 




+ 



+ 




(5.18) 



which is a linear set of equations for the anomalous dimensions. The threshold functions 

4f ■> m «i nii m ^p d an d m if 1 ^2 are a ^ so specified in appendix A. 

The flow equations (5.13), (5.16) — (5.18), constitute a coupled system of ordinary and 
partial differential equations which can be integrated numerically. Here we take the effective 
current quark mass dependence of hk, Z^^ and Z^^ into account by stopping the evolution 
according to eqs. (5.17), (5.18), evaluated for the chiral limit, below the pion mass m^. 
(For details of the algorithm used here see refs. [63, 68].) One finds for d = 4 that chiral 
symmetry breaking indeed occurs for a wide range of initial values of the parameters 
including the presumably realistic case of large renormalized Yukawa coupling and a bare 
mass rf%k 9 of order 100 MeV. Driven by the strong Yukawa coupling, the renormalized mass 
term u'(t, p = 0) decreases rapidly and goes through zero at a scale /c x sb not far below 
Here the system enters the spontaneously broken regime and the effective average potential 
develops an absolute minimum away from the origin. The evolution of of the potential 
minimum a 2 k = k 2 K,/2 turns out to be reasonably stable already before k ~ m n where it 
stops. We take this result as an indication that our truncation of the effective action 
leads at least qualitatively to a satisfactory description of chiral symmetry breaking. The 
reason for the relative stability of the IR behavior of the VEV (and all other couplings) 
is that the quarks acquire a constituent mass M q = ha ~ 300 MeV in the spontaneously 
broken regime. As a consequence they decouple once k becomes smaller than M q and the 
evolution is then dominantly driven by the light Goldstone bosons. This is also important 
for our approximation of neglecting the residual gluonic interactions in the quark sector of 
the model as outlined in subsection 5.3. 

Most importantly, one finds that the system of flow equations exhibits an approximate 
IR fixed point behavior in the symmetric regime [46]. To see this explicitly we study the 
flow equations (5.13), (5.16) — (5.18) subject to the condition (5.4). For the relevant range 
of p both u'(t,p) and u'(t,p) + 2pu"(t,p) are then much larger than ph 2 (t) and we may 
therefore neglect in the flow equations all scalar contributions with threshold functions 
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involving these large masses. This yields the simplified equations (d = 4, v 4 1 = 32tt 2 ) 

ju = -4u + (2 + m )pu'-^-/ F ^( 1 -ph'), 

Jt^ = V * k2m > (5.19) 

N c (F)4/ n x, 2 
V* = (0)Al , 

V<p = o . 

Of course, it should be clear that this approximation is only valid for the initial range of 
running below before the (dimensionless) renormalized scalar mass squared u ; (t,p = 0) 
approaches zero near the chiral symmetry breaking scale. The system (5.19) is exactly 
soluble. Using m 4 F ^ 4 (0) = 1 which holds independently of the choice of the IR cutoff we 
find 

h 2 (t) = Z?(t) = - % 

2 87T 1 (5.20) 
u(t,p) = e^uj^p^) - ^ dre~n^\\h\t)pe^) . 

(The integration over r on the right hand side of the solution for u can be carried out by 
first exchanging it with the one over momentum implicit in the definition of the threshold 
function Z^ 4 (see appendix A).) Here ui(p) = u(0,p) denotes the effective average poten- 
tial at the compositeness scale and h 2 is the initial value of h 2 at i.e. for t — 0. For 
simplicity we will use an expansion of the initial value effective potential ui{p) in powers 
of p around p = 

MP) = E (5-21) 

n=0 n - 

even though this is not essential for the forthcoming reasoning. Expanding also /q^ 4 in 
eq. (5.20) in powers of its argument one finds for n > 2 

^M) _ e 2(n-2 )t ^M N c (-lT(n-l)\ (F) r _ (n _ 2)<] 

For decreasing t — > — oo the initial values vJf^ become rapidly unimportant and /h 2n 
approaches a fixed point. For n = 2, i.e., for the quartic coupling, one finds 

^M^-^Sl (523) 
- 1 i-^?t (5 ' 23) 

leading to a fixed point value (u^ /h 2 )* = 1. As a consequence of this fixed point behavior 
the system looses all its "memory" on the initial values uf"- 2 ^ at the compositeness scale 
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This typically happens before the approximation u'(t, p), u'(t, p) + 2pu"(t, p) ^> ph 2 (t) 
breaks down and the solution (5.20) becomes invalid. Furthermore, the attraction to partial 
infrared fixed points continues also for the range of k where the scalar fluctuations cannot 
be neglected anymore. The initial value for the bare dimensionless mass parameter 



u'j(0) m 2 



A'.T 



hj k, 



(5.24) 



is never negligible. (In fact, using the values for m| and k<$> computed in ref. [44] one 
obtains rr%l^/k$ ~ 0.036.) For large hi (and dropping the constant piece w/(0)) the solution 
(5.20) therefore behaves with growing \t\ as 

, 2M 8tt 2 , 

h 2 (t) ~ -— , (5.25) 

u(t,p) * ^e-^^p-^^dre- 4 ^^^ 2 ^)^). 

In other words, for hi — > oo the IR behavior of the linear quark meson model will depend 
(in addition to the value of the compositeness scale /c$ and the quark mass m) only on one 
parameter, m 2 ,^ . We have numerically verified this feature by starting with different values 

for u^P(0). Indeed, the differences in the physical observables were found to be small. This 
IR stability of the flow equations leads to a perhaps surprising degree of predictive power! 
For definiteness we will perform our numerical analysis of the full system of flow equations 
(5.13), (5.16) — (5.18) with the idealized initial value ui{p) = u'j(0)p in the limit hj — > oo. 
It should be stressed, though, that deviations from this idealization will lead only to small 
numerical deviations in the IR behavior of the linear quark meson model as long as the 
condition (5.4) holds, say for hi ^ 15 [46]. 

With this knowledge at hand we may now fix the remaining three parameters of our 
model, m\ and m by using f w = 92.4 MeV, the pion mass M w = 135 MeV and 
the constituent quark mass M q as phenomeno logical input. Because of the uncertainty 
regarding the precise value of M q we give in table 7 the results for several values of M q . 
The first line of table 7 corresponds to the choice of M q and A/ = u"(k) which we will use 
for the forthcoming analysis of the model at finite temperature. As argued analytically 
above the dependence on the value of A/ is weak for large enough hi as demonstrated 
numerically by the second line. Moreover, we notice that our results, and in particular the 
value of j, are rather insensitive with respect to the precise value of M q . It is remarkable 
that the values for k§ and m k<6 are not very different from those computed in ref. [44]. 
As compared to the analysis of ref. [46] the present truncation of Y k is of a higher level 
of accuracy: We now consider an arbitrary form of the effective average potential instead 
of a polynomial approximation and we have included the pieces in the threshold functions 
which are proportional to the anomalous dimensions. It is encouraging that the results are 
rather robust with respect to these improvements of the truncation. 
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310 
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-(180) 3 
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64.4 
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-(174) 3 
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Table 7: The table shows the dependence on the constituent quark mass M q of the input parameters 
, m| / k% and j as well as some of our "predictions" . The phenomenological input used here besides 
M q is /tt = 92.4 MeV, = 135 MeV. The first line corresponds to the values for M q and A/ used in the 
remainder of this work. The other three lines demonstrate the insensitivity of our results with respect to 
the precise values of these parameters. 



Once the parameters k$, m£ and rh are fixed there are a number of "predictions" of 
the linear meson model which can be compared with the results obtained by other methods 
or direct experimental observation. First of all one may compute the value of rh at a scale 
of 1 GeV which is suitable for comparison with results obtained from chiral perturbation 
theory [127] and sum rules [128]. For this purpose one has to account for the running of 
this quantity with the normalization scale from k& as given in table 7 to the commonly 
used value of 1 GeV: m(lGeV) = A~ 1 rh{k^). A reasonable estimate of the factor A is 
obtained from the three loop running of rh in the MS scheme [128]. For M q ~ 300 MeV 
corresponding to the first two lines in table 7 its value is A ~ 1.3. The results for rh(l GeV) 
are in acceptable agreement with recent results from other methods [127, 128] even though 
they tend to be somewhat larger. Closely related to this is the value of the chiral condensate 



1 GeV) = -A< [UZZ% - 2m\ . (5.26) 

These results are quite non-trivial since not only f w and m£ enter but also the computed 
IR value Z^ k=0 . We emphasize in this context that there may be substantial corrections 
both in the extrapolation from k$ to 1 GeV and in the factor a q (see (D.l)). The latter 
is due to the neglected influence of the strange quark which may be important near k$. 
These uncertainties have only little effect on the physics at lower scales as relevant for our 
analysis of the temperature effects. Only the value of j which is fixed by m n enters here. 

A further more qualitative test concerns the mass of the sigma resonance or radial 
mode whose renormalized mass squared may be computed according to (5.9) in the limit 
k — > 0. From our numerical analysis we obtain Afc=o — 20 which translates into m a ~ 
430 MeV. One should note, though, that this result is presumably not very accurate 
as we have employed in this work the approximation of using the Goldstone boson wave 
function renormalization constant also for the radial mode. Furthermore, the explicit chiral 
symmetry breaking contribution to is certainly underestimated as long as the strange 
quark is neglected. In any case, we observe that the sigma meson is significantly heavier 
than the pions. This is a crucial consistency check for the linear quark meson model. A 
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Figure 20: The plot shows (solid line) and /„■ (dashed line) as functions of the current quark mass 
m in units of the physical value TO p h ys - 

low sigma mass would be in conflict with the numerous successes of chiral perturbation 
theory [3] which requires the decoupling of all modes other than the Goldstone bosons in 
the IR-limit of QCD. The decoupling of the sigma meson is, of course, equivalent to the 
limit A — > oo which formally describes the transition from the linear to the non-linear 
sigma model and which appears to be reasonably well realized by the large IR-values of A 
obtained in our analysis. We also note that the issue of the sigma mass is closely connected 
to the value of f^°\ the value of f n in the chiral limit rh = also given in table 7. To 
lowest order in (f n — f^)/ f n or, equivalently, in rh one has 

/. - /;°> = ^ttH = • (5 - 27) 

A larger value of m a would therefore reduce the difference between f^> and /„■. 

In figure 20 we show the dependence of the pion mass and decay constant on the 
average current quark mass rh. These curves depend very little on the values of the initial 
parameters as demonstrated in table 7 by f(°\ We observe a relatively large difference of 
12 MeV between the pion decay constants at rh — m p h ys and rh — 0. According to (5.27) 
this difference is related to the mass of the sigma particle and will be modified in the three 
flavor case. We will later find that the critical temperature T c for the second order phase 
transition in the chiral limit is almost independent of the initial conditions. The values of 
and T c essentially determine the non-universal amplitudes in the critical scaling region 
(cf. section 5.7). In summary, we find that the behavior of our model for small k is quite 
robust as far as uncertainties in the initial conditions at the scale are concerned. We will 
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see that the difference of observables between non-vanishing and vanishing temperature is 
entirely determined by the flow of couplings in the range < k ^ 3T. 

5.5 Thermal equilibrium and dimensional reduction 

The extension of flow equations to thermal equilibrium situations at non-vanishing temper- 
ature T is straightforward [55]. In the Euclidean formalism non-zero temperature results 
in ( ant i-) periodic boundary conditions for (fermionic) bosonic fields in the Euclidean time 
direction with periodicity 1/T [53]. This leads to the replacement 



I 4^ftf) - T E / SU?o(0 + f) (5-28) 



in the trace in (2.10) when represented as a momentum integration, with a discrete spec- 
trum for the zero component 

(i\ _ / 2l7lT for bosons , ■ 

qo[L) ~\ (21 + 1)ttT for fermions . [b * } 

Hence, for T > a four-dimensional QFT can be interpreted as a three-dimensional model 
with each bosonic or fermionic degree of freedom now coming in an infinite number of copies 
labeled by I G 2Z (Matsubara modes). Each mode acquires an additional temperature 
dependent effective mass term q\{V). ^ n a high temperature situation where all massive 
Matsubara modes decouple from the dynamics of the system one therefore expects to 
observe an effective three-dimensional theory with the bosonic zero modes as the only 
relevant degrees of freedom. In other words, if the characteristic length scale associated 
with the physical system is much larger than the inverse temperature the compactified 
Euclidean "time" dimension cannot be resolved anymore. This phenomenon is known as 
"dimensional reduction" [54]. 

The formalism of the effective average action automatically provides the tools for a 
smooth decoupling of the massive Matsubara modes as the scale k is lowered from k 3> T 
to < T. It therefore allows us to directly link the low-T, four-dimensional QFT to 
the effective three-dimensional high-T theory. The replacement (5.28) in (2.10) manifests 
itself in the flow equations (5.13), (5.16) — (5.18) only through a change to T-dependent 
threshold functions. For instance, the dimensionless functions l^iw; 77$) defined in eq. (5.14) 
are replaced by 

iHw t-t*) = ^-^v r dd ~ lq ( 1 dRk{q2) ) 1 

L ^k^ } - 4 Vd 'tJ^Y-A^ dt )[ P ^) + kM n+l 

(5.30) 

where q 2 = q% + q 2 and q = 2ttIT. A list of the various temperature dependent threshold 
functions appearing in the flow equations can be found in appendix B. There we also 
discuss some subtleties regarding the definition of the Yukawa coupling and the anomalous 
dimensions for T^O. In the limit k ^> T the sum over Matsubara modes approaches the 
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integration over a continuous range of go an d we recover the zero temperature threshold 
function l d (w; 77$). In the opposite limit k C T the massive Matsubara modes (/ 7^ 0) are 
suppressed and we expect to find a d — 1 dimensional behavior of l d . In fact, one obtains 
from (5.30) 

l d n {w,T/k- m ) ~ l d n {w- m ) for 

T (5 31") 

l d n {w,T/k- m ) ~ ^Z^ 1 ^;^*) for T > fc . 
fc ^ 

For our choice of the infrared cutoff function Rk, eq. (2.3), the temperature dependent 
Matsubara modes in l d l (w,T/k;r]^) are exponentially suppressed for T ^ k whereas the 
behavior is more complicated for other threshold functions appearing in the flow equations 
(5.13), (5.16) — (5.18). Nevertheless, all bosonic threshold functions are proportional to 
T/k for T ^> k whereas those with fermionic contributions vanish in this limit 41 . This be- 
havior is demonstrated in figure 21 where we have plotted the quotients lf(w,T/k)/lf(w) 
and l[ F) \w,T/k)/l[ F) \w) of bosonic and fermionic threshold functions, respectively. One 
observes that for k ^> T both threshold functions essentially behave as for zero temper- 
ature. For growing T or decreasing k this changes as more and more Matsubara modes 
decouple until finally all massive modes are suppressed. The bosonic threshold function if 
shows for k <C T the linear dependence on T/k derived in eq. (5.31). In particular, for the 
bosonic excitations the threshold function for tu C 1 can be approximated with reasonable 
accuracy by / 4 (w; 77$) for T/k < 0.25 and by (AT/k)l^(w; 77$) for T/k > 0.25. The fermionic 
threshold function Z^ 4 tends to zero for k T since there is no massless fermionic zero 
mode, i.e. in this limit all fermionic contributions to the flow equations are suppressed. On 
the other hand, the fermions remain quantitatively relevant up to T/k ~ 0.6 because of the 
relatively long tail in figure 21b. The transition from four to three-dimensional threshold 
functions leads to a smooth dimensional reduction as k is lowered from > T to C T! 
Whereas for k ^> T the model is most efficiently described in terms of standard four- 
dimensional fields $ a choice of rescaled three-dimensional variables $3 = <$>/y/T becomes 
better adapted for k <T. Accordingly, for high temperatures one will use a potential 

k k 
uz{t,pz) = yU(t,p) ; p 3 = -p . (5.32) 

In this regime corresponds to the free energy of classical statistics and r fc>0 is a 

classical coarse grained free energy. 

For our numerical calculations at non-vanishing temperature we exploit the discussed 
behavior of the threshold functions by using the zero temperature flow equations in the 
range k > 10T. For smaller values of k we approximate the infinite Matsubara sums 
(cf. eq. (5.30)) by a finite series such that the numerical uncertainty at k — 10T is better 
than 10~ 4 . This approximation becomes exact in the limit k <C 10T. 



41 For the present choice of Rk the temperature dependence of the threshold functions is considerably 
smoother than in ref. [55]. 
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Figure 21: The plot shows the temperature dependence of the bosonic (a) and the fermionic (b) threshold 
functions lf{w, T/k) and l^^^w, T/k), respectively, for different values of the dimensionless mass term w. 
The solid line corresponds to w = whereas the dotted ones correspond to w = 0.1, w = 1 and w = 10 
with decreasing size of the dots. For T ^> k the bosonic threshold function becomes proportional to T/k 
whereas the fermionic one tends to zero. In this range the theory with properly rescaled variables behaves 
as a classical three-dimensional theory. 
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5.6 Equation of state and the chiral phase transition 

In section 5.3 we have considered the relevant fluctuations that contribute to the flow of 
in dependence on the scale k. In a thermal equilibrium situation T k also depends on the 
temperature T and one may ask for the relevance of thermal fluctuations at a given scale k. 
In particular, for not too high values of T (cf. sect. 5.8) the "initial condition" for the 
solution of the flow equations should essentially be independent of temperature. This will 
allow us to fix from phenomenological input at T = and to compute the temperature 
dependent quantities in the infrared (k — > 0). We note that the thermal fluctuations which 
contribute to the r.h.s. of the flow equation for the meson potential (5.13) are effectively 
suppressed for T ^ k/A (cf. section 5.5). Clearly for T ^ k^>/3 temperature effects become 
important at the compositeness scale. We expect the linear quark meson model with a 
compositeness scale k§ ~ 600 MeV to be a valid description for two flavor QCD below a 
temperature of about 42 170 MeV. 

We compute the quantities of interest for temperatures T ^ 170 MeV by solving numer- 
ically the T-dependent version of the flow equations (5.13), (5.16) — (5.18) (cf. section 5.5 
and appendix B) by lowering k from k$ to zero. For this range of temperatures we use the 
initial values as given in the first line of table 7. This corresponds to choosing the zero 
temperature pion mass and the pion decay constant (f^ = 92.4 MeV for m n = 135 MeV) as 
phenomenological input. The only further input is the constituent quark mass M q which 
we vary in the range M q ~ 300 — 350 MeV. We observe only a minor dependence of our 
results on M q for the considered range of values. In particular, the value for the critical 
temperature T c of the model remains almost unaffected by this variation. 

We have plotted in figure 22 the renormalized expectation value 2<x of the scalar field 
as a function of temperature for three different values of the average light current quark 
mass m. (We remind that 2<x (T = 0) = /V.) For m = the order parameter a Q of chiral 
symmetry breaking continuously goes to zero for T — * T c = 100.7 MeV characterizing the 
phase transition to be of second order. The universal behavior of the model for small T — T c 
and small rh is discussed in more detail in section 5.7. We point out that the value of T c 
corresponds to = 80.8 MeV, i.e. the value of the pion decay constant for m = 0, which 
is significantly lower than f n = 92.4 MeV obtained for the realistic value m p h y s- If we would 
fix the value of the pion decay constant to be 92.4 MeV also in the chiral limit (m = 0), 
the value for the critical temperature would raise to 1 15 MeV. The nature of the transition 
changes qualitatively for m^0 where the second order transition is replaced by a smooth 
crossover. The crossover for a realistic m p h ys or m 7r (T = 0) = 135 MeV takes place in a 
temperature range T ~ (120 — 150) MeV. The middle curve in figure 22 corresponds to a 
value of m which is only a tenth of the physical value, leading to a zero temperature pion 
mass m n = 45 MeV. Here the crossover becomes considerably sharper but there remain 
substantial deviations from the chiral limit even for such small quark masses m ~ 1 MeV. 

42 There will be an effective temperature dependence of induced by the fluctuations of other degrees 
of freedom besides the quarks, the pions and the sigma which are taken into account here. We will 
comment on this issue in section 5.8. For realistic three flavor QCD the thermal kaon fluctuations will 
become important for T ^ 170 MeV. 
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The temperature dependence of m n has already been mentioned in the introduction (see 
fig. 18) for the same three values of rh. As expected, the pions behave like true Goldstone 
bosons for rh = 0, i.e. their mass vanishes for T < T c . Interestingly, m n remains almost 
constant as a function of T for T < T c before it starts to increase monotonically. We 
therefore find for two flavors no indication for a substantial decrease of m n around the 
critical temperature. 

The dependence of the mass of the sigma resonance m a on the temperature is displayed 
in figure 23 for the above three values of rh. In the absence of explicit chiral symmetry 
breaking, rh — 0, the sigma mass vanishes for T < T c . For T <T C this is a consequence of 
the presence of massless Goldstone bosons in the Higgs phase which drive the renormalized 
quartic coupling A to zero. In fact, A runs linearly with k for T ^ k/4 and only evolves log- 
arithmically for T ^ k/4. Once rh ^ the pions acquire a mass even in the spontaneously 
broken phase and the evolution of A with k is effectively stopped at k ~ m w . Because of the 
temperature dependence of (Jo,k=o (cf. figure 22) this leads to a monotonically decreasing 
behavior of m a with T for T < T c . This changes into the expected monotonic growth once 
the system reaches the symmetric phase for 43 T > T c . For low enough rh one may use 
the minimum of m a (T) for an alternative definition of the (pseudo-)critical temperature 
denoted as Table 6 in the introduction shows our results for the pseudocritical tem- 
perature for different values of rh or, equivalently, m n (T = 0). For a zero temperature pion 
mass 77V = 135 MeV we find Tj® = 128 MeV. At larger pion masses of about 230 MeV we 
observe no longer a characteristic minimum for m a apart from a very broad, slight dip at 

43 See section 5.4 for a discussion of the zero temperature sigma mass. 
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T ~ 90 MeV. A comparison of our results with lattice data is given in section 5.7. In fig. 24 
we show the renormalized quartic coupling A as a function of temperature for two fixed 
values of the average current quark mass. The upper curve corresponds to the physical 
value of m or, equivalently, m w (T = 0) = 135 MeV whereas the lower curve shows A for 
m — 0. One observes the vanishing of the renormalized quartic coupling in the chiral limit 
for T < T c as discussed above. The renormalized scalar $ 6 self interaction 

Us(T) = Z- ^ U ^ T \ p = 25g(T)) (5.33) 

assumes a small negative value for realistic quark masses in the temperature range T ^ 
35 MeV with 2U 3 (T)<j$(T) ~ -0.5 < A(T) and 2U 3 (T)<j%(T) ~ 8.0,8.5,1.5 for T = 
80, 120, 160 MeV. We display U%(T) in figure 25 for the chiral limit where one observes a 
discontinuity of U 3 (T) at the critical temperature T c . 

Our results for the chiral condensate as a function of temperature for different 

values of the average current quark mass are presented in figure 17 in the introduction. 
We will compare (T,m) with its universal scaling form for the 0(4) Heisenberg 

model in section 5.7. Another interesting quantity is the temperature dependence of the 
constituent quark mass. Figure 26 shows M q (T) for m = 0, m = m phys /10 and m = m phys , 
respectively. Its behavior is related to the temperature dependence of the renormalized 
order parameter a (T) = a 0:k=0 (T) and the renormalized Yukawa coupling h(T) = h k =o(T). 
The temperature dependence of h in the chiral limit can be found in fig. 27. Near the critical 
temperature one notices a characteristic dip. This results from the long wavelength pion 
fluctuations through a non-analytic behavior of the mesonic wave function renormalization 
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Figure 24: The plot shows the renormalized quartic scalar self coupling A as a function of temperature 
T for the physical value of rh (solid line) as well as for rh — (dashed line). 
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Figure 25: The plot shows the renormalized $ 6 scalar self coupling U3 as a function of temperature T 
in the chiral limit. 
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Fi gure 26: The plot shows the constituent quark mass M q as a function of T for three different values 
of the average light current quark mass m. The solid line corresponds to the realistic value m = m p h ys 
whereas the dotted line represents the situation without explicit chiral symmetry breaking, i.e., m = 0. 
The intermediate, dashed line assumes rh = m p h ys /10. 
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Figure 27: The plot shows the Yukawa coupling, h, as a function of temperature T in the chiral limit. 
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Figure 28: The plot shows the scalar (solid line) and quark (dashed line) wave function renormalization 
constants, Z${T) and Z^{T) ■ 10 -2 , respectively, as functions of temperature T for m = 0. 



Zq>(T) = Z<$>^ = q(T) which is displayed in figure 28. There we also present the temperature 
dependence of the fermionic wave function renormalization Z^(T) = Z^^ = q(T). Away from 
the chiral limit we take the effective quark mass dependence of hk(T), Z$ fc (T) and Z^ >k {T) 
into account by stopping their evolution when k reaches the temperature dependent pion 
mass. In this way we observe no substantial quark mass dependence of these quantities 
except for Z$(T), and consequently for h(T), in the vicinity of the critical temperature. A 
more complete truncation would incorporate field dependent wave function renormalization 
constants and a field dependent Yukawa coupling. 

Our ability to compute the complete temperature dependent effective meson potential 
U is demonstrated in fig. 29 where we display the derivative of the potential with respect 
to the renormalized field 0r = (Z<j>p/2) 1//2 , for different values of T. The curves cover a 
temperature range T = (5—175) MeV. The first one to the left corresponds to T = 175 MeV 
and neighboring curves differ in temperature by AT = 10 MeV. One observes only a weak 
dependence of dU (T) / <90_r on the temperature for T ^ 60 MeV. Evaluated for 4> R = a 
this function connects the renormalized field expectation value with m n (T), the source j 
and the mesonic wave function renormalization Z<$,{T) according to 

= *<>) = = ^o(T)ml(T) . (5.34) 

We close this section with a short assessment of the validity of our effective quark 
meson model as an effective description of two flavor QCD at non-vanishing temperature. 
The identification of qualitatively different scale intervals which appear in the context of 
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Figure 29: The plot shows the derivative of the meson potential U(T) with respect to the renormalized 
field 4> R = {Z^p/2) 1 / 2 for different values of T. The first curve on the left corresponds to T = 175 MeV. 
The successive curves to the right differ in temperature by AT = 10 MeV down to T — 5 MeV. 

chiral symmetry breaking, as presented in section 5.3 for the zero temperature case, can 
be generalized to T ^ 0: For scales below k$> there exists a hybrid description in terms of 
quarks and mesons. For k x ss < k ^ 600 MeV chiral symmetry remains unbroken where 
the symmetry breaking scale k x sB(T) decreases with increasing temperature. Also the 
constituent quark mass decreases with T (cf. figure 26). The running Yukawa coupling 
depends only mildly on temperature for T ^ 120 MeV (see fig. 27). (Only near the critical 
temperature and for m = the running is extended because of massless pion fluctuations.) 
On the other hand, for A; < 4T the effective three-dimensional gauge coupling increases 
faster than at T = leading to an increase of Aqcd(^) with T [17]. As A; gets closer to the 
scale Aqcd(7 1 ) it is no longer justified to neglect in the quark sector confinement effects 
which go beyond the dynamics of our present quark meson model. Here it is important 
to note that the quarks remain quantitatively relevant for the evolution of the meson 
degrees of freedom only for scales k > T/0.6 (cf. fig. 21, section 5.5). In the limit k <C 
T/0.6 all fermionic Matsubara modes decouple from the evolution of the meson potential 
according to the temperature dependent version of eq. (5.13). Possible sizeable confinement 
corrections to the meson physics may occur if Aqcd(7") becomes larger than the maximum 
of M q (T) and T/0.6. From fig. 26 we infer that this is particularly dangerous for small rh in 
a temperature interval around T c . Nevertheless, the situation is not dramatically different 
from the zero temperature case since only a relatively small range of k is concerned. We do 
not expect that the neglected QCD non-localities lead to qualitative changes. Quantitative 
modifications, especially for small rh and |T — T c | remain possible. This would only effect 
the non-universal amplitudes (see sect. 5.7). The size of these corrections depends on the 
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strength of (non-local) deviations of the quark propagator and the Yukawa coupling from 
the values computed in the quark meson model. 



5.7 Universal critical equation of state 

In this section we study the linear quark meson model in the vicinity of the critical tem- 
perature T c close to the chiral limit rh — 0. In this region we find that the sigma mass 



in 



-i 



is much larger than the inverse temperature T _1 , and one observes an effectively 
three-dimensional behavior of the high temperature quantum field theory. We also note 
that the fermions are no longer present in the dimensionally reduced system as has been 
discussed in section 5.5. We therefore have to deal with a purely bosonic 0(4)-symmetric 
linear sigma model. At the phase transition the correlation length becomes infinite and 
the effective three-dimensional theory is dominated by classical statistical fluctuations. In 
particular, the critical exponents which describe the singular behavior of various quantities 
near the second order phase transition are those of the corresponding classical system. 

Many properties of this system are universal, i.e. they only depend on its symmetry 
(0(4)), the dimensionality of space (three) and its degrees of freedom (four real scalar 
components). Universality means that the long-range properties of the system do not 
depend on the details of the specific model like its short distance interactions. Neverthe- 
less, important properties as the value of the critical temperature are non-universal. We 
emphasize that although we have to deal with an effectively three-dimensional bosonic 
theory, the non-universal properties of the system crucially depend on the details of the 
four-dimensional theory and, in particular, on the fermions. 

Our aim is a computation of the critical equation of state which relates for arbitrary 
T near T c the derivative of the free energy or effective potential U to the average current 
quark mass rh. The equation of state then permits to study the temperature and quark 
mass dependence of properties of the chiral phase transition. 

At the critical temperature and in the chiral limit there is no scale present in the theory. 
In the vicinity of T c and for small enough rh one therefore expects a scaling behavior of the 
effective average potential Uk [92] and accordingly a universal scaling form of the equation 
of state (cf. section 3.3). There are only two independent scales close to the transition 
point which can be related to the deviation from the critical temperature, T — T c , and 
to the explicit symmetry breaking through the quark mass rh. As a consequence, the 
properly rescaled potential can only depend on one scaling variable. A possible choice for 
the parameterization of the rescaled "unrenormalized" potential is the use of the Widom 
scaling variable [88] 

(T — T c ) /T c 

x ~ ( 2zf /T c r ' (5,35) 

Here (5 is the critical exponent of the order parameter <t in the chiral limit rh = (see 
equation (5.39)). With U'(p = 2<7q) = j/(2Wq) the Widom scaling form of the equation of 
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state reads [88] 




(5.36) 



where the exponent 5 is related to the behavior of the order parameter according to 
(5.41). The equation of state (5.36) is written for convenience directly in terms of four- 
dimensional quantities. They are related to the corresponding effective variables of the 
three-dimensional theory by appropriate powers of T c . The source j is determined by 
the average current quark mass rh as j = 2m| t rh. The mass term at the composite- 
ness scale, fn 2 .^, also relates the chiral condensate to the order parameter according to 



(ipipj = —2m\^{aQ — rh). The critical temperature of the linear quark meson model was 

found in section 5.6 to be T c = 100.7 MeV. 

The scaling function / is universal up to the model specific normalization of x and 
itself. Accordingly, all models in the same universality class can be related by a rescaling 
of cfo and T — T c . The non-universal normalizations for the quark meson model discussed 
here are defined according to 



We find D = 1.82 ■ 10~ 4 , B = 7.41 and our result for f3 is given in table 8. Apart from the 
immediate vicinity of the zero of f(x) we find the following two parameter fit (cf. section 
3.3) for the scaling function, 



to reproduce the numerical results for / and df/dx at the 1—2% level with 9 = 0.625 (0.656), 
A = —0.490 (—0.550) for x > (x < 0) and 7 as given in table 8. The universal properties 
of the scaling function can be compared with results obtained by other methods for the 
three-dimensional 0(4) Heisenberg model. In figure 30 we display our results along with 
those obtained from lattice Monte Carlo simulation [61], second order epsilon expansion 
[60] and mean field theory. We observe a good agreement of average action, lattice and 
epsilon expansion results within a few per cent for T < T c . Above T c the average action 
and the lattice curve go quite close to each other with a substantial deviation from the 
epsilon expansion and mean field scaling function 44 . 

Before we use the scaling function f(x) to discuss the general temperature and quark 
mass dependent case, we consider the limits T = T c and rh — 0, respectively. In these 
limits the behavior of the various quantities is determined solely by critical amplitudes and 
exponents. In the spontaneously broken phase (T < T c ) and in the chiral limit we observe 

44 We note that the question of a better agreement of the curves for T < T c or T > T c depends on the 
chosen non-universal normalization conditions for x and / (cf. eq. (5.37)). 




/(0) = D 



fi-B-W) = . 



(5.37) 



/fit (x) 



1.816 • 10~ 4 (1 + 136.1 x) 2 (1 + 160.9 9x) A 
(1 + 160.9 (0.9446 ^-i/( 7 -2-A) x y-2-A 



(5.38) 
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Figure 30: The figure shows a comparison of our results, denoted by "average action", with results 
of other methods for the scaling function of the three-dimensional 0(4) Heisenberg model. We have 
labeled the axes for convenience in terms of the expectation value Wq and the source j of the corresponding 
four-dimensional theory. The constants B and D specify the non-universal amplitudes of the model 
(cf. eq. 5.37). The curve labeled by "MC" represents a fit to lattice Monte Carlo data. The second order 
epsilon expansion [60] and mean field results are denoted by "e" and "mf", respectively. Apart from our 
results the curves are taken from ref. [61]. 
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V 


7 


S 


P 


V 


average action 


0.787 


1.548 


4.80 


0.407 


0.0344 


FD 


0.73(2) 


1.44(4) 


4.82(5) 


0.38(1) 


0.03(1) 


MC 


0.7479(90) 


1.477(18) 


4.851(22) 


0.3836(46) 


0.0254(38) 



Table 8: The table shows the critical exponents corresponding to the three-dimensional 0(4)-Heisenbcrg 
model. Our results are denoted by "average action" whereas "FD" labels the exponents obtained from 
perturbation series at fixed dimension to seven loops [98]. The bottom line contains lattice Monte Carlo 
results [129]. 



that the renormalized and unrenormalized order parameters scale according to 

2 < T ) , OE . 1/2 (T c -ry/ 2 



, / (5-39) 
2W (T) B (T C -T\ li 



respectively, with E = 0.814 and the value of B given above. In the symmetric phase the 
renormalized mass m = m n = m a and the unrenormalized mass fn = Z^ 2 m behave as 

^ - (eT ( W . 

in (5.40) 

where £ + = 0.270, C + = 2.79. For T = T c and non- vanishing current quark mass we have 

^ = D -Vs(-L\ 1/S (5.41) 

with the value of D given above. 

Though the five amplitudes E, B, £ + , C + and D are not universal there are ratios of 
amplitudes which are invariant under a rescaling of Wq and T — T c . Our results for the 
universal amplitude ratios are 

R x = C + DB & - 1 = 1.02 , 

= (£+)/V"£>V(«+i)£ = 0.852 , (5.42) 
t\ + E = 0.220. 



Those for the critical exponents are given in table 8. Here the value of 77 is obtained 
from the temperature dependent version of eq. (5.18) (cf. appendix B) at the critical 



106 



5 The equation of state for two flavor QCD 



temperature. For comparison table 8 also gives the results from perturbation series at 
fixed dimension to seven-loop order [98, 67] as well as lattice Monte Carlo results [129] 
which have been used for the lattice form of the scaling function in fig. 30 . 45 There are 
only two independent amplitudes and critical exponents, respectively. They are related by 
the usual scaling relations of the three-dimensional scalar 0(iV)-model [67] which we have 
explicitly verified by the independent calculation of our exponents. 



We turn to the discussion of the scaling behavior of the chiral condensate ypipj for the 

general case of a temperature and quark mass dependence. In figure 17 in the introduc- 
tion we have displayed our results for the scaling equation of state in terms of the chiral 
condensate 46 



as a function of T/T c = 1 + x(j/T c 3 /(x)) 1///3<5 for different quark masses or, equivalently, 
different values of j. The curves shown in figure 17 correspond to quark masses rh = 0, rh = 
m p hy S /10, rh = m p hys and rh = 3.5m phys or, equivalently, to zero temperature pion masses 
m n = 0, = 45MeV, m n = 135 MeV and m n = 230 MeV, respectively (cf. figure 20). 
One observes that the second order phase transition with a vanishing order parameter at 
T c for rh = is turned into a smooth crossover in the presence of non-zero quark masses. 

The scaling form (5.43) for the chiral condensate is exact only in the limit T — > T c , 
j — > 0. It is interesting to find the range of temperatures and quark masses for which \^ip} 
approximately shows the scaling behavior (5.43). This can be infered from a comparison 
(see fig. 17) with our full non-universal solution for the T and j dependence of (iptp^ 
as described in section 5.6. For m n = one observes approximate scaling behavior for 
temperatures T <; 90 MeV. This situation persists up to a pion mass of m n = 45 MeV. 
Even for the realistic case, m n = 135 MeV, and to a somewhat lesser extent for m n = 
230 MeV the scaling curve reasonably reflects the physical behavior for T <; T c . For 
temperatures below T c , however, the zero temperature mass scales become important and 
the scaling arguments leading to universality break down. 

The above comparison may help to shed some light on the use of universality arguments 
away from the critical temperature and the chiral limit. One observes that for temperatures 
above T c the scaling assumption leads to quantitatively reasonable results even for a pion 
mass almost twice as large as the physical value. This in turn has been used for two 
flavor lattice QCD as theoretical input to guide extrapolation of results to light current 
quark masses. From simulations based on a range of pion masses 0.3 ^ m n /m p ^ 0.7 and 
temperatures < T ^ 250 MeV a "pseudocritical temperature" of approximately 140 MeV 
with a weak quark mass dependence is reported [36] . Here the "pseudocritical temperature" 
T pc is defined as the inflection point of (ipip) as a function of temperature. The values of 



45 See also ref. [95] and references therein for a recent calculation of critical exponents using similar 
methods as in this work. For high precision estimates of the critical exponents see also ref. [89, 90]. 





(5.43) 
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the lattice action parameters used in [36] with N t = 6 were am = 0.0125, Q/g 2 = 5.415 and 
am = 0.025, Q/g 2 = 5.445. For comparison with lattice data we have displayed in figure 17 
the temperature dependence of the chiral condensate for a pion mass m w = 230 MeV. From 
the free energy of the linear quark meson model we obtain in this case a pseudocritical 
temperature of about 150 MeV in reasonable agreement with the results of ref. [36]. In 
contrast, for the critical temperature in the chiral limit we obtain T c = 100.7 MeV. This 
value is considerably smaller than the lattice results of about (140 — 150) MeV obtained 
by extrapolating to zero quark mass in ref. [36]. We point out that for pion masses as 
large as 230 MeV the condensate (ipip\ (T) is almost linear around the inflection point for 
quite a large range of temperature. This makes a precise determination of T c somewhat 
difficult. Furthermore, figure 17 shows that the scaling form of (^pip^ (T) underestimates 
the slope of the physical curve. Used as a fit with T c as a parameter this can lead to an 
overestimate of the pseudocritical temperature in the chiral limit. We also mention here the 
results of ref. [130]. There two values of the pseudocritical temperature, T pc = 150(9) MeV 
and T pc = 140(8), corresponding to am = 0.0125, Q/g 2 = 5.54(2) and am = 0.00625, 
Q/g 2 = 5.49(2), respectively, (both for N t = 8) were computed. These values show a 
somewhat stronger quark mass dependence of T pc and were used for a linear extrapolation 
to the chiral limit yielding T c = 128(9) MeV. 

The linear quark meson model exhibits a second order phase transition for two quark 
flavors in the chiral limit. As a consequence the model predicts a scaling behavior near 
the critical temperature and the chiral limit which can, in principle, be tested in lattice 
simulations. For the quark masses used in the present lattice studies the order and univer- 
sality class of the transition in two flavor QCD remain a partially open question. Though 
there are results from the lattice giving support for critical scaling [131, 132] there are also 
recent simulations with two flavors that reveal significant finite size effects and problems 
with 0(4) scaling [133, 134]. 

5.8 Conclusions 

So far we have investigated the chiral phase transition of QCD as described by the linear 
0(4)-model containing the three pions and the sigma resonance as well as the up and 
down quarks as degrees of freedom. Of course, it is clear that the spectrum of QCD is 
much richer than the states incorporated in our model. It is therefore important to ask to 
what extent the neglected degrees of freedom like the strange quark, strange (pseudo) scalar 
mesons, (axial) vector mesons, baryons, etc., might be important for the chiral dynamics 
of QCD. Before doing so it is perhaps instructive to first look into the opposite direction 
and investigate the difference between the linear quark meson model described here and 
chiral perturbation theory based on the non-linear sigma model [3]. In some sense, chiral 
perturbation theory is the minimal model of chiral symmetry breaking containing only the 
Goldstone degrees of freedom. By construction it is therefore only valid in the sponta- 
neously broken phase and can not be expected to yield realistic results for temperatures 
close to T c or for the symmetric phase. However, for small temperatures (and momentum 
scales) the non-linear model is expected to describe the low-energy and low-temperature 
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limit of QCD reliably as an expansion in powers of the light quark masses. For vanishing 
temperature it has been demonstrated recently [6, 7] that the results of chiral perturbation 
theory can be reproduced within the linear meson model once certain higher dimensional 
operators in its effective action are taken into account for the three flavor case. More- 
over, some of the parameters of chiral perturbation theory (L 4 , . . . , L 8 ) can be expressed 
and therefore also numerically computed in terms of those of the linear model. For non- 
vanishing temperature one expects agreement only for low T whereas deviations from chiral 
perturbation theory should become large close to T c . Yet, even for T < T c small quan- 
titative deviations should exist because of the contributions of (constituent) quark and 
sigma meson fluctuations in the linear model which are not taken into account in chiral 
perturbation theory. 

From [135] we infer the three-loop result for the temperature dependence of the chiral 
condensate in the chiral limit for N light flavors 

!jV 2 — IT 2 N 2 — 1 / T 2 \ 2 

3 ' ( 5 ' 44) 

+ ln£} +0 (T 8 ). 

The scale T\ can be determined from the .D-wave isospin zero tvtt scattering length and 
is given by Ti = (470 ± 100) MeV. The constant Fq is (in the chiral limit) identical to 
the pion decay constant F = ffl = 80.8 MeV. In figure 31 we have plotted the chiral 
condensate as a function of T/F for both, chiral perturbation theory according to (5.44) 
and for the linear quark meson model. As expected the agreement for small T is very 
good. Nevertheless, the anticipated small numerical deviations present even for T T c 
due to quark and sigma meson loop contributions are manifest. For larger values of T, 
say for T <; 0.8 ffl the deviations become significant because of the intrinsic inability of 
chiral perturbation theory to correctly reproduce the critical behavior of the system near 
its second order phase transition. 

Within the language of chiral perturbation theory the neglected effects of thermal quark 
fluctuations may be described by an effective temperature dependence of the parameter 
Fq(T). We notice that the temperature at which these corrections become important equals 
approximately one third of the constituent quark mass M q (T) or the sigma mass m a (T), 
respectively, in perfect agreement with fig. 21. As suggested by this figure the onset of the 
effects from thermal fluctuations of heavy particles with a T-dependent mass m^(T) is 
rather sudden for T ^ m^(T)/3. These considerations also apply to our two flavor quark 
meson model. Within full QCD we expect temperature dependent initial values at kq>. 

The dominant contribution to the temperature dependence of the initial values presum- 
ably arises from the influence of the mesons containing strange quarks as well as the strange 
quark itself. Here the quantity m| seems to be the most important one. (The temperature 
dependence of higher couplings like A(T) is not very relevant if the IR attractive behavior 
remains valid, i.e. if remains small for the range of temperatures considered. We 
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Figure 31: The plot displays the chiral condensate (ipip) as a function of T/fi°\ The solid line corre- 
sponds to our results for vanishing average current quark mass m = whereas the dashed line shows the 
corresponding three-loop chiral perturbation theory result for Ti = 470 MeV. 

neglect a possible T-dependence of the current quark mass to.) In particular, for three 
flavors the potential contains a term 

-^fc* ( det $ + det $t ) = -Vk»<P,$uu&M + ■■■ (5-45) 

which reflects the axial Ua(X) anomaly. It yields a contribution to the effective mass term 
proportional to the expectation value ($ ss ) = <p s , i-e. 

A ™L = -\ v k*<Ps • (5.46) 

Both, X7 fe<t and </? s , depend on T. We expect these corrections to become relevant only for 
temperatures exceeding ttik{T)/3 or M s (T)/3. We note that the temperature dependent 
kaon and strange quark masses, to^(T) and M S (T), respectively, may be somewhat differ- 
ent from their zero temperature values but we do not expect them to be much smaller. A 
typical value for these scales is around 500 MeV. Correspondingly, the thermal fluctuations 
neglected in our model should become important for T ^ 170 MeV. It is even conceivable 
that a discontinuity appears in tp s (T) for sufficiently high T (say T ~ 170 MeV). This 
would be reflected by a discontinuity in the initial values of the 0(4)-model leading to a 
first order transition within this model. Obviously, these questions should be addressed 
in the framework of the three flavor SU L (3) x SU R (3) quark meson model. Work in this 
direction is in progress. 

We note that the temperature dependence of V(T)ip s (T) is closely related to the question 
of an effective high temperature restoration of the axial Ua(X) symmetry [31, 32]. The r( 
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mass term is directly proportional to this combination, mi,{T) — m%(T) ~ ^T/(T)tp s (T) 
[104]. Approximate £7a(1) restoration would occur if tp s (T) or v(T) would decrease sizeable 
for large T. For realistic QCD this question should be addressed by a three flavor study. 
Within two flavor QCD the combination Vkip s is replaced by an effective anomalous mass 
term . The temperature dependence of (T) could be studied by introducing quarks 
and the axial anomaly in the two flavor matrix model of ref . [68] . We add that this question 
has also been studied within full two flavor QCD in lattice simulations [133, 136, 137]. So 
far there does not seem to be much evidence for a restoration of the Ua(1) symmetry near 
T c but no final conclusion can be drawn yet. 

To summarize, we have found that the effective two flavor quark meson model presum- 
ably gives a good description of the temperature effects in two flavor QCD for a temperature 
range T ^ 170 MeV. Its reliability should be best for low temperature where our results 
agree with chiral perturbation theory. However, the range of validity is considerably ex- 
tended as compared to chiral perturbation theory and includes, in particular, the critical 
temperature of the second order phase transition in the chiral limit. We have explicitly 
connected the universal critical behavior for small \T — T c \ and small current quark masses 
with the renormalized couplings at T = and realistic quark masses. The main quanti- 
tative uncertainties from neglected fluctuations presumably concern the values of ffl and 
T c which, in turn, influence the non-universal amplitudes B and D in the critical region. 
We believe that our overall picture is rather solid. Where applicable our results compare 
well with numerical simulations of full two flavor QCD. 
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6 Outlook 

The use of non-perturbative flow equations within the framework of the effective average 
action provides a practical means for non-perturbative studies within a field theoretic 
context. In this work its application culminated in the computation of the equation of state 
for two flavor QCD within the effective quark meson model. In particular, the method is 
free of infrared divergences and allows a precise description of critical phenomena in high 
temperature quantum field theory or in statistical physics. 

It is a valuable merit of the approach presented in this work to yield directly the whole 
equation of state of the theory under investigation. We have given a detailed quantitative 
picture of the non-analytic properties of the effective potential or the free energy for second 
order phase transitions. We have demonstrated that along the lines presented in section 
4 also the universal equation of state for weak first order phase transitions can be reliably 
computed. Furthermore, for k > the effective average action realizes the concept of a 
coarse grained free energy. We, therefore, were able to extract quantities like the "classical" 
interface tension used in the standard approach to phase conversion by spontaneous bubble 
nucleation. 

Apart from their phenomenological relevance and their intrinsic appeal, critical phe- 
nomena also provide a severe test for the reliability of a non-perturbative method. The 
present level of truncation used in this work is based on a lowest order derivative expansion 
and the most general potential term 47 . Where available we observe agreement with recent 
lattice results for the critical equation of state within a few per cent accuracy. So far the 
main drawback of the method is a precise error estimate which is a natural difficulty in a 
non-perturbative calculation with couplings that cannot be considered as small. However, 
the structure of the exact flow equation (2.10) is of considerable help in this context. For 
instance, for the O(N) symmetric iV-component model one easily verifies that the ansatz 
(cf. section 3 for notation) 

Y k = j d d x{u k (p) + ^<f> a Z k (p, -d»d v )d^ a + ^pY k (p, -d»d u )d,p}. (6.1) 

contains the most general terms which contribute to Y k when evaluated for a constant 
background field [92]. With the above ansatz one therefore obtains an exact flow equa- 
tion for the effective average potential U k containing the unknown functions Z k (p,q 2 ) and 
Y k (p,q 2 ), which play the role of field and momentum dependent wave function renormal- 
ization constants. Knowledge of these functions would yield an exact critical equation 
of state from the solution of the flow equation for U k . These functions involve, however, 
more complicated momentum dependences of 1PI vertices with three or four legs with 
non-vanishing momenta. Along the lines presented in section 3 recent calculations for 
the critical equation of state (N = 1) have been successively extended to include a wave 
function renormalization with a most general field dependence 48 [96]. These results seem 

47 For the scalar matrix models discussed in section 4 we kept the most general form of the potential 
only for the relevant p dependence. 
48 For N = 1 the Y k term is absent. 
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to improve the exponents and the equation of state within the expected accuracy of a few 
per cent. In conjunction with the smallness of the anomalous dimension (rj ~ 0.04) one is 
lead to expect only small corrections from the momentum dependence of the wave function 
renormalization. A crucial test for the stability of the results would be to change a given 
parameterization of the momentum dependence of the wave function renormalization in a 
reasonable range 49 . In this way we hope to obtain reliable error estimates for the critical 
equation of state, critical exponents etc. 

We also note that the approach is valid for arbitrary dimensions d. All evolution equa- 
tions presented in this work are formulated for general d. In particular, two dimensional 
field theories can be treated with this method. This allows, e.g., the comparison with 
known exact results. We have calculated the critical properties for the two dimensional 
O(N) model in the "replica limit" [138] N — > 0, which describes the critical swelling of 
long polymer chains in two space dimensions [139]. For instance, we find the exponent v 
(cf. section 3) to be given by v = 0.782 [140] in the two dimensional model which compares 
rather well to the exact result ^( ex act) = 0.75 [67]. We also note, e.g., that Grater and Wet- 
terich [141] have used the effective average action method in two dimensions to describe 
the Kosterlitz-Thouless phase transition by employing an expansion around the minimum 
of the effective average potential. The validity of the method for arbitrary dimensions d 
may also be used to get a deeper insight into the e-expansion approach [56]. 

The investigation of three dimensional scalar models for complex 2x2 matrices in 
section 4 can be extended in several directions. The generalization to complex N x N 
matrices for arbitrary N is straightforward. This opens the possibility to study the large- 
N behavior and to compare with 1/iV-expansions 50 . The most interesting generalizations 
in the context of chiral symmetry breaking are the systems with reduced SU (N) x SU (N) 
symmetry (cf. section 5). Another interesting application concerns Anderson localization. 
The problem of an electron in a random potential can be mapped onto a field theoretic 
matrix model with e.g. U(2N) symmetry [101]. To get insight into the involved replica 
limit in the context of flow equations, we have computed the free energy for the three 
dimensional O(N) model in the limit N — > [140] (cf. also section 3). Further studies 
are necessary to ensure the validity of this limit for matrix models within the approach 
presented in this work. 

Of course the approach is not restricted to the study of the universal behavior. The 
ability of the method to cover both the non-universal properties of the four dimensional 
theory at low temperature and the universal properties near T c and the chiral limit has 
been crucial for the investigation of the QCD chiral phase transition in section 5. We 
expect our results within the effective quark meson model to give a reliable description 
of the equation of state for two flavor QCD. The employed method can be generalized to 
the (more) realistic description including the strange quark. The relevant generalizations 

49 One expects a very mild momentum dependence of Zk(p = po,q 2 ) for q 2 between zero and k 2 due to 
the presence of the infrared cutoff. For q 2 3> k 2 the anomalous dimension should determine the behavior 
near the transition, Zk{p = po,q 2 ) ~ (q 2 )~ 7 ^ 2 - 

50 Sec refs. [142, 143] for a renormalization group study within the large-A^ limit of hermitean matrix 
models. 
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have been outlined in section 5 and the investigation is work in progress [74]. We expect 
our approach to be very well suited to the non-perturbative investigation of the three 
flavor case. At non-zero temperature we should be able to reliably answer a number 
of urging questions, like the order of the phase transition or if one has to deal with a 
crossover phenomenon. Connected to this, we are well prepared to answer the question 
of a large scale that could lead to a disoriented chiral condensate (cf. section 1). The 
three flavor results will also give details about the speculative "effective restoration" of 
the axial Ua{1) around T c . If the approximate partial infrared fixed point, which we have 
observed for the two flavor case, is realized for Nf = 3 our equation of state will develop a 
considerable degree of predictive power. Of course, this concerns also the zero temperature 
low-energy properties of QCD. In particular, these results will be very interesting for the 
scalar spectrum which is experimentally not very well understood. Further investigations, 
extending the description in terms of quarks and the low-lying scalar and pseudoscalar 
mesons, have to include fields for the lightest vector and pseudovector mesons. 

We finally note that a computation of the equation of state within QCD, i.e. as a 
function of a s (lGeV) and the quark masses, depends on a reliable computation of the 
scalar mass term m 2 : (k$) / 'k% and the Yukawa coupling h(k$) at the compositeness scale. 
Integrating out the gluons completely for the determination of the quark meson model at 
k$ is, however, a difficult task. Present investigations in this direction within the QCD 
framework for the effective average action proposed in ref. [45] confirm the relatively simple 
ansatz at kq> for the resulting quark propagator we have used in section 5 [144] . 
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A Threshold functions 

In this appendix we list the various definitions of dimensionless threshold functions ap- 
pearing in the flow equations and the expressions for the anomalous dimensions for T = 0. 
They involve the inverse scalar average propagator 51 P(q) 

and the corresponding fermionic function P F which can be chosen as [46] 

P F ( q ) = P( q )= q *(l + r F (q)f . (A.2) 

We abbreviate 

F <*> ^ • s s p w • w t p s rlt p • (A - 3) 



etc., and use the formal definition 

9 _ 1 dR k d 2 P F d[Z^ k r F \ d 



dt Z $ifc dt OP Z^l + rp dt dP F 
The bosonic threshold functions read 



(A.4) 



Jo \Z^ k dt M > 



2 

= -h 2(n ^- d dxx^^^P + Wl k 2 Y^ (P + W 2 k'< 



-n 2 



(A.5) 

where n,rii,n 2 > is assumed. For n^O the functions l d may also be written as 

l d n (w- m ) = --k 2n ~ d I™ dxxi- l ^-(P + wk 2 )~ n . (A.6) 

Zi JO Uv 

The fermionic integrals l^ d (w; rj^) = l^ d (w) — r]^l^ d {w) are defined analogously as 



= (» + s n , ) k 2 - d r dxxi- 1 ^- p * 8 [ y (p + w ey (n+1) 

Jo Z/^k J- + Tp at v ' 



(A.7) 



51 We note that the definition of Pk(q) in section 4.4 obtains from the definition used here by a rescaling 
with the factor Zk- 
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Furthermore one has 

€2>i> ^ Vl>, V*) = G\ d (wi, w 2 ) - r^ffK, w 2 ) - WKVi, w 2 ) 

= _}_ k 2( ni +n 2 )-d [°° J TT i'l^_\ \ 

2 Jo dt\[P F (x)+k 2 w 1 ] n ^[P(x)+k 2 w 2 ] n - 

^i una (wi,w 2 ;ri 9 ) = m* u ^(w 1 ,w 2 )-ri*m* ltna (w 1 ,W2) 

2 Jo dt\ [P(x) + fc 2 wj ni [P(x) + k 2 w 2 } n2 

(F)d( \ (F)d, \ . x 



"CVi, ^2; = W2) - ^>™S2Vi ? w 2 ) - ri9m£3?(w 1 ,w 2 ) 

= _l A ,2(n 1+ n 2 -l)- d rdxxi-{ 1 + ^ ^ ] 

2 7o 9t \ [Pf(x) + A; 2 W!]™i [P(x) + £; 2 w 2 ]™ 2 J 

(A.8) 

The dependence of the threshold functions on the anomalous dimensions arises from the 
t-derivative acting on and within and Z^^f, respectively. We furthermore 
use the abbreviations 

£(ifc) = J*(0;ifr) , e )d (%)^4 F)d (0;%) 
^W^feo) , tf = £(0;0) 

etc. and note that in four dimensions the integrals 

4(0, 0) = if )4 (0, 0) = /g B)4 (0, 0) = mf )4 (0) = mf 2 B)4 (0, 0) = 1 (A.10) 
are independent of the particular choice of the infrared cutoff. 



B Temperature dependent threshold functions 

Non-vanishing temperature manifests itself in the flow equations (5.13), (5.16) — (5.18) only 
through a change to T-dependent threshold functions. In this appendix we will define these 
functions and discuss some subtleties regarding the definition of the anomalous dimensions 
and the Yukawa coupling for T^0. The corresponding T = threshold functions can be 
found in appendix A where also some of our notation is fixed. 

The flow equation (5.13) for the effective average potential involves a bosonic and a 
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B Temperature dependent threshold functions 



fermionic threshold function whose generalization to finite temperature is straightforward 



\n+l 



Vd 



x TV / dxx^Z^ PfM 



re 

Jo 



where T = T/k and 



,tz 7o ~~ " 1 + Mvf) [P F (y F ) + k 2 w] n+1 

y = x + (2/ttT) 2 , 



(B.2) 



y F = x + (2/ + l)VT 2 , 

The computation of the anomalous dimensions ?7$, 77^, and the flow equation for the 
Yukawa coupling h at non-vanishing temperature requires some care. The anomalous di- 
mensions determine the IR cutoff scale dependence of Zq>^ and Z^^ according to 77$ = 
— d t In Z^^, r]^ = —dflnZ^k with t = \nk/kq>. It is important to realize that for a compu- 
tation of the scale dependence of the effective three-dimensional Zq>^ and Z^^ for T ^ 
momentum dependent wave function renormalization constants of the four-dimensional 
theory are required. This is a consequence of the fact that in the three-dimensional model 
each of the infinite number of different Matsubara modes of a four-dimensional bosonic or 
fermionic field <p{Q) corresponds to a different value of Qo = 2-kIT or Q = (21 + 1)ttT, 
respectively, with Q 2 = Qq + Q 2 and I e ZZ. We will therefore allow for momentum 
dependent wave function renormalizations, i.e. for a kinetic part of 1^ of the form 

t ** = / ^{^V^^ (B.3) 
in momentum space. 

In the 0(4)-model the evolution equation for Zq>^{Q) may then be obtained by con- 
sidering a background field configuration with a small momentum dependence, 

(x) = ip5 n + (5 V e- lQx + 5<p*e iQx ) 5 j2 ; j = 1 , . . . , 4 (B.4) 

supplemented by 

^ = ^ = 0; a = 1,2. (B.5) 

Expanding around this configuration at the minimum of the effective average potential £4 
we observe that 5ip corresponds to an excitation in the Goldstone boson direction. The 
exact inverse two-point function turns out to be block-diagonal with respect to scalar 
and fermion indices for this configuration. It therefore decays into corresponding matrices 
T^l and T^l acting in the scalar and fermion subspaces, respectively. The scale dependence 
of the scalar wave function renormalization for non-vanishing T is obtained from (2.10) 
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and (B.3) for the configuration (B.4) as 



— Z<$,k{Q) = -=!-[ li m ttt— ; — ri-Tr 

dt * Q2 I s <pS<p'^o 5{5<^S<f*) 1 2 



Tr 



(2) 



(B.6) 



In the three-dimensional theory there is now a different scalar wave function renormaliza- 
tion Z$ t k{Qo, Q) for each Matsubara mode Qq. As in the four-dimensional model for T = 
we neglect the momentum dependence of the wave function renormalization constants and 
evaluate Z<$>^ for Q = for each Matsubara mode. We will furthermore simplify the trun- 
cation of the effective average action by choosing the Matsubara zero-mode wave function 
renormalization constant for all Matsubara modes, i.e., approximate 



Z^ k (f) = Z^ k (Ql = 0,Q 2 = 0) . 



(B.7) 



This is justified by the rapid decoupling of all massive Matsubara modes within a small 
range of k for fixed T as discussed in section 5.5. This results in the expression (5.18) for 
rjq> but now with temperature dependent threshold functions (T = T/k) 



1 



2(ni+n 2 -l)-d+l 

2 (d-l)v d 



dv d -i 



P{y) 



P{y) 



dt \ [P(y) + fc 2 wi] ni [P(y) + k 2 w 2 ] n2 



mf )d (w,T;r)ip) = m [ { >a (w,T) - r] i ,rh\ r >a (w,T) 



AF)d 



(F)d, 



(d - l)v d 



TJ2 



d=i d ( d 1 +r F (y F ) 
axx 2 yFTr 



dt \dx Pf(vf) + k 2 w / 



(B.8) 

For further technical details we refer the reader to ref . [] . 

The fermion anomalous dimension and the flow equation for the Yukawa coupling can 
be obtained by considering a field configuration 



<p8ji ; j = 



1,...,4, 



1,2 . 



(B.9) 



The derivation follows similar lines as for the scalar anomalous dimension discussed above. 
For computational details we refer the reader to ref. [46]. An important difference as 
compared to Z$ ifc relates to the fact that there are no fermionic zero modes. It would 
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therefore be inconsistent to define Z^^(T) or h k (f) at Qq = if Q denotes the external 
fermion momentum. Yet, we will again resort to the approximation of using the same 
wave function renormalization constant and Yukawa coupling for all fermionic Matsubara 
modes. For the same reason as for Z^ >k (T) we will use for this purpose the mode with the 
lowest T-dependent mass, i.e. define 

Z^,k{T) = Z^ k {Ql = n 2 T 2 ,Q 2 = 0) , 

(-D.10) 

h k (T) = h k (Q 2 = n 2 T 2 ,Q 2 = 0), 

where we have neglected a possible dependence of h k on the external scalar momentum of 
the Yukawa vertex. This yields the expressions (5.17) and (5.18) for the flow of h 2 and r]^, 
respectively, but now with the T-dependent threshold functions 

m[ F 2 B)d (w l ,w 2 ,f;ri i ,,ri< s> ) = m ( { 2 B)d (w 1 ,w 2 ,f) 

~ (FB)d( rf,\ ~ (FB)d, 

- r]^m\ )2 (w 1 ,w 2 ,T) -r)^m\ 2 (w 1 ,w 2 ,T) 

2 (d-l)v d ££Jo 

9/ l + r F (y F ) P{y) \ 

x ■ 



dt { [P F {y F ) + k 2 Wi} ni [P(y) + k 2 w 2 ] n * J ' 

- viGl d (vi, w 2 , f) - vJZT^ f ) 

2 v d Jo 



(B.ll) 



d_ r 1 

di \ [P F {y F ) + k 2 Wl }^[P{y) + k 2 w 2 } n * 

O L J. LJ 1 i UJ I « J ! | | u ■ jd 

The integrals 



C Pole structure of the 1% integrals 



l» = -nfdyy^^iyil + riy^+u;}-^, 

Tl f°° d 

#M = 2J0 d yy* r (y)ly( 1 + r (y)) 



(C.l) 



with 

r( y ) - e ~ V dr ^ - e ' V (C 2) 

exhibit for d < 2(n + 1) a singularity at uo = —1. The massless dimensionless average 
propagator y{l + r(y)) is a monotonic function of y that takes on its minimum at y = 
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with lim y _+ y(l + r(y)) = 1. We define new variables 5 and z, 

5 = uj + 1 , z = y(l + r(y))-l (C.3) 



and substitute in (C.l), 



poo 

l d (S) = / dzG d (z)(z + 5)'^ +1 \ 
Jo 

l d n (5) = / dzG d n (z)(z + S)~^ 
Jo 



(C.4) 



with 



G*{z) 



G d n {z) 



±+i dr (y) 
n v 2 o — 



1 • 'I'/) • / 



d V (C.5) 
ny*r(y) 



and y = y(z). For d < 2{n + 1) the integrals (C.4) have a pole at 5 — 0. (The singularity 
becomes logarithmic in 5 for d = 2(n+l)). In this case for 5^0 the dominant contribution 
to the integral comes from the region y ~ or equivalently z ~ 0. To find an approximate 
expression for l d and l d near the pole we expand the regular part of G d and G d around 
z = 0. With 

r(y) = l -( l -\y + -&y 2 + °(y A ))i 



= -7? U- +°(y ) 



(C.6) 



one obtains 



dy y 2 V 12' 

G*(z) = n 2/ ^ 1 (2-^ + ly 2 + 0(y 3 )), 
Gfc) = ^- 1 (2-^ + gy 2 + 0(y 3 )). 



(C.7) 



The inversion of z(y) given by (C.3) can be done by expanding 

z = \y+Y 2 y 2 + °(y A ) ( c - 8 ) 

and we find 

y = 2z- 2 -z 2 + ^ + 0{z i ). (C.9) 
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D The quark mass term 



Insertion of (C.9) in (C.7) yields 

1 (d 
3 I 2 



+ 1 \z + 



G d n (z) = 2n(2z)^ 1 (l 
G d n {z) = n{2z)^(l- l -(± + ^z + 



z 2 + 0(z 3 ) 
z 2 + 0(z 3 ) 



(CIO) 

We consider for d = 3 and n > 1 the zeroth order expression for and that obtains from 
the first term in (CIO) and the exchange of summation and integration in (C.4). Near the 
pole one finds 

l 3 n (5) ~ 2 3 / 2 n / dzz 1 / 2 (z + 5)~^ 

Jo (C.ll) 

SO*) = 2 Z M- 

The leading contributions to Zf, ^ and /f are therefore given by 

11(5) = 2 1 / 2 n5~ 1 / 2 , 



11(5) = 2- 1 / 2 n5- 3 ' 2 , 
11(5) = 2~ 5 / 2 3n5- 5 / 2 . 



(C.12) 



We have verified this numerically. 



D The quark mass term 

In this appendix we determine the source j = diag(j u , j^, . . .) as a function of the average 
current quark mass m. In this context it is important to note that the source j does not 
depend on the IR cutoff scale k. Since j is determined by the properties of the quark 
meson model at the compositeness scale and also enters directly the value of the pion 
mass, which is determined at k — 0, this relation provides a bridge between the short and 
long distance properties of the quark meson model. This will allow us to compute the 
chiral condensate ^VV') or the parameter B of chiral perturbation theory [3]. (We expect, 
however, sizeable corrections when going from two to three flavors. They arise because of 
the relevance of strange quark physics at scales near kq> ). In a more general context we 
need the proportionality coefficient a q between the source j q and the current quark mass 
m q , q — u,d, . . ., taken at the renormalization scale 52 \i = k&, 

J q = ^=^a q m q . (D.l) 

52 We will occasionally use the notation m(/u), m q {p) or (ipip) (/j,) = (tpip) fc _ (/x) (not to be confused 
with (i>ip) k = (i } i') k (a 1 = fc*)) in order to indicate the renormalization scale \x. If no argument is given 
H = fc$ is assumed. 
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For a computation of the coefficient a q we need to look into the details of the introduction 
of composite meson fields in QCD [44, 122]. Let us assume that at the scale kq> a part of 
the QCD average action for quarks [ip] factorizes in the quark bilinear 



/dp 

such that 

r k M = -FkM + r' k M- ( D - 3 ) 

We can then introduce meson fields by inserting the identity 

N J £>$exp(-F fc , + = 1 (D.4) 

into the path integral which formally defines r^[V']- (Here AMs a field independent nor- 
malization factor.) This effectively replaces in (D.3) the term —F kiS> [x] by 53 

-FM + ^ + *] = °^* M + + (D ,) 

The original multi-quark interaction — [x] is canceled by the lowest order term in the 
Taylor expansion in $. Instead, we have substituted mesonic self-interactions i*fc$[3>] and 
interactions between mesons and quarks corresponding to the terms in the expansion which 
contain powers of x an d In particular, we may specialize to the case where the derivative 
terms in are small and consider a local form F^ = J d 4 xfk^(x)- A quark mass term is 
linear in x an d translates into a source term for $ 



Z 



^Tr(x j m + m f x) — ► -^±Tr(&m + m f $) = ~Tr(^j + /$) (D.6) 



g g 2 

where m = diag(m u ,m rf , . . .). A factorizing four fermion interaction yields 

m^Tr X t x — ?^Mt$ + m^T^X + *t$) . ( D . 7 ) 

The second term corresponds to the Yukawa interaction with /i fc<s> = m 2 k<j> g. We can therefore 
extract a q from eq. (D.6) as 

a q = 2fnl. (D.8) 

We note that only the terms linear and quadratic in x influence the value of a g . We could 
either restrict the composite fields from the beginning to the ones contained in the 0(4) 
symmetric linear a-model or work with all the fields contained in a complex 2x2 matrix $. 
In the latter case the anomaly term would contribute to both the masses and the Yukawa 
coupling. The net result is the same with m\ denoting the relevant mass term for the 0(4) 

53 The summation over internal indices as well as the integration over momenta has been suppressed. 
For complex Xabisi) similar terms have to be supplemented in the expansion. See ref. [44, 122] for a more 
detailed description. In our Euclidean conventions one has xlb ~ +9*' t l ; Rb 1 l J La- 
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D The quark mass term 



vector. For our conventions with h k<s> = 1 we have to normalize with g = m k 2 . Finally an 
eight fermion interaction becomes 

(Tr^x) 2 — \ A fe ,(Tr$t$) 2 + A fc ,Tr$t$ Tr(*t$ + $t % ) + . . . . (D . 9) 

We see here the appearance of terms quadratic in the quarks involving higher powers of 

There is an alternative, equivalent way of understanding the relation between j and 
m q : The quark masses in the picture with mesons must be equal at the scale kq> to the 
current quark mass m q (kq>). Let us consider an 0(4) symmetric fermionic interaction 
m^Trx^X + f^fcsC^X^x) 2 which leads to a meson potential 

U k , = m 2 ^ Tr $t$ + lx k , (Tr $t$) 2 . (D.IO) 

In the mesonic picture the quarks acquire masses through the Yukawa coupling to $ 

M k = p^(l + =^Tr($t )fe ($) } ($)fc (D.ll) 

where the second term arises from the higher order coupling in (D.9). Here ($) fc = 
diag(y? n , ipa, ■ ■ .) is the expectation value at the coarse graining scale k in the presence 
of the source term and M k = diag(M u , Ma, ■ ■ ■)■ It is sufficient to specify the dependence 
of U k<s> on real diagonal fields § qq . Then the ip q are determined from the condition 



dU, 



kit 



<9$ 



W = 2 [ml + A*. £ <4 V, = ^ • ( D - 12 ) 



9' 

Identifying M k=k<t in (D.ll) with m(k<f>) one has 

1 + =T£ = £ = 2< + 2 A fc . £ (D.13) 

and we recover (D.8) or, in our normalization with = 1, h k<6 = 1, 

j = 2m^m. (D.14) 

It is remarkable that higher order terms (e.g. ~ X k ^) do not influence the relation between 
j and m. Only the quadratic term rn\ enters which is in our scenario the only relevant 
coupling. This feature is an important ingredient for the predictive power of the model as 
far as the absolute size of the current quark mass is concerned. An appearance of higher 
order couplings in a q would make it very hard to compute this quantity. We emphasize that 
within our formalism there is no difference of principle between the current quark mass and 
the constituent quark mass. Whereas the current quark mass m q (k^) at the normalization 
scale /i = kq, corresponds to M q>k at the compositeness scale fc$ the constituent quark mass 
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is M 9j fe =0 . As k is lowered from k$> to zero one observes that the quark mass increases, 
similarly to the running current quark mass. Once chiral symmetry breaking sets in at the 
scale k x sB there is a large increase in the quark masses, especially for M u and 

The formalism of composite fields also provides the link [44] to the chiral condensate 
ipip^ since the expectation value ($) is related to the expectation value of a composite 
quark- ant iquark operator. For A = one has [122] 

+ ] )k = -4r (k<s>) + m q (k*) + mjM (D.15) 

with (tptp) {kq>) a suitably regularized operator normalized at ji — kq>. 
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